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ABSTRACT 


The solution to the radiative transfer equation for a 
non-gray absorbing atmosphere in local thermodynamic 
equilibrium has been programmed for an IBM 7094 com- 
puter to calculate theoretically the outgoing infrared radi- 
ance in five wave number spectral intervals. The spectral 
range covered is from 5 to 20 microns. The theoretical 
outgoing radiances may be utilized in several ways: (1) to 
assist in selection of the desired spectral characteristics 
for infrared satellite experiments, (2) to test the validity 
of balloon and satellite observations, and (3) to aid in the 
analysis of planetary infrared radiation data. 

Theoretical thermal emission spectra have been com- 
puted for several representative model atmospheres for 
Earth and Mars. These spectra illustrate that a consid- 
erable amount of information is available from radiation 
data of five wave number resolution. This information al- 
lows the recovery of the atmospheric and surface param- 
eters from measured thermal emission spectra. 
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THEORETICAL COMPUTATIONS OF THE OUTGOING INFRARED 
RADIANCE FROM A PLANETARY ATMOSPHERE 


by 

Virgil G. Kunde 
Goddard Space Flight Center 


I. INTRODUCTION 


The increasing spectral resolution and accuracy of satellite measurements of the outgoing 
radiation from a planetary atmosphere require a corresponding improvement in the theoretical 
determination of the outgoing radiation in order to properly interpret the measurements. Examples 
of the improved instrumentation are the Satellite Infrared Spectrometer (SIRS) (Hilleary, et al., 1966) 
and the Michelson type Infrared Interferometer Spectrometer (IRIS) (Hanel and Chaney, 1964, 1965). 
The primary purpose of this investigation was to develop a program for theoretical calculation of 
the outgoing spectral radiance as an aid in interpretation of the IRIS measurements. The IRIS 
measures the thermal emission spectra in the spectral range of 500(20^) to 2000(5^) wave numbers 
with a spectral resolution of 5 wave numbers. Thus, the theoretically determined spectral radiances 
must be of equal or greater spectral resolution and cover the same spectral range. The thermal 
emission spectrum of a planet, measured at the top of its atmosphere, depends on many atmospheric 
and surface parameters. The most important atmospheric parameters are the types of optically 
active gases present, abundance and distribution of these gases and the temperature profile, while 
the most important surface parameters are temperature, pressure, composition and structure. 

These atmospheric and surface parameters are recovered from the measured thermal emission 
spectra by using a synthetic approach. A secondary purpose of the program is to assist in select- 
ing the desired spectral characteristics for infrared satellite experiments and to test the validity 
of balloon and satellite observations. 

II. SUMMARY OF SPECTRAL LINE SHAPES 


A monochromatic absorption line may be broadened by one or more of many different line 
broadening mechanisms, both physical and geometrical. The types of mechanisms which are 
applicable at any given time depend on the physical and geometrical conditions under which the 
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absorption takes place. In the infrared, for conditions representative of planetary atmospheres, 
only two physical mechanisms are of importance: 

1. Collisional or pressure broadening, due to collisions with neighboring atoms or molecules, 
which is characterized by an absorption coefficient with a deep, narrow core and broad 
wings. 

2. Doppler broadening, due to random molecular thermal motions, which is characterized by 
a Gaussian shaped absorption coefficient. 

Pressure broadening produces a line profile which may be symmetric or asymmetric with respect 
to the unbroadened line; the Doppler line profile is symmetric. 

A. Collisional Broadening 

The absorption or emission line of an atom or molecule under the influence of collisions by 
neighboring particles can be characterized by certain properties: 

1. Line intensity, 

2. Frequency profile or line shape, 

3. Asymmetries in the line shape, 

4. Line shift. 

These properties relate to fundamental properties of the emitting or absorbing particle, and of the 
interaction between these particles and the perturbing particles. To determine the fundamental 
properties of the particle or of the interaction during collision from the above listed properties 
of the absorption or emission line knowledge is needed of the basic physical laws governing the 
absorption and emission of radiation in the presence of a perturbing force. An extensive literature 
on collisional line broadening exists as it represents an area of application for basic physics, 
primarily quantum mechanics, and also is a means of obtaining information on the energy levels 
and inter molecular forces of atomic or molecular systems undergoing collisions. 

It seems worthwhile and appropriate to briefly review the theoretical background concerning 
molecular line broadening. It is to be emphasized that the only purpose of this review is to place 
the various line shapes and the physical assumptions on which they are based in the proper per- 
spective, in a theoretical sense. The details of the mathematical treatment and experimental 
verification are found in more extensive reviews by Townes and Schawlow (1955), Ch T en and 
Takeo (1957), and Breene (1961). The work of Margenau and Watson (1936), Foley (1946), Bloom 
and Margenau (1953), Unsold (1955), Benedict, et al. (1956a), Breene (1957), and Traving (1960) 
has also been considered in this review. 

(1) Introduction 

Historically, the theory of collisional broadening has been approached from two viewpoints. 
These are the phase-shift approach, where the spectral lines of the emitting and absorbing particle 
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are broadened by collisions with other particles, and the statistical or static approach, where the 
line broadening is due to the electric field of the surrounding particles. These approaches 
are limiting cases of a more general theory. The phase shift theory is considered valid for low 
pressures and for frequencies near the line center (1-3 cm’ 1 - Spitzer (1940), Goody (1964a)) 
while the statistical theory has its application for high pressures and for frequencies in the wing 
of the line. As the pressure increases, the line broadening changes from collisional broadening 
to statistical broadening. Holstein (1950) has shown that the statistical theory is a limiting case 
of the collisional theory. 

The general theory of collisional broadening considers the intensity distribution of dipole 
radiation from a time-dependent charge distribution. The following development has been given 
by Woolley and Stibbs (1953). A convenient starting point is the classical expression for the 
instantaneous power radiated by an accelerating electron (non-relativistic) 


G 



(1) 


with the displacement x being considered complex. The total radiated power G T is the time 
integral 


g t = 



( 2 ) 


Equivalently, g t can be considered as the integral over frequency of the emitted spectral intensity 
distribution 


J o 


I e (go) d O) 


Equating Equations 2 and 3 and substituting in Equation 1 


f T e M , f 2 e 2 .. ... 

I I e (go) dco = I x x* 

Jo Jo 3 c 2 


dt > 


(3) 


(4) 


where ^ is angular frequency. 

The acceleration x of the radiating electron may be determined by solving the equation of 
motion for x(t) . For an undamped oscillator, corresponding to monochromatic radiation from a 
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sharp line, the wave train is infinite and of constant amplitude. The wave train of a damped oscil- 
lator produces a broadened line and consists of radiation of different frequencies and amplitudes. 

The amplitude of the wave train can be represented by the Fourier integral 


x(t) = 




a (co) e l0}t dco , 


(5) 


where a (") is the spectral amplitude of the emitted radiation. The corresponding Fourier integral 
for a (w)is 


r 

a (co) =— x (t) e" 

^ J-a, 

Substituting Equation 5 in Equation 4, we can write 


dt . 


( 6 ) 


f I e (co) d co= 2 e2 

f 

r 

f 

J 0 3 c 3 277 

J-co 

J-»- 03 

Jo 


f f f <a' 2 w 2 a(<u)a*(o>')e i(,i> '“ )t d£oda>'dt , 

J- co J- cc Jo 


(7) 


f 

J o 


I e (co) dco = 


4 e 2 


f 

J o 


co A a (co) a* (co) dco 


Equating the integrands of Equation 7, we can define l e (w) as 


4 P 2 qA 

V («) = , 3 a (co) a* (co) 

3 c- 3 


(8) 


Using Equation 6 for a(oi) ? the above equation for the intensity distribution can also be written as 
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i e (*o 


2 e 2 cJ* 
3 7TC* 


l 


x (t) e“ ia,t d t 


0 ) 


With the definition of the dipole moment for the radiating system as 


M(t) = e x (t) . 


( 10 ) 


Equation 9 can be rewritten as 


I e (OJ) 


2 co 4 
3 7TC 3 



jjl( t) e" ia,t dt 
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( 11 ) 


Either Equation 8, 9 or 11 represents the general formulation for the intensity distribution 
of a radiating system under the influence of collisional broadening. An equation analogous to 
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Equation 11 can be derived from quantum mechanics, using time dependent perturbation theory. 
Terms in the Hamiltonian due to collisions with neighboring particles and to incident electro- 
magnetic radiation on the particle are considered perturbations. References to the literature on 
the quantum mechanical derivation of line broadening are found in Bohm (1960). The general case, 
too complex to be of practical use, is solved only for special limiting cases. A summary of the 
reduction of the general formulation to its limiting cases of phase-shift and statistical broadening 
can be found in Bohm (1960). 

For very low pressures, damping by collisions is negligible and natural broadening oc- 
curs due to radiation damping. The spectral intensity distribution, corresponding to radia- 
tion damping, can be derived classically by solving the equation of motion for a radiation 
damped oscillator and substituting this expression into Equation 9. For conditions representative 
of planetary atmospheres, the natural line width is negligibly small compared to the Doppler and 
collision line half-widths. 

Two basic assumptions concerning the particle interaction are necessary to reduce Equation 11 
to either the phase shift or static approximation. They are (1) adiabatic collisions — no transitions 
are induced by the collision interaction, and (2) the classical path approximation. The adiabatic 
collision assumption requires that there be no energy loss from the emitted wave train due to the 
collision; the amplitude of the radiated wave must remain constant. Thus, the collision is con- 
sidered to affect only the phase of the radiated wave. Quantum mechanically, with the exclusion of 
effects of the radiation process, adiabaticity requires the state of the system to remain the same 
during the collision, i.e., the collisions cannot cause a shift in the energy levels. 

Criteria for determining if a collision is adiabatic or diabatic have been discussed by Breene 
(1961), Margenau and Lewis (1959), Spitzer (1940) and Townes and Schawlow (1955). In general, the 
adiabatic approximation starts breaking down in the infrared and microwave region as the transi- 
tion energies become about equal to the interaction energy. 

The classical path approximation consists of two assumptions: (1) the path of the perturbing 
particle can be described classically; (2) the path of the perturbing particle can be described by a 
straight line 


r * (t) = (b' 2 + v r 2 t 2) 1 / 2 , 


(12) 


where r ' ( t ) is the distance between the radiating particle and perturbing particle as a function of 
time, b' is the impact parameter and v r is the velocity of the perturbing particle relative to the 
emitting particle. This approximation is considered valid for encounters where the uncertainty 
in the velocity and position of the perturbing particle is small compared to the velocity and posi- 
tion of the perturber. From the uncertainty principle, Margenau and Lewis (1959) have shown that 
the classical path approximation is satisfied for large values of the angular momentum of the 
perturbing particle. 


5 



(2) Phase -Shift Theory 

The intensity distribution in the phase shift approximation can be obtained by considering the 
instantaneous frequency change Az7 in the emitted wave train, due to one collision, as 


m 

Av (t) =C/(r'(t)) n > 


(13) 


where C and n"' are constants depending on the properties of the radiating and perturbing particle 
(Ch’en and Takeo, 1957) and r' (t) is given by Equation 12. Values of n" for different interaction 
forces can be found in Townes and Schawlow (1955). The total phase change VC'O due to all colli- 
sions during the time interval 0-r is given by the integral 


t?' CO 


f 


A*/ (t) dt , 


(14) 


where Ao/(t) is the instantaneous phase shift due to all perturbing particles undergoing a collision 
at time t . If the duration of collision is considered to be short, only binary collisions are of im- 
portance and the integrand of Equation 14 can be considered as a summation of the phase shifts 
of each individual collision, 


A^'(t) (t) = £ C./(r! (t))”" , . (15) 

i i 

the subscript i denoting different values of the impact parameter. The exclusion of multiple col- 
lisions restricts the application to low pressures. The frequency of the emitted wave train can 
now be considered as 


°o =w o + 7 ?' 5 


(16) 


where co 0 is the natural angular frequency and the corresponding amplitude equation is 


x (t) = x 0 e 


(17) 


With the radiating system described by Equation 17, the intensity distribution (Equation 9) for the 
phase shift approximation, as given by Weisskopf (1933), is 


V (u) = 


2 e 2 ci/* 

3 77 c 3 


f 


iC^o + v (*) 


dt 


(18) 
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(2a) Impact theory: 

A limiting case of Equation 18 occurs for a type of interaction classified as a strong encounter. 
In a strong encounter, it is assumed that the emitted wave train is completely terminated by the 
collision and that no correlation exists between the phase of the wave train before and after the 
collision. Between collisions, radiation of the natural frequency of the oscillator is emitted. The 
main assumption which denotes an impact theory is that the time during collision can be neglected. 
This is equivalent to neglecting the phase-shifted radiation (V = 0) emitted during the collision. 
With the above assumptions, Equation 18 becomes 


I e (co) = 


2 e 2 co 4 

3 77C 3 


JT n /2 

f 

-Tp/2 


i(^0 -^)t 


dt 


(19) 


where T B is the time between collisions. The impact treatment involves a Fourier analysis of 
only the unperturbed radiation emitted between collisions. Equation 19 integrates to 


r (a,) 


2 e 2 
3ttc 3 



( 20 ) 


a line shape first obtained by Michelson (1895). The line broadening denoted by Equation 20 is due 
only to the finite extent of the wave train which is limited by the time between collisions. 

For random collisions, kinetic theory gives the probability distribution for the time between 
collisions as 


<D (v) = 1 /T Bo (v) 


-T b (v)/T b (V) 


( 21 ) 


where T B (v) is the mean time between collisions for a given molecular velocity v. Averaging the 
line shape of Equation 20 over the distribution function of Equation 21 yields the emission intensity 
per oscillator as 


1 __ 

(v 0 -T 7)2 + ( 1/2 77 T Bo )2 

The volume emission intensity is 


I e (V) = — (27r ^ X 2 1 

3 ttc 3 2 77 2 T 

B o 


( 22 ) 


i; O') = N' f ' I e (v) , 


(23) 
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where N' is the number of absorbing particles per unit volume and f ' is the oscillator strength, 
which represents the number of oscillators per absorbing particle. Equation 23 is normalized by 

i;<v)dv = i; , (24) 



and the volume emission intensity can now be written as 


I? 1/2ttT 0 (v) 


r (V) — 

77 (?* - V ) 2 + (l/2 77T B (v)) j 


(25) 


For local thermodynamic equilibrium, Van Vleck and Margenau (1949) have shown that the absorp- 
tion and emission line shapes are identical, thus we can write the mass absorption coefficient 
as 


k m (V) = S m /TT 


1/2ttT 0 (v) 

( (*-v o ) 2 + (l/2vT Bo (v)) 2 


(26) 


where S"" is the integrated absorption coefficient with dimensions cm 2 gm' 1 sec -1 . The above 
equation has been derived in the framework of the Weisskopf phase shift theory for impact colli- 
sions. However, Equation 26 is known as the Lorentz line shape as it was originally derived in a 
similar form by Lorentz (1906) by consideration of absorption of radiation by particles. The half- 
width of the spectral line a L (v), defined by k(z7) equals half its maximum value, is given by 


a L (v) = l/2 7rT Bo (v) 


(27) 


Goody (1964a) has pointed out that the absorption coefficient should be averaged over all molecular 
velocities, but that usually it is assumed that a mean half-width a L can be associated with the 
mean time between collisions. The expression from kinetic theory for the mean time between 
collisions is 


a L = 1/27T \ =Nv r A e > 


(28) 


where v r is the mean relative velocity of the colliding particles, N the number of colliding 
particles per unit volume, and A e is the effective collision cross-section. 

The mean relative velocity v r is 



( 29 ) 
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The line half-width can now be written as 


a 


L 



+ b i ) 2 



(30) 


where D is the optical collision diameter and M is the molecular mass. The i subscript repre- 
sents the type of perturbing particle. If the optical collision diameter is known, the line half-width 
can be calculated from Equation 30. In general, it is found that the optical collision diameter is 
larger than the gas kinetic collision diameter by a factor ranging from two to a hundred [Ch’en 
and Takeo (1957), White (1934)] . For a binary mixture of absorbing gas and nonabsorbing broad- 
ening gas, a L takes the form (Burch, et al., 1962) 

P 

e 



where the equivalent pressure P e is given by 


P e = P + (B - l) p a 


(32) 


The constant B is the self-broadening coefficient, a 0 is the half-width at reference pressure 
P 0 , P is the total pressure, and p a is the partial pressure of the absorbing gases. 

Weisskopf (1932) has arbitrarily defined an impact collision to have occurred when the phase 
shift V of the emitted wave is greater than unity. The emitted wave train is thus terminated at 
77 ' = 1. This approximation neglects two types of collisions: (1) close collisions with V > 1 , and 
(2) distance collisions with 77 ' < 1 . An optical collision diameter can be determined from Equation 
14 by equating 77 ' = 1 and solving for the corresponding impact parameter b^; the line half-width 
is then formulated in terms of as 


2>,<* [ 2k ‘ T (s*K)] , ! ' <33> 

Characteristic of the impact theory is the linear dependence of the half-width on the number 
of perturbers. As the phase-shifted portion of the emitted wave train is neglected during the time 
of collision, the impact theory does not directly predict a value for the line half-width in terms of 
the interaction forces. Only indirectly is this accomplished in the Weisskopf theory through the 
definition of a collision by 77' = 1. 
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The derivation of the Lorentz line shape using phase shift theory implicitly involves an as- 
sumption that a L and v - v 0 are smaller than v 0 (Goody, 1964a). Correcting for this assumption, 
the line shape is 


k m (v) 


S ra v 


|_( v - v o) 2 +a l O + v o ) + a lj 


(34) 


Further modification of the above theory is necessary to avoid contradicting a theory of nonresonant 
absorption by Debye (1929). If the orientation of the displacement and velocity vector of the oscil- 
lator after the collision is assumed to follow a Maxwell -Boltzmann distribution, the line shape 
becomes (Van Vleck and Weisskopf, 1945) 

a . cl 

7 r + 2 ’ ( 35 ) 

( V-V o) +°L ( V+V 0 ) +°L_ 

For a L « v 0 , which is satisfied in the infrared, the Van Vleck- Weisskopf line shape reduces to 
the Lorentz line shape. 

(2b) Line shift: 

Two of the assumptions involved in the impact theory were: (1) collisions due to distant 
perturbers are neglected as only collisions which cause quenching are considered, and (2) the 
time of collision is small compared to the time between collisions, which is the same as neglecting 
close collisions which give phase changes greater than unity. Consideration of these assumptions 
has indicated that the first is connected with line shifts and the second with line asymmetries. 
Lindholm (1941) considered the problem, including phase shifts due to distant collisions, and ob- 
tained the Lorentz type line profile 



k m (v) 


S m 

17 (^-^ 0 + ^) 2+a L 


(36) 


which is symmetric with respect to v 0 but with the line center shifted to the red of by the 
distance p. The line half- width is the same as the Weisskopf half-width (Equation 33) with the 
exception of a numerical constant. For van der Waals broadening, the relationship 




2.75 


(37) 


may be derived. The inclusion of the time during collision is in the realm of statistical broadening, 
which is discussed in the next section. 
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(3) Statistical Theory 

The statistical theory considers both the absorbing or emitting particles and the perturbing 
particles to be at rest. Through distortion of the energy levels, each spatial configuration 
produces a shift in frequency of the emitted radiation with respect to the frequency of the radia- 
tion which would be emitted if no perturbers were present. The line intensity in the angular 
frequency range dco is taken to be proportional to the probability of occurrence of the spatial 
configuration which yields a frequency shift in the range of dco. 

Considering van der Waals forces and several perturbing molecules, Margenau (1935) finds 
the intensity distribution to be 


l(v) 



A co < 0 


N 


Aco 3//2 


-in - n 2 
e 9 Aw 


A co > 0 


where the half-width is given by 


a = 0.822 tt 3 y N 2 


and the line shift (to the red) is 


A v = 



N 2 • 


(38) 


(39) 


(40) 


In the above equations y is a force constant, N the number of perturbers per unit volume, and 
A co = co ~ co 0 . For large Aoj, the red wing intensity approaches 


I (co) ~ AoT 3/2 . 


(41) 


The "statistical" theory in considering the time during collision yields an asymmetrical line 
shape. 


Work by Spitzer (1939a, 1939b, 1940) and Holstein (1950) has indicated the regions of appli- 
cability for the phase- shift theory and the statistical theory. Spitzer has shown that the phase - 
shift theory is valid near the center of the line for frequencies satisfying 


Aco << 


(2 77 C) 


1 / n 


(42) 
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Correspondingly, Holstein has shown the statistical theory is valid for frequencies satisfying 


A co >> - — - — _ 

t / m 

(2 77 C) /n “ 1 


( 43 ) 


which occurs in the wings of the line. 


(4) Lindholm r s General Theory - Line Shift and Asymmetry 

Lindholm’s general theory (Lindholm, 1945) takes into account the large phase shifts which 
occur during collision. Because the collision time is included, the theory is no longer considered 
an impact theory. Lindholm also included the small phase shifts, due to distant collisions, which 
lead to a line shift. The interaction force was assumed to be van der Waals (n = 6). The general 
formulation reduces to the statistical theory for high pressures and for the line wings and to the 
impact theory for low pressures and near the line center. As an approximation, Lindholm has 
determined the following expressions for the red and violet wings 


k m (v ) = .933 


[ m '(^ 0 ’ 3/2 ( red ) , 


M' 2 S m a 


k m (v) = .638 — [m'(^ - v 0 )] 7/ ^ 2 (violet) , 


(44) 

(45) 


where the notation of Plass and Warner (1952a, 1952b) has been used. The above equations are 
valid for \v - v 0 | > 3M* 1 . A value of M' of approximately 3 cm is typical for the earth T s atmos- 
phere (Plass and Warner, (1952a). From Lindholm T s general theory, Curtis and Goody (1954a) 
have derived wing equations similar to those of Equations 44 and 45 plus an equation representing 
the center of the line by assuming only that the line half-width is smaller than the frequency repre- 
senting the transition frequency between the impact and statistical theory. The effect of the 
Lindholm line shape on atmospheric transmission has been considered by Plass and Warner (1952) 
and Curtis and Goody (1954a; 1954b). The Lindholm theory has been experimentally verified for 
sodium lines broadened by argon (Kleman and Lindholm, 1945). The Lindholm theory is not con- 
sidered to be applicable for molecular lines in the infrared in general because adiabaticity is as- 
sumed, and in the particular case of n = 6, because close collisions cannot be described by a van 
der Waals force (Plass, 1954). 

(5) Anderson r s Diabatic Theory for Molecular Lines 

Anderson (1949) has generalized the theory of collisional broadening by including the effect 
of transitions between quantum states due to collisions. The starting point for Anderson is the 
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quantum-mechanical analog of the Fourier integral intensity expression given by Equation 11, with 
the subsequent treatment involving the assumptions of the classical path and of binary collisions. 
Anderson further assumes that the duration of the collision is small compared to the time between 
collisions, thus restricting the treatment to an impact type treatment. With the exception of a 
shift of the center frequency S A , a line shape similar to the Lorentz shape is obtained: 


k m ( V ) 


77 <?-*o - S a ) 2 +* 2 


(46) 


where 


a 


i^N,CDACO (I* A} 


1/2 


(47) 


The collision cross-section A(T) is given by 


A(T)= Z] M T)A J 2 (T) ’ 


(48) 


with p being the Boltzmann distribution for the perturbing particle of quantum state J 2 , and Aj 
being a partial collision cross-section for the transition between states i and f of the absorber, 
with the perturber in state J 2 . The partial collision cross-section is 


A (T) 
J 2 


■ r 

J o 


2tt b' S (b, J ) db' 


(49) 


with S (b ' , J 2 ) representing the probability that the perturber in state J 2 at impact parameter dis- 
tance b' will induce a transition from state i or f of the absorber. Following Ch’en and Takeo 
(1957): 


S(b'J 2 )=I 


L 


L f 2 




2 ,M 2 | P 2 | f.m f .J 2 ,M 2 > 


(2J f + 1) (2 J 2 + 1) 


L 


<i."> i .J 2 .M i |P 2 |i,m., J 2 ,M 2 ) 
(2 J, + 1) (2 J 2 + 1) 


(50) 


where J is the rotational quantum number, 2j f + 1 is the number of degenerate states, m f is the 
degenerate index for the state f , and M 2 is the magnetic quantum number of the perturbing particle. 
A typical matrix element is given by 

r 00 

(m' | P | n’) =1 (m' |Vj (t) | n') exp (i co m . n . t) dt , ( 51 ) 

J-a> 
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with V t (t) being the perturbing potential due to collisions, and the transition frequency 

between states m' and n\ Formulation of a for several types of molecular interactions is given 
by Townes and Schawlow (1955) and Tsao and Curnutte (1962). An important result of Anderson’s 
theory is the prediction of the line half-width and its dependence on the rotational quantum number 
and temperature. 

Benedict and Kaplan (1959) have applied the Anderson theory to calculate the half -widths of 
H 2 0 lines broadened by N 2 collisions. The main limitations of their calculations are: (1) the 
molecular interaction is due only to the dipole of H 2 0 and the quadrupole of n 2 and (2) all the 
half- widths depend on the calibration of one observed H 2 0 microwave line. Half- widths for the 
pure rotational H 2 0 spectrum were computed for several temperatures. In the Anderson theory 
the collision cross-section A(T) depends on temperature through the Boltzman distribution for the 
perturbing molecule and the temperature dependence of the partial collision cross-sections (Smith, 
Lackner, and Volkov, 1955). The partial collision cross-section temperature effect is due to the 
dependence of A T on the relative velocity between the absorbing and perturbing particles. In the 
range 220°-300°K, Benedict and Kaplan found the temperature variation of the half-width could be 
described for a constant pressure as 


a (T) = a (T 0 ) 



(52) 


where n" varied from 0.756 to -0.045. For 300°K, the line half-widths ranged from .111 to .032 
cm" 1 atm" 1 . For collisions where Aj 2 does not depend on temperature, n" is equal to 0.5. In 
addition, half-widths for lines in the v 2 and 2 v 2 vibration -rotation band were computed. The re- 
sults showed a change in half- width ranging from +2% to -18% with respect to the pure rotational 
half-widths. Later calculations by Benedict and Kaplan (1964) for self-broadened half-widths for 
the pure rotational H 2 0 lines gave half-widths which varied from .566 to 0.61 cm" 1 atm' 1 and a 
temperature exponent n" which varied from 1.04 to -.218. 

Comparison of Anderson’s theory with experimental values of the half- width, which is limited 
to a few cases, in general indicates fair agreement [Townes and Schalow (1955), Goody (1964a)]. 
However, Sanderson and Ginsburg (1963) have found experimentally, for three far infrared water 
vapor lines, self-broadened and nitrogen-broadened half-widths larger by a factor of 1.5 than the 
theoretical values of Benedict and Kaplan. 


(6) Benedict Modification of Lorentz Line Shape 


As the Lorentz line shape is theoretically valid only for frequencies near the line center, 
several experimental studies have been made to determine the shape in the wing of a collision- 
broadened line. For the HCI fundamental, Benedict, et al. (1956b) determined that for v - v 0 in 
the range of 1.5 to 40 cm' 1 the line wing could best be described by 


k (i/) % 


1 


o - V 0 T 


(53) 
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where m" is between 1.7 and 1.85, and v represents wave number. The observed wing gives greater 
absorption than does a Lorentz wing (m" = 2). 

From absorption by the fundamental and first overtone of CO, Benedict, et al. (1962) have 
determined experimentally that the Lorentz line shape should be modified by an exponential mul- 
tiplying factor dependent on \v - v 0 \ for \v - v 0 \ > dcm" 1 . The line shape is 


k ( v ) = — - 


(v- v 0 ) 2 + a l 


I V - v 0 I < d 


(54) 


k(v) 


-a # [|v- v 0 | - d] 


(V-V Q ) 2 + 0,2 


| v - v 0 I > d 


(55) 


where the units of k(v) and S are cm -1 atm" 1 and cm" 2 atm" 1 , respectively. The transformation 
from mass absorption coefficient to absorption coefficient per unit length per unit pressure k(v) 
will be discussed in section IIIA. For CO, the constants were found to be d = 4 cm" 1 , b" = 1 , and 
a ' ~ 0.015. The line shape of Equation 55 is denoted as the Benedict modification of the Lorentz 
line shape. Due to the exponential factor, the Benedict modification drops off much faster in the 
wings than does the pure Lorentz line. 

Bignell, et al. (1963) have observed less absorption than that given by the Lorentz line for the 
15^ band of C0 2 . 

Winters, et al. (1964) have made an analysis of absorption in the spectral region short of the 
4.17^ R-branch band head of the 4.26 C0 2 band. The absorption for wavelengths less than 4.17/i 
is due to the wings of strong lines in the R-branch. The line shapes determined from these measure- 
ments were of the same form as the Benedict modification of Equation 55. For self-broadened C0 2 
the appropriate constants are d = 5 cm" 1 , a 1 = 0.08, and b" = 0.8. In addition, measurements 
were made for carbon dioxide pressurized with nitrogen and also with oxygen. No detailed reduc- 
tion exists for the foreign gas broadening; however, Winters, et al. suggest that values in the 
range b" = 0.46 and a' = 0.46 would best fit the Benedict modification for atmospheric conditions. 

As the lowest C0 2 partial pressure of Winters, et al. was .25 atm, experimental evidence for con- 
ditions representative of the earth T s atmosphere is not available. 

Kyle, et al. (1965) have obtained numerous solar spectra in the 4.3^ C0 2 band at altitudes 
from 10 to 30 km along slant paths in the atmosphere. Comparison of theoretical slant path 
transmittances with experimental results gave satisfactory agreement at high altitudes for the 
Lorentz shape and at low altitudes for the Benedict modification. 

From absorption measurements on the high wave number side of the head of the 3 v 3 C0 2 band, 
Burch, et al. (1965, 1966) have determined that the observed absorption is less than that calculated 
for a Lorentz line. The authors found that the line shape varied with the type of broadening gas and 
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that the shape in this region (~7000 cm" 1 ) differed from that found by Winters, et al, in the region 
of 2400 cm" 1 . 

Bignell, et al. (1963) determined the Lorentz line shape giving the best agreement with H 2 0 
extinction measurements of the 8-12 p window region. These observations tested the line shape 
very far into the wings (~1 00-600 cm" 1 ) of the 6.3 /j. vibration -rotation band and the pure rotational 
lines of H 2 0. From sky emission measurements, Bolle (1965) has derived water vapor continuum 
absorption coefficients for window regions in the 7.5-26/x region. Comparison of theoretical con- 
tinuum absorption coefficients with the experimental values indicated a Lorentz line shape. 

The Lorentz line shape and the Benedict modification in the wings of a Lorentz line are the 
most applicable of the collisional line shapes in the infrared, with the Anderson theory being the 
most applicable method for calculating the line half-width. To date, no theoretical basis has been 
established for the Benedict modification. Further theoretical and experimental work is needed 
to obtain the applicable line shape for the wings of lines. 

B. Doppler Broadening 

The motion of an atom or molecule, through the Doppler effect, causes the frequency (V) of 
emitted radiation to shift from the monochromatic frequency (z7 Q ) of emission of a stationary 
atom or molecule. The derivation of the absorption coefficient for Doppler broadening assumes 
(1) the velocities of the emitting atoms or molecules in the line of sight have a Maxwellian 
distribution, and (2) the intensity of the spectral line at frequency v is proportional to the number 
of atoms or molecules radiating at frequency v. Using the Doppler principle, the absorption 
coefficient can then be derived (Aller, 1953) 


k m 




(56) 


where v p is the most probable Maxwellian speed for a given temperature. Substituting for the 
most probable speed 


V 


2 

P 



(57) 


the Doppler absorption coefficient in terms of frequency becomes 


k m (v) 


_ S m / 

= ^ V 


2 77k b T 


exp 


2 k b T 


j -v 2 


(58) 
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In wave number units, 


kM= ^/ 


2 77k b T 


exp 


2 / V - V. 


2k b T \ 


(59) 


The Doppler half-width, at the point where the absorption coefficient equals one-half its maximum 
value, is 


“d = ]/ 


2k K T 


ln2vl 


(60) 


which reduces to 


a D = 3.58x , 

where M is the molecular weight. The corresponding equations in wavelength units are 


(61) 


k(\) 


= x / 


2vk h T 


exp 


2k b T\ \ 0 


V 2 


a D = 3-58 x 10-/f\ 0 


^ • 


(62) 


(63) 


Figures 1 and 2 show the Doppler half-width graphed with T/M as a parameter in terms of wave 
number and wave length units, respectively. 



Figure 1 — The Doppler half-width in wave number units 
is given over a range of wave numbers for different 
values of T/M. Special cases of T/M are given for 
H 2 0 and C0 2 at temperatures of 300 and 3000 °K. 



Figure 2— The Doppler half-width in wavelength units 
is given over a range of wavelengths for different values 
of T/M. Special cases of T/M are given for H 2 0 and 
C0 2 at temperatures of 300 and 3000 °K. 
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C. Collisional-Doppler Broadening 


For physical conditions, as in the earth’s stratosphere, where both the collisions of neighbor- 
ing particles and the thermal motion of the emitter contribute significantly to the broadening of a 
line it is necessary to consider simultaneously the collisional and Doppler line shape. The absorp- 
tion coefficient for a mixed line can be written (Aller, 1953) 


k O) 


, -t ' 2 

k y_ r dt , 

n J-„ y ' 2 + (X' - t ') 2 


(64) 


where 


k 


o 




(65) 



V - V 

x' = ° /Wi . (67) 

a D 


In the above expression, the collisional line shape is assumed to be Lorentzian. The integration 
over t ' represents the integration over the different velocities of the absorbers in the line of 
sight, assuming the distribution of absorbers in terms of velocity is Maxwellian. The integration 
over velocity thus accounts for all of the radiation which is Doppler shifted to wave number v. 
Equation 64 can be rewritten as 


k O') =k 0 H (y\x') 


( 68 ) 


with 


H (y', x') =i_ 
n 


— — 

J.« y' 2 + (x' - t') 2 


dt' . 


(69) 


As H (y\ x') cannot be evaluated analytically, various methods of evaluation for different ranges of 
x' and y' have evolved (Penner, 1959); however, no method is suitable for the entire range of x' 
and y' occurring in a planetary atmosphere. Young (1965) has generated a method and developed 
it into a Fortran program for the evaluation of H(y', x'), covering a large range of x' and y', for 
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routine use on a digital computer. The evaluations of the mixed line shape considered in this paper 
have been made utilizing Young 1 s program. 


The quantity k 0)/S for a mixed line is graphed as a function of wave number for a range of 
pressures for a D = 0.00014 cm' 1 in Figure 3 and for a D = 0.016 cm' 1 in Figure 4. Also shown 
are the k(v)/S profiles for pure collisional and pure Doppler broadening. From Figures 3 and 
4 it is evident that Doppler broadening determines the shape of the core and collisional 
broadening determines the shape of the wings of the mixed line shape. The calculations on 
which Figures 3 and 4 are based indicate that the mixed line shape should be used for pressures 
smaller than 100 mb and for distances from the line center less than 2.5 cm" 1 . The exception to 
this is for very low pressures, where only pure Doppler broadening need be considered. Outside 
these ranges, only collisional broadening need be taken into account. The above limits of 

100 mb and 2.5 cm" 1 are valid provided the 



Figure 3— The ratio k(v)/S is given as a function of 
v -v Q for the mixed, the Lorentz and the Doppler line 
shapes for various values of a L / a D . The Doppler half- 
width is 1.432 x 10' 4 cm' 1 . 


Figure 4— The ratio k(v)/S is given as a function of 
v “ v 0 for the mixed, the Lorentz and the Doppler line 
shapes for various values of a \_/ a t> The Doppler half- 
width is 1.601 x 10~ 2 cm’ 1 . 


III. BASIC EQUATIONS 

(Aj Infrared Radiative Transfer Equation 

The outgoing spectral radiance from a planetary atmosphere is computed for a slant path 
through concentrically stratified atmospheric layers. As shown in Figure 5, the atmospheric 
path along s is specified by the geometrical parameter r 0 . The integration is from the top of 
the atmosphere to the level R s which represents the lower boundary of the path. For a clear 
atmosphere, R s denotes the planetary surface. Only the case where the atmospheric path intersects 
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Figure 5— Slant path geometry through the atmosphere. 


atmospheric slant path is shown in Figure 6. For 
dynamic equilibrium, the outgoing spectral radian 



Figure 6— Illustration of radiation components con- 
tributing to outgoing planetary radiance. 


the lower boundary surface is considered in 
Figure 5. A more detailed description of the 
a pure absorbing atmosphere in local thermo- 
:e I„ + (S T ) at level S T (Chandrasekhar, 1960) is 


K ( s t) 


- r s T k m 

ij (0) e lo 


( v >p P 


f 


B v ( 8') e 




(v) p ds 


( v )P t ds' 


(70) 


where k m (v) is the mass absorption coefficient, B v is the Planck function, and o f , is the density of 
the absorber. To conform to the system of units generally in use with respect to the absorption 
coefficient, it is necessary to transform the mass absorption coefficient k m (v) into an absorption 
coefficient per unit length per unit pressure k(i'). The transformation is made by defining 


k (?) = 


k™ ( v ) p r 
Pr 


(71) 


and 


du = - p r 



ds , 


(72) 


where p r is the partial pressure of the absorber, u is the optical path length in cm atm, and p r is 
the absorber density corresponding to reference conditions, usually taken at NTP. The absorp- 
tion coefficient k (v) has units of cm" 1 atm" 1 . With Equations 71 and 72 the product k m (^)p g ds 
may be written as 

k m (v) ds = - k (v) du . (73) 

The absorption coefficient transformation also requires that the integrated absorption coefficient 
be redefined. The transformed integrated absorption coefficient S in units of cm' 2 atm" 1 is 
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related to S m (Equation 26) 


s m 


P r c 


(74) 


The factor c enters the conversion because of the change from frequency to wave number units 
in the absorption coefficient. Making the substitution of Equation 73 into Equation 70, the radiative 
transfer equation becomes 


= 0) = I> s ) e 


r u s 

-J n k(V)du 
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B v (u') e 


-r 

J o 


k(v) du 


k(V) du' 


(75) 


where the coordinate convention is as indicated in Figure 6. The transmittance between the levels 
u' and u" in the atmosphere is given by 


T (u", u') = e 




k(^) du 


(76) 


where u" < u' . The radiative transfer equation can be written in terms of the spectral trans- 
mittance as 


K (T = 1) = I*(T®) T (0, u s ) 


f 


B„ (T„) dT„ (0, u') 


(77) 


or after an integration by parts as (Plass, 1956) 


I* (T = n = 1+ 05) T (0, u.) + B (T = 1) - B„ (T®) T (0, u ) 


J 

X(T„=1) 


(0, u') dB (T ) 


(78) 


Assuming the underlying surface does not radiate as a blackbody, the outgoing radiance of the 
surface consists of two terms 


K ( 1 ?) 


(79) 


where r® is the reflectivity for a perfectly diffuse reflector, and nF~ s is the downward monochro- 
matic radiant emittance. Because of the non-blackbody character of the surface, not all of the 
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radiation incident on the surface will be absorbed; some of the radiation will be reflected back up 
into the atmosphere. The reflected component of the outgoing radiance from the surface is given 
by the second term of Equation 79. The reflected component can be neglected for spectral regions 
where the downward atmospheric emission is small as in "window" regions, and where the surface 
emissivity is close to unity, as the surface reflectivity is then close to zero. 

Including the pressure and temperature dependence of the absorption coefficient, the trans- 
mittance (Equation 76) can be written as 


(P, T, 


U') rr e 


S'], k(P,T t V)du 


( 80 ) 


The evaluation of this equation requires an integration along the variable pressure and tempera- 
ture atmospheric slant path. Assuming an isothermal slant path ( T = T t ), an approximation by 
Curtis (1952) and Godson (1953) allows the transmittance over a variable pressure path to be de- 
termined by using a mean value for the pressure. With the mean pressure given by 



the transmittance simplifies to 


(81) 


T/P, T r u" 


-k(P, T v) (u' -u") 
u ) = e 1 


(82) 


The Curtis-Godson approximation thus equates the transmission along the variable pressure 
slant path with the transmission of a homogeneous path of length u' - u" at constant pressure 
p. In the weak and strong line limit, the Curtis-Godson approximation gives the same absorption 
as the exact expression. The Curtis-Godson approximation is more transparent at the line center, 
less transparent in the line wings, and gives a larger equivalent width than the corresponding exact 
expression (Goody, 1964a, Drayson, 1966). As one would expect, for a given atmospheric slab the 
error in the Curtis-Godson approximation reaches a maximum for intermediate values of Su. The 
maximum error also increases as the pressure difference over the atmospheric slab increases, 
becoming about 5% for end point pressures of 50 and 1000 mb (Drayson, 1966). In general, the con- 
clusion is that the Curtis-Godson approximation is fairly accurate for carbon dioxide and water 
vapor transmission studies, but is not very good for ozone. These results have been discussed by 
Kaplan (1959), Goody (1964a, 1964b), and Drayson (1966). 


The measured outgoing radiance for a spectral interval Av is 


[ Av < T V = 1) = f K < T V = »<**= f <fi v I ; (T l) T V (0, u s ) dv + f f 0 V B V (T V ) dT„ (0, u')dv, (83) 

•'Av -t® 
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with <P U describing the spectral response of the measuring instrument. We can define the average 
transmittance for Av 


T <h X 

I V V 

J A V 

J 


dv 


<t> v \ dv 


(84) 


where represents I„ + and B v of the first and second terms of Equation 83, respectively. Assum- 
ing l„ + and B v are constant over Av, 



Equation 83 can be written in terms of the average transmittance, as defined by Equation 84, 


lAv(T Av = 1) = T Av (0, 


u s^ 


I. 


K ( T Av> ^ + 


f f *,»„< 

d— s •'A v 


(T Ai )ch, dT Al/ (0, u'-) 


Av 


( 86 ) 


(B) Optical Path Length 

From Equation 72 the optical path length in cm atm for NTP conditions is 

du = - p g d S/ /p NTp . (87) 

The number of cm atm of the absorbing gas in a slant path can be found by taking the number of 
absorber molecules per unit area in the slant path, and determining the length per unit area which 
they occupy under normal temperature and pressure conditions. For C0 2 at NTP, one cm atm 
equals 0.00197 gm/cm 2 . In the atmosphere the density of the absorbing gas can be represented 
as 


p, = % p* . < 88 ) 

where q m represents the mass fraction of the absorber and p a , the total atmospheric density. With 
the ideal gas law and Equation 88 the optical path length becomes 


du = - 




(89) 
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where P and T represent the atmospheric pressure and temperature and M a and M g are the 
molecular weights of the atmosphere and the absorbing gas, respectively. Alternatively the 
above equation can be written in terms of the volume fraction of the absorber, q v = M a /M g q m , as 

du = - q v P/P ntp ds . (90) 

In the Earth's atmosphere, carbon dioxide has an average value — average with respect to 
season, latitude and altitude — of q m = .0477 gm/kg (Bolin and Keeling, 1963). For water vapor 
and ozone, q m is highly variable with respect to altitude. 

The water vapor content is usually specified in terms of the mixing ratio q^, which is the ratio 
of the density of water vapor to the density of dry air. The optical path length is considered in 
units of precipitable centimeters (pr cm) - the length of a column of liquid water formed by con- 
densing all the water vapor in a column along the slant path - the unit of conversion being 1 pr cm 
= 1245 cm atm. The number of pr cm is equal to the number of gm/cm 2 in the path. The optical 
path length for water vapor becomes 


du = 


P T ntp 


1245 Mjj P NTp T 


ds 


(91) 


From Equation 87, the optical path length for o in cm atm may be written as 


du = - q ds , 

3 


(92) 


where q Q is in cm atm/km. 

From the geometry of Figure 5, the relationship 


ds = r(r 2 - r 2 )" 1/2 dr 


(93) 


may be obtained with r 0 = sin 9 (R p + h), where R p is the radius of the planet. The integrated 
optical path length from the top of the atmosphere to the point r is 
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(96) 


u(cm atm, NTP) 
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ds (km) 


Assuming hydrostatic equilibrium for the atmosphere gives the pressure as a function of 
height, 


P(r) = P s exp 

where P s is the surface pressure, g is gravity, and R is the universal gas constant. 

The next section is devoted to a description of the methods used to determine the average 
transmittance over a finite spectral interval. In section V, the development of the computer logic 
involving the equations discussed in this section will be covered. 


/' 

•'R. 


g(r) M 0 10 s dr (km) 
RT 


(97) 


IV. AVERAGE TRANSMITTANCE FOR A FINITE SPECTRAL INTERVAL 

The average transmittance for a finite spectral interval, as defined by Equation 85, is 

T Av = 1/Av f T v dv , (98) 

where has been set equal to unity. The average transmittance may be determined semi- 
theoretically by using a band model, which can be evaluated analytically, to represent the vibration- 
rotation lines or by direct integration across the spectrum considering each vibration -rotation 
line individually. Which approach to use depends on the particular problem. Where high resolu- 
tion is not demanded, the band models offer a practical solution which can be computed very fast. 
However, as the width of the spectral interval decreases to the order of five wave numbers or less, 
the computational time required by a band model such as the quasi-random model begins to approach 
that of the direct spectral integration. Serious limitations of the band models are (Drayson, 1966): 

1. Representation of line strengths, positions, and collision broadened half-widths is not 
adequate. 

2. Of necessity, the Curtis-Godson approximation must be used for the vertical integration. 

3. Instrumental response functions and the mixed line shape cannot be simply taken into 
account. 


The utilization of the band model approach, in particular the quasi-random model, and the 
direct spectral integration approach will be discussed in more detail in the following sections. 
First, however, a brief digression will be made to consider the equivalent width of a single line. 
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A. Equivalent Width for a Single Line 

The equivalent width for a single line is 

+oo 

(1 - dv . (99) 



With 


k(i/) = Sb(V) , 


( 100 ) 


where b(v) is the relative line shape subject to the normalization condition 



b(v) Av = 


1 


the equivalent width is 


w = 



(1 _ e' Sub ( v )) Av . 


( 101 ) 


( 102 ) 


The equivalent width reduces to 


w = Su (103) 

for small path length and/or weak lines (Green and Wyatt, 1965). The above expression, known 
as the weak line approximation, is independent of the line shape and half-width. For large values 
of Su/a, the equivalent width for a Lorentz line is 


W = 


(104) 


which is referred to as the square root approximation. 

Equation 99 has been evaluated numerically for the Lorentz line shape and for the Benedict 
modified Lorentz line shape with the results shown in Figure 7. As one would anticipate, the 
exponential drop-off of the wings of the Benedict line shape causes the difference between the 
Lorentz and Benedict modification equivalent widths to increase as the pressure and absorber 
concentration increase, with the Benedict modification equivalent width always being smaller than 
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the corresponding Lorentz value. Equation 99 also has been evaluated numerically for the mixed 
line shape. The equivalent widths of a Lorentz line with = 1 x 10" 4 cm -1 are compared with 
the corresponding values for several mixed line shapes (Figure 8). The Doppler half-widths 



Figure 7— Curve of growth for Lorentz line shape and 
for Benedict modification of Lorentz line shape. 


for the mixed line shapes are assumed to be 
1.432 x 1(T 4 and 1.601 x 10' 2 cm" 1 . 



Figure 8— Curve of growth for Lorentz line shape and 
for mixed line shape. 


B. Quasi-Random Model for Molecular Line Absorption 

The various models (Elsasser, statistical, random Elsasser) that have been proposed to repre- 
sent molecular band absorption, approximations to these models, and their regions of validity have 
been discussed by Plass (1958, 1960). Of these models, only the statistical model is of interest in 
the present study as it is essential to the development of the quasi-random model. The quasi- 
random model will be presented in some detail as the author has programmed the quasi- random 
model for the IBM 7094, thus allowing a check on the accuracy of the transmittance calculations 
of Wyatt, Stull and Plass (1962b, 1963). 

(1) Theory and Computational Technique 

The statistical model (Goody, 1952; Plass, 1958; Howard, Burch and Williams, 1955; Goody, 
1964a) considers a spectral interval of width Av containing n total lines, with each individual line 
having intensity S i and wave number v. occurring within Av. Let n(v 2 ••• v. ••• v n ) dv x ••• dv. ••• dv n 
represent the probability that the i th line occurs between v. and v. + dv., and similarly P(S i )dS i 
is the probability that the i th line has an intensity between S. and S. + dS. . If v x ••• v n , Sj ••• S n 
are mutually independent, the probability of finding the set of n lines with the distribution v x ••• 

Si ». s n is 


22 N(^) PCS,) dv. d Si . (105) 
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The average transmittance at wave number v is found by averaging the transmittance at v over 
the probability distribution of the set of n lines 


/. /. /*00 /.CD n 

T(v)= L ... N(v.) P(S.) exp t-Si ub(v, v.)] dv. dS. . (106) 

*Av rAv *'o •'o 


Assuming that the line positions are distributed at random in Av and that all the lines have the same 
intensity, 


T(zA = 



~S u b(v , v. ) 


Av 


- dv . 


(107) 


The average transmittance T(v) is usually evaluated for the wave number v c occurring at the 
center of Av, with the assumption that T(v c ) is equivalent to the average transmittance t Av 


,-/■ 


T(V) dv 


(108) 


of the interval. The quasi-random model calculations of Wyatt, Stull and Plass (1962b, 1963) used 
the above assumption. The validity of this assumption will be discussed later. 

The quasi-random model differs from other band models in several ways (Wyatt, Stull and 
Plass, 1962a): 

a. Method of spectral division — To determine the average transmittance for the interval A*, 
a subdivision of A into smaller intervals 8 is made. The vibration -rotation lines occurring in an 
interval 8 are assumed to have a random distribution of line positions. The average transmittances 
for the 8 * s within A are arithmetically averaged to obtain the average transmittance of A. In 
principle, the exact line positions are obtained by allowing s to approach zero. 

b. Intensity distribution — The lines in a subinterval 8 are divided into subgroups, each sub- 
group covering one intensity decade. Wyatt, Stull and Plass have found that only the five strongest 
intensity decades need be retained in a subinterval. The intensity distribution is then simulated by 
the number of lines and the average intensity S for each decade. Considering the five intensity 
decades, Equation 106 gives the average transmittance at v due to the lines in s as 
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(109) 
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where the total number of lines in s . is 

j 

n. 

j 

c. Wing effects — The contribution to the transmittance of a given subinterval from the 
wings of lines in adjacent subintervals is taken into account. 

The total transmittance at wave number v in subinterval 8 k is 

CD 

= Z] ’ ( 11Q ) 

i=l 



where T. (v) is the contribution to T k (v) from the n. lines in interval 8. . The direct contri- 
bution to T k (v) occurs due to the n k lines in subinterval S k when j = k. The wing contribution 
is due to all remaining intervals for which j / k. Combining Equations 109 and 110, 
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-S . u b(v , v . ) 
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(111) 


From the above equation, the average transmittance for each subinterval 8 can be evaluated. 
Determination of the quasi-random average transmittance for interval A is accomplished by 
arithmetically averaging the T (v)’s of the subintervals contained in A. 

To avoid errors introduced by the arbitrary division of the spectrum into subintervals S, the 
quasi-random transmittance is calculated for an unshifted and shifted mesh. The two 
meshes are offset from each other by s/2 cm" 1 . For an interval 8 between v Q and v Q + S, the 
resultant transmittance at v c - + s/2 is then averaged from the two meshes 

T k (v c ) = 1/3 [f J (v 0 ) + T“ <v c ) + T* (u 0 + S)] , (112) 


where the superscripts s and u represent the shifted and unshifted meshes, respectively. 

The main term to be evaluated in Equation 111 is the expression for the average transmittance 
for a single line, 

rr 1/8 f e S i ub ^ ’ . (113) 

•'s. 

1 


29 



Introducing the following symbols 


y ~ v i - - s /2 , 

(114) 

z ~ v - v Q - S /2 , 

(115) 

"i = S. u/tto^ , 

(116) 

p ~ 2 a/S , 

(117) 

V = 2 y/S , 

(118) 

e = 2 z /8 , 

(119) 


and assuming the Lorentz line shape, 

T.(v) = 1/2 
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exp 


-p 2 ^ 

(e - n ) 2 + p 2 


drj . 


( 120 ) 


The interval S is considered to be between v 0 and v Q + 


In addition, is defined by 


^i 0 = s i p o /”*0<r > . 


( 121 ) 


where a 0 (T) is the half-width at pressure P 0 and temperature t, thus allowing Equation 116 to 
be rewritten as 


<fi = u/P . ( 122 ) 

X x 0 

The quantity and the number of lines n. have been determined semi-theoretically by Wyatt, 
Stull and Plass (1962b and 1963) for each of the five intensity decades for C0 2 and H 2 0 at 300°K, 
250°K and 200°K and listed in coefficient tables. The tabulations are every 2.5 cm” 1 for sub- 
intervals of 8 = 5 cm” 1 . The spectral range included is 505-10,000 cm” 1 for C0 2 , and 700-10,000 
cm” 1 for H 2 0. From Equation 120, Wyatt, Stull and Plass (1962a) have evaluated the direct contri- 
bution to the transmittance in analytical form as 


T.(v) = 1/2 (1 - e) O + 1/2 (1 + e) Q(i f., £') , 


(123) 
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where 


£ = p/ 1 - e 

(124) 

£' = p/1 + € 

(125) 

n(<fO = exp [-a 2 /(l + £ 2 )] -I^exp (-IfJ Jl 0 (l 


00 

+ 2 £ exp (- -i £■) 2^ X n (1^) sinn ^ + £f e xp (--ijjfj 

n ~ 1 

1 t MH * I .*,(|^)]^. U26) 

n- 1 


where \p = 2 tan” 1 £ , and I n is the Bessel function of imaginary argument and order n . Young 
(1964) has shown that the direct contribution may be evaluated numerically by dividing the integra- 
tion interval of Equation 120 into seven subintervals and applying seven-point Legendre- Gauss 
quadrature formula to each subinterval. The seven subintervals were chosen as 

j] = 0 . 0 , 0 . 001 , 0 . 005 , 0 . 01 , 0 . 05 , 0 . 1 , 0 . 5 , 1 . 0 . 

In this work, Young T s method was used with one exception, i.e., eight-point rather than seven-point 
quadrature formula was applied. The numerical and analytical integration agree exactly to six 
significant figures. 

The average transmittance T.(V) has been evaluated from Equation 120 for the center v - v c 
and for the lower boundary v - of the interval with the results given in Table 1. The large 
difference between T(^ c ) and TO 0 ) indicates that T(v) should be averaged with respect to 


Table 1 


Average Transmittance at Center and Boundary of Subinterval s . 


p'Z i 

p 

T(v c ) 

T(v 0 ) 

0.001 

0.05 

0.9724 

0.9860 

0.01 

0.05 

0.8454 

0.9202 

0.1 

0.05 

0.5413 

0.7464 

1.0 

0.05 

0.0897 

0.3544 

10.0 

0.05 

0.0 

0.1108 


frequency as discussed previously (Equation 108). In effect, this averaging with respect to v is 
accomplished in a fashion by Wyatt, Stull and Plass (1962a) by the averaging over the two meshes. 
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The wing contribution is obtained by assuming a 2 < < (y - v Q ) 2 . Equation 120 then becomes 
for the Lorentz line shape 


's'-if 

J -1 


exp 


- P 2 Z; 


which Wyatt, Stull, and Plass (1962) have evaluated as 


dn , 


(127) 


T.(vl =± «((e + 1) exp 


jfrd - (€ - ” " p [^rh]} -5”' /ia ' /! ["' (^) ■ *" (£9. 


, (128) 


where A = p 2 . For the Benedict modification, the wing expression is 




exp 


’-•f.p 2 exp [-a' (S/2) b "U - n| b " + a' | d | b "T 


di7 


(129) 


It should be noted that the above formulation of the Benedict modification used by Wyatt, Stull and 
Plass (1962a, 1962, 1963) is not the same as appears in Equation 55. The Wyatt, Stull and Plass 
formulation yields an absorption coefficient about 5-10% too large. In this investigation the Wyatt, 
Stull and Plass formulation has been used. The wing contributions of Wyatt, Stull and Plass 1962b, 
1963) were calculated using the Lorentz line shape for H 2 0 and the Benedict modification for C0 2 . 
As Equation 129 cannot be evaluated analytically Wyatt, Stull and Plass have assumed the trans- 



Figure 9— Wing transmittance due to lines occurring in 
adjacent subinterval (e = -2). 


mittance at the center of the interval (77 = 0) is 
representative of the average transmittance 
over the interval. The integration of Equation 
129 can then be avoided. In this investigation, 
8-point Legendre-Gauss quadrature formula 
was applied to the interval 77 = -lto ?7 = +l 
for Equations 127 and 129. This technique 
allows the average transmittance to be deter- 
mined to an accuracy of 5-6 significant figures. 
The wing transmittance due to lines in the 
nearest adjacent interval (e = -2) is shown in 
Figure 9 for both the Lorentz and Benedict 
modified line shapes. The average transmittanc< 
at the center of the interval for the Benedict 
modified line shape, as determined by the 8- 
point quadrature formula, is indicated in 
Figure 9 by the small dots. Comparison of 
the average transmittance with the transmittance 
for v = 0 gives an indication of the error intro- 
duced by Wyatt, Stull, and Plass in assuming 
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that the average transmittance can be approximated by the transmittance at 77 - 0. This compari- 
son is given in Table 2. The differences are significant enough to warrant inclusion of the 
numerical averaging. In their work, Wyatt, Stull and Plass used values of a' = 0.0675 and b" = 0.7 
for the Benedict modification parameters, which differ slightly from the values of a' =0.08 and 
b" = 0.08 given by Winters, et al. (1964). 


Average Wing 
in 

Table 2 

Transmittance T(v) due to Lines 
Adjacent Subinterval. 

pH, 

77 = 0 

Quad. Avg. 

.001 

.9998 

.9997 

.01 

.9977 

.9968 

.1 

.9772 

.9691 

1.0 

.7939 

.7455 

10.0 

.0995 

.1336 


The wing contribution for a given interval is due to wings of lines occurring in intervals both 
to the high wave number side (e = +2, 4, 6, • • • ) and the low wave number side (e = -2, -4, -6, . . .) 
of the given interval. For a gas at a specified temperature, pressure and optical path length, an 
interval at some maximum epsilon € M may be defined such that the contributions of lines in intervals 
with + e > + e M and with - e < - e M will have a negligible effect on the wing contribution of the given 
interval at e =0. Considering only the strongest line intensity decade ( i = 1) for 77 = 0, the wing 
transmittance may be written as 


T(v) - 



(130) 


Setting t(v) = 0.999, e M is defined by 


e 2 = ~ ni p2 ^ 1 

M In (0.999) ’ 


which can be simplified to 


e 2 = const n x 


fj p u . 


(131) 


(132) 
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Selecting the maximum combination of nj <f 0 j from the coefficient tables and inserting into Equa- 
tion 132 gives 


ejj = const Pu . (133) 

With P and u specified, e M can be determined* The wave number v H corresponding to e M can be 
determined using the relations in Equations 115-119. The parameter e M thus specifies the number 
of intervals on either side of a given interval that must be considered to insure that all lines 
which could possibly contribute. to the wing contribution of the given interval have been taken into 
account. The above expression is valid only for a Lorentz line shape; a similar expression can 
be derived for the Benedict modification. 

2 . Apparent Error in Wyatt , Stull and Plass Quasi-Random Transmittance Tables 

Using the numerical evaluation techniques discussed in the previous section, the quasi-random 
model was programmed on the IBM 7094 digital computer. Dr. Plass kindly provided magnetic 
tapes containing the coefficient tables of £ . and n i for H 2 0 at 300, 250 and 200°K and for C0 2 at 
300 and 200°K. As a side result, this investigation allows an independent check on the accuracy of 
the Wyatt, Stull and Plass (WSP) transmittance tables. In all cases, agreements with the WSP 
transmittance tables should be exact to at least three significant figures since the same input data 
and the same formulation is used. 

Exact agreement was found for H 2 0, with one exception. For values of Pu greater than ap- 
proximately unity, the transmittance values calculated in this investigation are lower than the 

corresponding WSP values, the disagreement 
increasing with increasing values of Pu . This 
is illustrated for the case v = 1002.5 cm" 1 , 

P = 1, T = 300°K in Figure 10. 

For C0 2 , disagreement was found for all 
the cases that were considered. In general, the 
disagreement correlates with the amount of wing 
contribution with the disagreement increasing 
as the wing contribution becomes a larger frac- 
tion of the total transmittance. In Figure 11 the 
transmittance for v - 605 cm" 1 , P - 1, and 
T = 300°K is graphed as a function of u. The 
transmittance values of this calculation are 
more transparent than the corresponding values 
of WSP, with the difference between the two 
values increasing with increasing u. A second 
comparison of C0 2 transmittances is illustrated 
in Figure 12 for the R branch lines of the 4.3 /x 
band for P = 1, T = 300°K and for several values 



Figure 10— Comparison of water vapor transmittances for 
V = 1002.5 cm- 1 and T = 300°K. 
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of u . From this figure the correlation of the difference between the two sets of transmittances 
with the amount of wing contribution is very evident. 



Figure 11— Comparison of carbon dioxide transmittances 
as a function of u for v - 605 cm" 1 ; P = 1 and T = 300° K. 

The reason for the disagreement for the 
H 2 0 and C0 2 calculations is not known at the 
present time. However, there seems to be a 
strong possibility that the C0 2 transmittances 
of WSP may be in error. Stull* has indicated 
that instead of using the wing contribution for a 
given interval, the wing contributions of the two 
adjacent intervals may have been used, thus 
overestimating the absorption. A few checks by 
error into the present calculations would give m 
WSP. 



Figure 12— Comparison of carbon dioxide transmittances 
as a function of wave number for various values of u 
with P= 1 and T = 300° K. 

the author have indicated that inducing this type of 
uch better, although not identical, agreement with 


3 . Comparison of Quasi-Random Transmittances to Experimental Data 

As the quasi-random transmittances are normalized to the low spectral resolution experi- 
mental data by the normalization constants C v v (T) (Wyatt, Stull and Plass, 1962b, 1963), one would 
expect that the quasi-random transmittances would in turn give good agreement with the experimental 
data. As a check, the quasi-random integrated absorptance has been calculated for several cases 
for H 2 0 and C0 2 . 

For the 6.3 /x H 2 0 band Burch, et al. 1962, have measured the integrated absorptance for the 
spectral intervals 1250 to 1590 cm" 1 and 1590 to 2100 cm” 1 for a range of pressures and H 2 0 con- 
centrations. The largest path length was 0.077 pr cm. A comparison of the calculated and 
measured integrated absorptance is given in Table 3. The fourth and seventh columns represent 


*V. R. Stull. Private communication, 1965. 
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Table 3 

Comparison of Quasi-Random and Experimental Integrated 
Absorptances for the 6.3/., H 2 0 Band.- 



J Ad v 1250-1590 cm" 1 

/ Adv 1590-2100 cm" 1 

Sample 

Pe 

u 

Exp. 

a L — 

a L = * 16 

Exp. 

a L = .100 

a L = .115 


atm 

pr cm 







i 

.0184 

.0017 

12.3 

9.6 

15.0 

12.7 

13.4 

15.4 

2 

.0232 

.0018 

14.9 

10.9 

17.1 

14.4 

15.4 

17.6 

3 

.0345 

.0019 

17.0 

13.4 

20.9 

16.4 

18.9 

21.6 

4 

.0530 

.0019 

21.1 

16.2 

25.0 

22.4 

22.9 

26.1 

5 

.1167 

.0017 

30.0 

21.3 

32.6 

30.3 

30.3 

34.5 

6 

.3540 

.0019 

50.9 

35.1 

51.8 

50.9 

50.3 

56.7 

7 

1.0171 

.0018 

79.4 

48.9 

69.0 

79.4 

71.2 

79.1 

8 

.0184 

.0034 

18.2 

13.7 

21.4 

13.6 

19.1 

21.9 

9 

1.0171 

.0036 

102 

70.1 

94.8 

114. 

101 

112 

10 

.0447 

.0041 

30.6 

22.3 

34.2 

31.3 

31.4 

35.7 

11 

.0580 

.0047 

38.0 

26.8 

40.6 

39.2 

37.7 

42.7 

12 

.0755 

.0048 

45.5 

30.3 

45.6 

48.8 

42.7 

48.3 

13 

.1382 

.0045 

56.7 

37.6 

55.7 

59.5 

53.4 

60.2 

14 

.0447 

.0081 

47.3 

31.5 

47.4 

53.1 

43.9 

49.7 

15 

.0580 

.0092 

54.0 

37.4 

55.6 

59.8 

52.2 

59.0 

16 

.0755 

.0095 

63.9 

42.3 

62.3 

69.9 

59.4 

66.9 

17 

.1382 

.0089 

76.6 

52.3 

75.3 

80.7 

74.0 

82.8 

18 

.0447 

.0161 

63.7 

43.9 

64.6 

70.5 

60.8 

68.4 

19 

.0580 

.0183 

72.1 

51.9 

75.2 

80.2 

72.0 

80.7 

20 

.0755 

.0190 

86.6 

58.7 

83.9 

96.5 

81.7 

91.3 

21 

.1382 

.0176 

100 

71.5 

99.4 

117. 

101. 

112. 

22 

.0447 

.0322 

78.9 

60.5 

86.4 

99.5 

83.1 

92.8 

23 

.0580 

.0364 

90.4 

70.8 

99.2 

108. 

97.6 

108 

24 

.0755 

.0378 

110. 

79.5 

110 

125. 

110. 

122 

25 

.1382 

.0352 

131. 

96.0 

128 

153. 

135 

148 

26 

.0487 

.0045 

33.8 

24.3 

37.1 

38.1 

34.1 

38.8 

27 

.0491 

.0089 

48.1 

34.3 

51.4 

55.3 

47.8 

54.1 

28 

.0496 

.0179 

69.7 

48.2 

70.4 

79.2 

66.7 

74.9 

29 

.0500 

.0359 

93.0 

66.4 

93.8 

109. 

91.3 

102 

30 

1.0658 

.0102 

153 

113 

141 

181 

164. 

176 

31 

1.0632 

.0198 

190 

144 

172 

226. 

211 

224 

32 

1.0605 

.0389 

220 

178 

204 

273 

263 

275 

33 

1.0592 

.0770 

256 

212 

234 

319 

315 

325 


36 


the experimental values of / Adz;, with the fifth and eighth columns representing the theoretical 
calculations for a L = 0.10, which is the value used by WSP. It is evident from the comparison 
that the theoretical values are systematically too small, especially for the 1250-1590 cm -1 
interval. Trial and error calculations with different values of a L have indicated that values of 
a L = 0.16 and 0.115 for the intervals 1250-1590 cm -1 and 1590-2100 cm -1 respectively give fair 
agreement with the experimental data. The above half-widths were used only for the direct con- 
tribution with the half -width for the wing contribution kept at 0.10 cm" 1 . The calculated 
integrated absorptances for the half-widths of .16 and .115 are given in columns six and nine of 
Table 3. The half-widths given should not be considered as representing the actual average half- 
width of the respective position of the absorption band, but rather as a parameter adjusted to 
normalize the theoretical and experimental transmittances. In addition to the above comparison, 
the integrated absorptance for the entire 6.3^ band was compared to experimental measurements 
of Burch, et. al. (1955) which covered a wider range of optical path length. Using the half-widths 
of 0.16 and 0.115, satisfactory agreement was found. 


Table 4 

Comparison of Quasi-Random and Experimental Integrated 
Absorptances for the 15/^ C0 2 Band. 


p 

e 

u 

(cm- atm) 


J Ad^ (cm 

- 1 ) 

(atm) 

Experimental 

Theoretical 

(Young) 

Theoretical 
(This Calculation) 

0.00051 

0.58 

3.37 

5.8 

2.4 

0.00137 

1.53 

8.86 

10.1 

6.3 

0.00496 

5.56 

21.7 

26.7 

20.1 

0.00827 

9.20 

32.6 

39.3 

30.4 

0.0205 

5.73 

34.6 

43.3 

34.4 

0.0397 

5.73 

43.8 

54.0 

42.8 

0.0837 

5.73 

54.7 

67.3 

53.6 

0.2052 

5.73 

69.5 

83.6 

68.1 

0.3999 

5.73 

81.9 

94.0 

79.2 

1.0089 


95.3 

103.9 

93.2 


Values of J Adv for the 15/x band of C0 2 are given Table 4. Column 3 represents the experi- 
mental measurements of Burch, et. al. (1962). The integrated absorptances of column 4 are theo- 
retical calculations of Young (1964) using the quasi-random model but different line strengths and 
column’5 represents the results of the present calculation for a half-width of 0.60 cm" 1 . For the 
range of conditions considered in Table 4, the agreement between the present results and the ex- 
perimental values can be considered satisfactory. 
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4. Modification for Doppler Broadening 

Including the effect of the mixed line shape, the average transmittance for a single line (Equa 
tion 113) becomes 


T. O) 


4 1 


-k„ H(y', x ')u 


dv. 


(134) 


The above expression is evaluated using 8 point Legendre -Gauss quadrature formula for each of 
seven subintervals in the same manner as for the direct contribution of the pure Lorentz line. As 
discussed in Section He, the mixed line shape is used for pressures less than 100 mb and for dis- 
tances from the line center less than 2.5 cm' 1 . A comparison of the average transmittance for a 
Lorentz line and for a mixed line is presented in Table 5. 


Table 5 

Quasi-Random Transmittance for Lorentz 
Line Shape and for Mixed Line Shape, T = 300°K. 





Transmittance 

v (cm 1 ) 

Pressure 

Optical Path Length 


1 

Lorentz 

Mixed | 




585 

100 mb 

100 cm-atm 

.7722 

.7722 

640 

100 mb 

100 cm-atm 

.0605 

.0605 

665 

100 mb 

100 cm-atm 

.0005 

.0005 

585 

10 mb 

10 cm-atm 

.9745 

.9738 

640 

10 mb 

10 cm-atm 

.7702 

.7681 

665 

10 mb 

10 cm-atm 

.3528 

.3501 

585 

10 mb 

100 cm-atm 

.9073 

.9059 

640 

10 mb 

100 cm-atm 

.4169 

.4146 

665 

10 mb 

100 cm-atm 

.0709 

.0706 

585 

1 mb 

10 cm-atm 

.9902 

.9864 

640 

1 mb 

10 cm-atm 

.9175 

.9043 

665 

1 mb 

10 cm-atm 

.6791 

.6578 

585 

.001 mb 

10,000 cm-atm 

.9894 

.9142 

640 

.001 mb 

10,000 cm-atm 

.9172 

.8522 

665 

*001 mb 

10,000 cm-atm 

.6789 

.6197 
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C. Direct Spectral integration 


Drayson (1966) has developed a computer program to determine atmospheric transmittances 
by direct spectral integration across the absorption band, using theoretically calculated line posi- 
tions and strengths. Other features of his program are: (1) the mixed line shape is used, and 
(2) the inclusion of the integration of the absorption coefficient along the atmospheric slant path, 
thus avoiding the use of the Curtis-Godson approximation. The initial calculations by Drayson 
have been made for the 15 /x C0 2 absorption band. 

Drayson has kindly provided the author with a copy of his program. With only slight modifi- 
cation the program has been made compatible with the IBM 7094. In the development of his pro- 
gram, Drayson has divided the 15 /x C0 2 absorption lines into 982 strong and 1008 weak lines. This 
portion of the program was changed so that all the lines were considered as strong lines. The re- 
sults obtained with this program will be covered in a later section. 


V. PROGRAMMING LOGIC 

The logical development of the computer program for determining the outgoing spectral 
radiance, using the quasi-random transmittance, will be discussed in this section. The program 
is modular in form, consisting of three basic sections: 

1. Atmos 1 — For an input model atmosphere, Atmos 1 calculates the optical path length 
and the average pressure from a given level in the atmosphere to the top of the atmosphere as 
a function of height for each absorbing gas considered. The input model atmosphere specifies 
the temperature, the water vapor distribution and the ozone distribution as a function of height. 

The atmospheric slant path is specified by 0 and h (Figure 5). Initially, Atmos 1 numerically 
integrates the hydrostatic equation for a zenith angle of 0° from the surface to the top of the 
atmosphere using 0.1 km intervals. As the integration is carried out, the quantity ds (Equation 
93) plus the temperature and pressure at the mid-point of each 0.1 km interval is computed 
and stored. Using the above stored quantities along with the vertical gas distributions, the 
optical path length (Equations 94, 95, 96) is integrated from the top of the atmosphere to 
the lower boundary surface. At the same time the average effective pressure (Equation 81) 
is computed as a function of height. A matrix of height vs optical path length vs average 
effective pressure is written on magnetic tape for each absorbing gas, considered for a number 
of pre-selected heights, to serve as an input data tape for the second section of the program. The 
number of height levels in the matrix is arbitrary, subject only to the restriction that the product of 
the number of spectral intervals times the number of atmospheric levels cannot exceed 13000 for 
H 2 0 or 15000 for C0 2 . 

2. Atmos 2 — Using the Curtis-Godson approximation, Atmos 2 calculates the quasi- random 
transmittance for each spectral interval of interest for each atmospheric level of the input height 
matrix produced by Atmos 1. The widths of the spectral intervals were chosen as A = S = 5 cm -1 , 
which means that the average transmittance for A is based on the assumption of random line 



positions in A. In this sense the quasi-random model is not really being used, as the quasi-random 
model differs from a pure random model in the assumption of random line positions only for A > S . 
The width of A = 5 cm" 1 was chosen so that the coefficient tables of Wyatt, Stull, and Plass (1962b, 
1963) could be used. In principle, of course, A can be chosen to be any value. Atmos 2 requires 
the center wave number of the initial v x and final v F spectral interval to be specified. 

The transmittances are evaluated, considering the integrated absorption coefficients to. be at 
an isothermal temperature along the atmospheric slant path. The isothermal temperature which 
best represents the temperature variation along the slant path must be chosen from 200 and 300°K 
for C0 2 and from 200, 250 and 300°K for H 2 0. With the optical path length, the isothermal tempera- 
ture and the average effective pressure known, the quasi-random transmittance can be evaluated 
using the formulation given in section IVB1. 

The quasi-random transmittances for a given temperature, pressure, and optical path 
length are calculated as follows. To insure that the contributions of all lines are taken 
into account, all intervals from v 1 - v u to v F + v u must be considered. For the above 
intervals, the coefficients <f 0 and r}. corresponding to temperature T t for both the un- 

i 

shifted and shifted meshes, are read from magnetic tape and put in core storage. The direct con- 
tribution is then computed for each 5 cm" 1 interval from to v Y using the numerical quadrature 
of Equation 120. Next, the wing contribution is computed starting with the interval at v x - v M and 
subsequently considering all intervals until the interval v 1 + is reached. For a given interval, 
the associated coefficients were used to calculate the wing contribution to the transmittances at the 
center of all the other intervals from v r to v F . The wing contribution was computed for higher 
and lower wave numbers until the contribution exceeded 0.999. In each interval of interest, the 
wing contribution was multiplied into the values previously calculated. When the computation of 
the wing contribution was completed, the direct and wing contribution for each interval of interest 
were multiplied together. The same procedure was then used for the shifted mesh. The quasi- 
random transmittance is then obtained by averaging the transmittances of the unshifted and shifted 
meshes according to Equation 112. The transmittances are computed for each wave number range 
of interest starting with the level in the matrix representing the ground, and subsequently evaluating 
the transmittances for each level in the matrix until the top of the atmosphere is reached. The out- 
put of Atmos 2 consists of a height vs transmittance matrix for each 5 cm -1 interval. This data 
is put on magnetic tape to serve as input data for the third section of the program. 

3. Atmos 3 — Solves the radiative transfer equation for each 5 cm -1 interval. Atmos 3 reads 
the height vs transmittance matrix for a given spectral interval from the input magnetic tape. The 
outgoing radiance for this interval is then determined by integrating the radiative transfer equation, 
in the form of Equation 77, from the top of the atmosphere to the lower boundary surface. The 
integration is done numerically for atmospheric layers of 0.1 km thickness. The transmittance for 
each atmospheric level is determined by interpolation from the height vs transmittance matrix. 

The spectral response and spectral surface emittance must be specified as input parameters for 
Atmos 3. If more than one absorbing gas must be taken into account, the total transmittance is 
taken as the product of the transmittances for each individual gas. 
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Each of the above sections represents a separate Fortran IV program. The three sections 
cannot be combined into a single program as Atmos 2 utilizes all the available core storage of the 
IBM 7094. An IBM 7094 listing for the Atmos 1, Atmos 2 and Atmos 3 programs may be found in 
Appendix C, Appendix D and Appendix E, respectively. 

The description given above has been for transmittances as determined from the quasi-random 
model. However, the modular form allows one section to be modified or even completely replaced 
without affecting the remaining sections. This is extremely beneficial as Atmos 1 and Atmos 3 can 
be considered relatively permanent programs whereas Atmos 2 will undergo frequent changes as 
improved techniques for computing atmospheric transmission become available and are imple- 
mented. For example, the atmospheric transmission of Drayson, discussed in Section IVC, repre- 
sents a considerable improvement over the quasi-random band model and the Curtis-Godson ap- 
proximation. For the 15^ C0 2 band, this program essentially replaces Atmos 1 and Atmos 2 and 
it has been modified to write on magnetic tape the height vs transmittance matrix required by 
Atmos 3. 

VI. OUTGOING INFRARED RADIANCES FOR THE EARTH AND MARS 

The outgoing infrared radiances have been computed at 5 cm" 1 intervals for several typical 
model atmospheres of the Earth and Mars in the spectral range 500 to 2000 crti" 1 . The transmittances 
for this spectral range have been assembled from a variety of sources. The high spectral resolu- 
tion transmittances discussed in the previous sections were used in the manner described in 
Table 6. 


Table 6 

Sources of High Spectral Resolution Transmittances. 


Gas 

Band 

Spectral Range 

Source 

co 2 

15 m 

500 - 860 cm” 1 

Drayson (1966), direct sp. integ. 

co 2 

5.2, 9.4, 10.4/* 

865 - 2000 cm" 1 

Plass (1963), quasi-random model 

h 2 o 

6.3/u, 

1000 - 2000 cm" 1 

Plass (1962b), quasi-random model 


As a temporary measure, the transmission of the 9.1 and 9.6^ bands of ozone and the rotational 
lines of water vapor has been included at a spectral resolution much lower than 5 cm -1 . As the 
higher spectral resolution data become available for these spectral regions, they will be incorporated 
into the radiance calculations. These transmittances were determined using generalized absorp- 
tion coefficients. For the 9.1 and 9.6 y. ozone bands and the water vapor rotational lines, the gener- 
alized absorption coefficients of Elsasser (1960) and Moller and Raschke (1963), respectively, were 
used. 
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A. Earth 


The model atmospheres chosen to represent the Earth’s atmosphere in this investigation have 
been described previously by Hanel, Bandeen and Conrath (1963). The models were modified 
in several ways: (1) the water vapor distributions were assumed to be constant in height above 
the tropopause at the value occurring at the tropopause, (2) the 1962 U. S. Standard temperature 
profile was substituted for the 1959 ARDC Standard temperature profile, and (3) the 1959 ARDC 
Standard water vapor distribution was scaled so that the total amount of H 2 0 in the vertical was 
approximately 1 pr cm. Several characteristic properties of the models are listed in Table 7. 


Table 7 

Surface Temperatures and Total Gas Content for One Air Mass. 


Atmosphere 

Surface 

Water vapor 

Carbon dioxide 

Ozone 

Temperature 

u(pr cm) 

u(cm-atm) 

u(cm-atm) 

Tropical 

305.30°K 

4.63 

250.11 

0.24 

1962 U.S. Std 

288.15°K 

0.97 

250.04 

0.36 

Arctic Winter 

246.00°K 

0.14 

249.93 

0.34 


Hovis (1965, 1966a, 1966b, 1966c*) has made total reflectance measurements of some common 
surface minerals, iron oxides, and different types of igneous rocks in the wavelength interval 0.5 to 
22a. Emissivities derived from these measurements are shown in Figure 13 for selected materials 



Figure 13— Spectral emissivity of surface materials 
(Hovis, 1965, 1966a, 1966b, 1966c*). 


*W. A. Hovis, Jr. Private communication, 1966. 


from each of the above categories. The solid 
curve represents soil from the Mojave Desert of 
California; the dashed curve, limonite of particle 
size less than 0.105 mm from Caracas, Venezuela; 
the dotted curve, granite of particle size less than 
0.038 mm; and the dot-dash curve, sand from 
Daytona Beach, Florida, with particle sizes in the 
range 0.25 to 0.50 mm. The most pronounced 
infrared spectral features are the residual ray 
reflections in the 8 to 10 /x and 18to 23^ region. 

Less pronounced is the water of hydration feature 
at 6.3/x. The features for wavelengths less than 
5^ are not relevant to this investigation. The 
measurements of Hovis indicate that the emis- 
sivity depends strongly on the particle size of 
the sample. 
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The computed theoretical thermal emission spectra for the 1962 U. S. Standard, Tropical and 
Arctic Winter model atmospheres are shown in Figures 14, 15, and 16 respectively for the case of 
vertical viewing. The vertical scale for all three figures is brightness temperature, the black- 
body temperature which is necessary to produce the outgoing radiance in a 5 cm" 1 interval. The 
surface emissivity for the three models was assumed to be blackbody. In addition, the thermal 
emission spectra were computed for the 1962 U. S. Standard model considering the surface emis- 
sivities of Mojave Desert soil and of Daytona Beach sand. These spectra are also presented in 
Figure 14. No instrumental effects have been considered in these computations. 



Figure 14— Theoretical thermal emission spectra, 1962 U. S. Standard Model Atmosphere. The effect 
of surface emissivity on the spectrum is illustrated. 
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Figure 15 — Theoretical thermal emission spectrum. Tropical Model Atmosphere. 
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Figure 16 —'Theoretical thermal emission spectrum, Arctic Winter Model Atmosphere. 
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In Figures 14, 15, and 16 the most pronounced spectral features are the 6.3m H 2 0, 9.6m 0 3 , 
and 15.0 m C0 2 absorption bands and the H 2 0 rotational lines which dominate most of the spectral 
range from 5 to 20m- Minor absorption due to the 5.2 m C0 2 , 9.1m 0 3 , 9.4m C0 2 , and 10.4m C0 2 ab- 
sorption bands is also evident in the spectra. The minor constituents CH 4 , with an absorption band 
at 7.7m, and N 2 0, with absorption bands at 7.8m, 8.5m and 17.0m, have not been included in the present 
calculations. The ground transmission spectra, corresponding to the three model atmospheres, 
are given in Figure 17. The transmittances in the 9 to ^region are somewhat inaccurate as the 
"continuum” H 2 0 absorption has not been considered and because the generalized absorption co- 
efficients used for 0 3 overestimate the absorption of the 9.6m 0 3 band (Elsasser, 1960). 



Figure 17— Ground transmission spectra of model atmospheres. 


Considerable spectral structure can be seen in the 6.3m H 2 0 and 15.0m absorption bands. Simi- 
lar structure is not evident in the 9.6m 0 3 band and the H 2 0 rotational structure because the low 
spectral resolution absorption data were used. The outgoing radiance was computed for 5 cm’ 1 
spectral intervals with the center wave number of each interval being 5 cm’ 1 apart (500 cm’ 1 , 

505 cm’ 1 , 510 cm" 1 ,* • •)• Because of this procedure some spectral resolution is lost when the 
strongest absorption features do not coincide with the center wave number of an interval. For 
example, this has occurred for the 618 cm’ 1 C0 2 band in the 15.0/^ region. The 618 cm’ 1 and the 
721 cm" 1 C0 2 bands are approximately equal in strength yet do not appear so in the spectra. The 
15^ C0 2 band is shown on an expanded wave number scale in Figure 18. In this presentation the 
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Figure 18— Theoretical thermal emission spectrum of 15 m C0 2 absorption band, 1962 U. S. Standard 
Model Atmosphere. The brightness temperature has been computed every 1 cm* 1 . 
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outgoing radiance was computed for 5 cm' 1 intervals every 1 cm' 1 . The 15ft C0 2 band is due 
mainly to the v 2 fundamental and about fourteen overtone and combination bands of C 12 O 16 O 16 . 
The lower and upper state and transition wavelength (Young, 1964) is denoted for each of these 
bands in Figure 18. 

As mentioned previously, the effect of surface emissivity is illustrated in Figure 14. The 
solid curve represents a blackbody surface; the dashed curve, Mojave Desert soil; and the dotted 
curve, Daytona Beach sand. The major effect is in the 8 to 9m "window” where the residual ray re- 
flection has decreased the brightness temperature by about 4°K for the Mojave Desert soil surface 
and by 4 to 14°K for the Daytona Beach sand surface. In the 11 m "window" the brightness tempera- 
ture has been decreased by about 2°K. 

B. Mars 

Two Mars model atmospheres with different surface pressures (5 and 15 mb) were considered. 
The models include about 45 m-atm of C0 2 (Cann, Davies, Greenspan and Owen, 1965), 14 m of pre- 
cipitable water (Kaplan, 1964) and a surface temperature of 230°K. Current estimates of the total 
amount of C0 2 are somewhat higher. From the + 2v 2 + 3v 3 C0 2 band at 1.05m in the Martian 
spectrum, Belton and Hunten (1966) have derived a C0 2 abundance of 60 ± 26 m-atm while Gray 
(1966), from a personal communication by Schorn, has given 60 to 80 m-atm as the best current 
estimate of the Martian C0 2 abundance. In these calculations the value of 45 m-atm has 
been used. From near infrared spectroscopic observations Kuiper (1952) has set an upper 
limit of 0.05 cm-atm for 0 3 in the Martian atmosphere. Prabhakara and Hogan (1965) have 
used photochemical theory to determine the total amount of Martian 0 3 as approximately 0.03 
cm atm. In addition, Prabhakara and Hogan determined that this amount of 0 3 was small 
enough so that absorption by the 9.6m 0 3 band need not be considered. From recent ultra- 
violet spectral measurements of Mars from a rocket, Evans (1965) has determined an upper 
limit of 0.004 cm atm for 0 3 (Rea, 1965). In these calculations the effects of 0 3 have been 
neglected. The temperature profiles, based on the above parameters, were calculated* under 
the assumption of radiative and convective equilibrium and are shown in Figure 19. 

The computed spectra are presented in Figure 20 for a blackbody surface with surface pres- 
sures of 5 and 20 mb. With the exception of the 15 m C0 2 band, the Martian spectrum is essentially 
featureless in the 5-20m region. Because of low H^O abundance, the 6.3m H 2 0 band shows only 
weak absorption. The interval from 875 to 2000 cm -1 is shown in Figure 21 with the temperature 
scale expanded. The parallel structure of the upper state bands of C0 2 at 9.4 and 10.4m and the 
perpendicular structure of the 3 v 3 band at 5.2 f± can be clearly seen. Under the assumed condi- 
tions, water vapor shows a maximum absorption of about ten per cent or in terms of brightness 
temperature not exceeding several degrees. Consequently it will be difficult to identify H 2 0 in 


*R. A. Hanel and F. Bartco. Private communication, 1965. 
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Figure 1 9— Temperature profiles for Martian atmosphere 
(Hanel and Bartko, 1965). 


the Martian spectrum and even more difficult 
to determine its abundance. Fortunately the 
band spreads over 600 to 700 cm 1 and shows con- 
siderable rotational structure. This allows the 
use of statistical methods to obtain information 
about H 2 0 from the observed spectrum. Atmos- 
pheric H 2 0 has been identified using a cross- 
correlation technique* in which the H z O ground 
transmission is cross-correlated with the com- 
puted brightness temperature spectrum. A 
power spectrum analysis* of the ground trans- 
mittance of H 2 0, corresponding to the spectra 
of Figure 21, indicates a periodicity of about 
17 cm" 1 in the rotational structure. More will 
be discussed later concerning the 17 cm -1 peri- 
odicity. The 15 fj. C0 2 band is shown on an ex- 
panded wave number scale in Figure 22 for the 
Mars model atmosphere with a surface pressure 
of 5 mb. 

On comparison, several distinct differences 
can be discerned in the Earth (Figures 14 and 18) 
and Mars (Figures 20 and 22) spectra. The most 
significant difference is due to the low abundance 
of H 2 0 and lack of 0 3 on Mars. A second dif- 
ference can be seen in the Q branch of the 15 /x 
CO 2 band. The strongly absorbing Q branch 
absorbs fairly high in the atmosphere and thus, 
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Figure 20— Theoretical thermal emission spectra for Mars. The effect of surface 

pressure is illustrated. 


*B. J. Coarath. Private communication, 1967. 
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Figure 21— A portion of the spectra (700 - 2000 cm -1 ) shown in Figure 20 with an expanded 
brightness temperature scale to better exhibit the water vapor absorption features. 
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Figure 22— Theoretical thermal emission spectrum of 15 >la C0 2 absorption band for Mars model with 
surface pressure of 5 mb. The brightness temperature has been computed every 1 cm -1 . 

for the Earth, most of the absorption and subsequent re-emission at the Q branch occurs in the 
stratosphere where the temperature is increasing. The re-emission at higher temperatures leads 
to an increasing brightness temperature. For the Mars model, "seeing” higher in the atmosphere 
corresponds to cooler temperatures, and the brightness temperature decreases. Thus qualitatively 
the general shape of the temperature profile in the region of Q branch absorption can be determined 
by observing the shape of the Q branch in the spectrum. 

Due to the high concentration of C0 2 , the absorption by the C0 2 bands in the 5 to 20 m region is 
pressure dependent. Thus the procedure of using the two limiting absorption laws (weak band 
strong band) for determining both the atmospheric C0 2 abundance and surface pressure cannot be 
made in this spectral region. However, several techniques utilizing curves of growth analyses are 
available to determine both surface pressure and atmospheric C0 2 abundance. These include the 
measurement of the equivalent widths of two different C0 2 bands (e.g., the 5.2, 9.4, 10.4 and 15^ 
bands) or of resolvable spectral intervals within a given C0 2 band (e.g., the 15^ band). 

For the Mars model of 5 mb surface pressure, the effect of surface emissivity on the thermal 
emission spectra has been considered. The computed spectra are shown in Figure 23 with the 
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Figure 23~-Theoreticai thermal emission spectra for Mars model with surface pressure of 5 mb. The 

effect of surface emissivity is illustrated. 


solid curve representing the blackbody surface; the dashed curve, limonite; and the dotted curve, 
granite. The limonite and granite samples were chosen mainly to study the effects of the strong 
water of hydration feature at 6.3/x on the thermal emission spectrum. 

Whether limonite and/or granite is representative of certain areas of the Martian surface is 
unsettled. The main evidence for believing that the Martian desert areas are composed of 
limonite is: 

(1) Color — Sharonov (1961) has found that limonite best matches both the brightness and 
color index of the Martian bright areas. 

(2) Monochromatic albedo — Sagan, Phaneuf and Ihnat (1965) have compared Russell-Bond 
monochromatic albedos for Mars with laboratory samples in the spectral region 0.3 to 4 m. 
The best correlation was obtained with a pulverized limonite. 

(3) Polarimetric — Dollfus (1961) has attempted to match a curve of polarization versus 
phase angle, obtained from measurements of Martian desert areas in the visible part 
of the spectrum, with laboratory samples. The best comparison was obtained with the 
limonite sample. The reduction of the observational data assumed a surface pressure 
of 90 mb. The recent measurements of surface pressure by Kaplan, Munch and Spinrad 
(1964), Evans (1965) and Kliore, Cain, Levy, Eshelman, Fjeldbo and Drake (1965) have 
indicated a surface pressure less than 25 mb. Considering the change in surface pres- 
sure alters the polarization curve so that surface materials other than limonite will 
give the observed polarization (Younkin, 1966). 

Arguments against the existence of large amounts of limonite on the Martian surface are: 

(1) Mean density — Van Tassel and Salisbury (1964) reason that it is unlikely that the surface 
composition of a terrestrial planet would be mainly hydrated iron oxides, and that the low 
mean density of Mars compared to that of Earth indicates that the metal content of Mars 
is low. 
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(2) Near-infrared spectral reflectance — Laboratory reflectance data of powdered mixtures 
of geothite and hematite have been obtained by Draper, Adamcik and Gibson (1965) which 
indicate an absorption feature around 0.9/x with an increasing reflectance toward longer 
wavelengths. From reflectance measurements of the complete Martian disk, Younkin 
(1966) finds no indication of this absorption feature or of the increasing reflectance. 

(3) Limb darkening at 4920 and 6430 A — Coulson (1965) has compared the limb darkening of 
a planet with a limonite surface with Mars observations at 4920 and 6430 A. For angles 
near the limb, considerable disagreement was found. 

Sinton and Strong (1959) obtained Martian spectra in the 8 to 13^ region. From these spectra the 
C0 2 bands at 9.4, 10.4 and 12.6/x were identified, and it was also concluded that less than 20% of 
the Martian surface materials contained silicates. The conclusion concerning silicates was based 
on the absence in the spectrum of residual ray features in the 8 to 9/x region. However, Van Tassel 
and Salisbury (1964) have shown that fine-grained silicates and powdered limonite emit like gray- 
bodies with little indication of the residual ray features. Low thermal conductivities, derived 
from the diurnal variation of surface temperature, suggest the Martian surface material is small 
grained (Leovy, 1965). 

The spectra of Figure 23 corresponding to the limonite and granite samples do not match the 
observed spectra of Sinton and Strong (1959). The granite surface does not compare because of 
the emissivity features in the region of 850 cm" 1 . The limonite surface comparison fails because 
the emissivity features in the region of 9.4 and 10.4yu enhance the absorption by the 9.4 and 10.4^ 

C0 2 absorption bands with respect to that of the 12.6/x band. Thus the strengths of the three bands 
would not be in the observed ratio. Smaller grained samples of limonite and granite, different 
surface materials or combinations of surface materials probably would give a better comparison 
with observations. However, the purpose here is not to reproduce the observed spectra of Sinton 
and Strong, but to illustrate the magnitude of the effect of surface emissivity on the spectrum. It 
can be seen that the surface emissivity is the dominant parameter in determining the shape of the 
spectrum and that the brightness temperature may be considerably lower than the radiating tem- 
perature of the surface. 

For the 6.3 f± H 2 0 band, the main effect of the surface water of hydration features of limonite 
and granite is to lower the predicted brightness temperatures without seriously affecting the band 
contour or the 17 cm” 1 periodicity of the rotational structure. The band contour and rotational 
structure is fairly preserved because the surface emissivity changes only slowly with wave 
number, as it exhibits no rotational structure. To identify H 2 0 in the atmosphere of Mars from 
measurements in this spectral region requires that any surface water of hydration be taken into 
account. One possible method for identifying atmospheric H a O under these conditions depends on 
the rotational structure of the atmospheric H 2 0 which is not present in the surface water of hy- 
dration. Qualitatively, the 17 cm" 1 periodicity of the rotational structure should allow a separation 
of the effects of atmospheric H 2 0 and surface water of hydration. This also has been accomplished* 
using cross correlation techniques. 

*B. J. Con rath . Private communication, 1967. 
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VII. CONCLUSIONS 


The main object of this investigation has been to theoretically calculate the outgoing spectral 
radiance from a planetary atmosphere for spectral intervals of 5 cm" 1 . The theoretical radiances 
provide an aid both in the interpretation of high spectral resolution planetary radiation data and in 
instrument design. The high spectral resolution spectra presented for the model atmospheres of 
Earth and Mars show that considerable information is available to recover atmospheric and sur- 
face parameters. 

Other more detailed results are: 

(1) An apparent error in the transmittance tables of Wyatt, Stull and Plass (1962b, 1963) was 
found. Disagreement with the WSP tables is mainly for the C0 2 calculations. It was also 
found that the theoretically calculated integrated absorptances for the 6.3m absorption band 
were systematically smaller than the corresponding experimental values. 

(2) For the Earth, the effect of surface emissivity in the 10.5 to 11.5m "window” is to lower the 
brightness temperature about 2°K. 

(3) Identification of the minor constituent H 2 0 in Mars T atmosphere from the weak absorption 
of the 6.3^ band will be very difficult. The difficulty increases due to possible surface 
water of hydration. 

The direction of future efforts is very clear. The atmospheric absorption by all molecules of 
importance should be calculated using the techniques of Dr ay son (1966), thus eliminating the band 
model approximations and Curtis-Godson approximation. The pressure-broadened line half-width 
should be considered as a function of rotational quantum number and in addition it is important to 
determine the appropriate pressure-broadened line shape. 
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Appendix A 

List of Symbols 


a (ao) - Spectral amplitude of emitted radiation. 

a' - Constant in Benedict modification of Lorentz line shape. 
b(v) - Relative line shape. 
b 1 - Impact parameter, 
bg - Weisskopf impact parameter. 

b" - Constant in Benedict modification of Lorentz line shape, 
c - Velocity of light, 2.997925 x 10 10 cm sec -1 . 
d - Constant in Benedict modification of Lorentz line shape, 
e - Charge of electron, 1.60210 x 10" 19 coul. 
f - Final quantum state of transition, 
f ' - Oscillator strength, 
g - Gravitational acceleration, 
h - Height above planetary surface. 

-ft - Planck constant/277. 
i - Initial quantum state of transition . 

k - Absorption coefficient per unit length per unit pressure, cm" 1 
k m - Mass absorption coefficient, cm 2 /gm. 
k b - Boltzmann’s constant, 1.38054 x 10" 16 erg deg -1 . 
k 0 - Mixed line parameter, Equation 64. 
m - Mass. 

m' - Quantum number. 

m" - Coefficient, Equation 53. 

m f - Degeneracy index for state f . 


atm 1 . 
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n 


Number of lines. 


n. - Number of lines in i th intensity decade. 

n. - Total number of lines in subinterval s . 
n' - Quantum number, 
n" - Temperature coefficient, Equation 52. 
n M - Constant, Equation 13. 
p - Partial pressure. 

q v - Volume fraction, ratio of absorber pressure to total pressure. 
q m - Mass fraction, ratio of absorber density to total density, 
q* - Mixing ratio, ratio of density of water vapor to density of dry air. 
q 0 - Vertical distribution parameter for 0 3 , cm-atm/km. 
r - Geocentric distance. 

r' - Distance between radiating and perturbing particle. 
r Q - Geometrical parameter, 
r* - Spectral reflectivity of planetary surface, 
s - Slant path for integration of optical path length, 
s' - Dummy integration variable, 
t - Time. 

t' - Dummy variable. 
u - Optical path length. 

u\ u" - Optical path length at arbitrary level in atmosphere. 

u s - Optical path length from top of atmosphere to planetary surface, 
v - Molecular velocity. 

v r - Relative velocity of perturbing particle. 
v r - Mean relative velocity. 

v p - Most probable speed of Maxwellian distribution. 
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x' - Mixed line shape parameter, Equation 67. 
x - Complex displacement of electron. 


X 


0 


y - 


y' “ 


A - 
B - 


q'cn - 

c - 

D - 



H(y', x') - 
I e (cS) - 


i;c«) - 

*'t - 

I>) - 

Km - 


a^cta, = 1) - 


i 

n 


J - 


J 2 - 


Acceleration along x. 

Displacement at t = 0, 

Change of variable parameter, Equation 114 
Mixed line shape parameter, Equation 66. 

Change of variable parameter, Equation 115. 

Change of variable parameter. 

Self-broadening coefficient. 

Planck function. 

Normalization constant. 

Constant, Equation 13. 

Optical collision diameter. 

Downward monochromatic radiant emittance. 
Instantaneous power radiated by an accelerating electron. 
Total radiated power. 

Relative mixed line shape. 

Spectral emission intensity per oscillator. 

Spectral emission intensity per unit volume. 

Integrated emission intensity. 

Outgoing spectral radiance at level s. 

Outgoing spectral radiance at level u. 

Outgoing spectral radiance in interval Av at level where 
Bessel function of imaginary argument and order n. 
Rotational quantum number. 

Total quantum state of perturbing particle. 




M' - 



N' - 
N(v. ) - 

P - 
P - 

e 

P - 

s 

P o“ 
P - 
P(s.) - 

R - 
R - 

p 

R - 

s 

s T - 
s - 
s m - 
S(b',J 2 ) - 

s - 

T - 



Molecular weight. 

Molecular weight of atmosphere. 

Molecular weight of absorbing gas. 

Constant, Equation 44. 

Magnetic quantum number of perturbing particle. 

Number of colliding parti cles/unit volume. 

Number of absorbing particles/unit volume. 

Probability distribution function for line position v. . 

Total pressure. 

Effective pressure. 

Surface pressure. 

Reference pressure. 

Mean pressure over variable pressure path. 

Probability distribution function for line of strength s. . 

Universal gas constant, 8.3143 x 10 7 erg deg" 1 mole" 1 . 

Planetary radius. 

Lower integration boundary. 

Level representing top of atmosphere. 

Integrated absorption coefficient, cm" 2 atm" 1 , wave number units. 
Integrated absorption coefficient, cm 2 gm" 1 sec" 1 , frequency units. 
Probability of induced transition. 

Average integrated absorption coefficient. 

Temperature. 

Time between collisions. 

Mean time between collisions. 

Reference temperature. 
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Tj - Isothermal temperature. 

VjCt) - Perturbing potential due to collision. 

W - Equivalent width, 
x - Dummy variable. 

a - Mean half -width, collisional broadening. 
a L - Mean half-width, for Lorentz line shape. 
a D - Half -width of Doppler line shape at half- maximum. 
a L (v) - Half-width of Lorentz line for molecular velocity v . 
a w - Weisskopf half-width. 
a Q - Half-width at reference conditions. 
j3 - Line shift. 
y - Force constant. 

8 - Spectral sub- interval. 

S A - Line shift in Anderson theory, 
e - Change of variable parameter, Equation 119. 

e M - Specifies number of spectral intervals to be considered on either side of interval 
of interest. 

e® - Spectral emissivity of surface. 

L V - Change of variable parameter, Equations 124, 125. 
v - Change of variable parameter, Equation 118. 

? 7 ' - Total phase change. 

Q - Zenith angle at level s. 
k - Wavelength. 

- Natural wavelength. 
fj. - Dipole moment. 
v - Wave number, cm” 1 . 
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Center wave number of starting interval S . 

Center wave number of final interval S . 

Natural wave number. 

Frequency. 

Natural frequency. 

Frequency increment. 

Change of variable parameter, Equation 116. 

Defined by Equation 121. 

Change of variable parameter, Equation 117. 

Total atmospheric density. 

Absorber gas density. 

Boltzmann distribution for perturbing particle in state J 2 
Absorber density at reference condition. 

End time of time interval. 

Spectral response of measuring instrument. 

Change of variable parameter. 

Angular frequency. 

Natural angular frequency. 

Transition frequency between states m # and n \ 

Dummy variable. 

Perturbed angular frequency. 

Angular frequency increment. 

Instantaneous total phase shift. 

Spectral division. 

Collision cross-section. 


Effective collision cross-section. 



Partial collision cross-section. 


A 


J 


2 


T v (u\ u') - 

- 

T(v) - 
T.(v) - 

V-) - 

- 

KM - 

- 

$ _ 
o “ 


Spectral transmittance between level u" and u' * 

Spectral transmittance at lower boundary surface. 

Average transmittance at v. 

Contribution to transmittance at v from n. lines in S.. 

} j 

Total average transmittance at v in interval S k , 

Average transmittance over interval . 

Average transmittance at v in interval S k , shifted mesh. 
Average transmittance at v in interval S k , unshifted mesh. 
Probability distribution for time between collisions. 

Defined by Equation 126. 




Appendix B 

Theoretical Brightness Temperatures and Weighting Functions for the Tiros VII and 
Nimbus II 15 M Carbon Dioxide Absorption Band Channels 

Both the Tiros VII and Nimbus II meteorological satellites carried a radiometric experiment 
which responded to thermal radiation in the 15^ C0 2 absorption band. The spectral response of 
the instrument on each satellite can be found in the "Tiros VII Radiation Data Catalog and Users’ 
Manual" (1964) and the "Nimbus n Users’ Guide" (1966). 

The radiometers are calibrated in terms of brightness temperature T BB , which is related to 
the measured effective radiance N through the expression 

f l/=CO 

^VT BB )d^, (Bl) 

v=0 

where 4> v is the instrumental spectral response. The outgoing radiation measured by the instru- 
ment can be written as 

N(Top) = ^ , (B2) 

i = l 

where N Av is the effective radiance from the i th spectral interval Av and the summation covers 
the spectral range over which the instrument responds. 

The effective radiance for each spectral interval can be obtained from Equation 86 of the main 
text in terms of height along the vertical h as 

N Av (Top) = T Av (h = 0) f 4> v Ij(h = 0) dv + f P f 4> V B V (h) dv 9TAl/ (H - 9s £ g) dh. (B3) 

‘'Av *n~0 ‘ttv ' ’ 1 


The above equation can be rewritten as 


N Av (Top) = T a ^(0) 



<^r(0) dv + 



</>' (h) dh 


(B4) 


where \p‘( h)is given by 


<P‘ (h) = 


^ Av (h) 3s(ht g) \ r 

3s 3h / -L 




dv . 


(B5) 
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The quantity 0 # (h) represents the thermal emission contribution of each atmospheric layer to the 
observed outgoing radiation. A normalized weighting function 0(h) which can be defined by 

0(h) = 0'(h)/N(Top) (B6) 


represents the contribution of each atmospheric layer relative to the total outgoing radiation. 

The effective radiance N(Top) and the weighting function 0(h) have been calculated for a set of 
model atmospheres for the Earth using the direct spectral integration program of Drayson (1966). 
The temperature profiles, the weighting functions for vertical viewing and the weighting functions 
for a ground zenith angle of 70° are shown in Figures B1 and B2 for the Tiros VII and Nimbus II 
radiometers respectively. Brightness temperatures for these cases are presented in Table Bl. 



Figure Bl — Temperature profiles used in the calculations 
are shown at left. The weighting functions are shown 
for vertical viewing and for a ground zenith angle of 
70° for the Tiros VII 1 5/a absorption band channel. 



Figure B2— Temperature profiles used in the calculations 
are shown at left. The weighting functions are shown 
for vertical viewing and for a ground zenith angle of 
70° for the Nimbus II 15/a absorption band channel. 


Table Bl 

Brightness Temperatures for Tiros VII and Nimbus II 15/a Absorption Band Channel. 


Model 

Tiros VH 

Nimbus II 


<9=0° 

II 

o 

0 

0 = 0 ° 

II 

O 

0 

High lat., summer 

235°K 

236°K 

237°K 

236°K 

Tropical 

220°K 

224°K 

228°K 

226°K 

1962 U.S. Std. 

223°K 

224°K 

227°K 

225°K 

High lat., winter 

212°K 

211°K 

217°K 

214°K 
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Appendix C 

Atmos 1 - IBM 7094 Fortran IV Program for Computing Optical Path Length 


SID JOIT V.K. OUTGOING ATMOS RADIATION SECT 1-OPTICAL PATH 

SPAUSE 

SDATE 012766 

SEXECUTE I B JOB 

S I B JOB GO 

SIBFTC ATMOS 1 M94.XR7 .NODECK .NOLI ST .NOREF 

COMMON GO.XMINH.P1* RAD .CAPK *C 1 .C2.DELR.DELHO.HT .XMAXH *HC .TO ♦ I GAS . I 
1 ATM/ZA/ I XA » DATA (200 )/ZB/IXB» DATB ( 100J/ZC/I XC .DATC ( 100 ) 
COMMON/ZD/HT3. LAYER 1 .LAYER2 . THLAY 1 . THLAY2 . THLAY3 
DIMENSION PLKA < 4 ) .PLKB(4) .PLK1<4) .XT A (7*90 *6) .WORDS! 11) 

DIMENSION PRESS! 100 ) .TEMP ( 100 ) »HH ( 100 ) »HM I N ( 100 ) .PRESH! 100) .PREHH! 
1100 ) »XP ( 1200) .XT1( 1200 ) .XDEL ( 1200 ) . X < 1200) .XI! 1200) 

DIMENSION XT2! 1200) .XPBOT! 1200) 

1 FORMAT (3 13) 

2 FORMAT ( 1H0.10X»F8.2>6X.F12.9.6X.F7.2»6X.F12.7,6X.F12.7*6X.F11.7.6X 
1 .F12.4/27X.F7.2 . 8X . F 12 . 6 . 6X . 1 PE12 . 5 . 5X > 1PE 1 2 . 5 . 5X .0PF8 . 5 . 
26X.0PF12.4/27X.1PE13.6.5X.1PE13.6) 

3 FORMAT (13) 

4 FORMAT ( 1 0 ( F5 • 2 ) ) 

5 FORMAT (51F3.1.F7.5) ) 

6 FORMAT ( 1H .10X.4HK = F 1 1 . 7 . 3X . 7HK ( 1 ) = F4.2 »3X .7HM 2 ) = F4.2) 

7 FORMAT (I3*I3»I3*I3»(E13.7) ) 

8 FORMAT (11A6) 

9 FORMAT ( 1H » 10X . 8HANGLE = F6 . 2 .3X . 7HXMI NH= F8» 2 * 3X »26H LOWER INTEGRA 
IT I ON LIMIT = F7.2) 

10 FORMAT ! F 5 • 2 ) 

11 FORMAT ( 1H .415) 

12 FORMAT (F6.2.I3. (F10.6) ) 

13 FORMAT (3E12.6.F6.2.E12.6.F6.2 ) 

14 FORMAT ( 1 H . 10X « 12HCASE NUMBER 13) 

15 FORMAT ( 1H0.4( 1PE13. 6 .5X )/ ) 

16 FORMAT ( 1H0.4HPLK=1PE13. 6) 

17 FORMAT ( 1H0.5HSF = F7.3) 

18 FORMAT (1H »10X»F12.6.3X.F12.6.3X.F12.6»3X.F12.6»3X.F12.6.3X»F12.6. 
13X.F12.6 ) 

19 FORMAT ( 1H1 .8X.16H HT-BTLAY .15H AVET(H20) .15H AVEPIH2 

10) ♦ 15H DEPTH ( H20 ) »15H AVET ( C02 ) »15H AVEP(C02) »15H 

2 DEPTH! C02) ) 

23 FORMAT (F8.6.F4.2.F4.2.I1) 

24 FORMAT (1H1.8HM I D-HT . 9H M ID-TEMP . 12H MID-PRESS. 8H RHOl . 1 1H 

1 PEFF1 .8H AVET1 ,11H AVEP1 ,12H DEPTH1-BLAY.7H PRAT .11H 

2 PEFF2 . 8H AVET2 .11H AVEP2 .12H DEPTH2-BLAY/50X. 13H!WATER 

3 VAPOR) .29X.17H! CARBON DI-OXIDE)//) 

25 FORMAT ! 1H .10X.4HKK= F11.4.3X * 7HKK ( 1 ) = F4. 2 .3X »7HKK I 2 ) = F4.2) 

26 FORMAT ( 1H . F7 .3 . 2X . F6. 2 »4X . F8 . 3 . 3X .F6 . 3 * IX . F8 . 3 .3X .F5 . 1 »4X . F8. 3 
1.1X.F9.6. 4X *E9«3*3X»F6.1.3X*F5.1*3X*F8«3*3X»F11.4) 

28 FORMAT (5 (3X.F5.2.2X.F9.4) ) 

29 FORMAT ( 1H1 . 10X . 62HOUTGO I NG ATMOSPHERIC RADIATION SECT 1 OPTIC 
1AL PATH LENGTH) 

31 FORMAT ( 6E 12 .5 ) 

32 FORMAT I12F6. 2) 

33 FORMAT !5X.I4.F10.2.E15#5.E12»4.E10«5.E10.5) 

34 FORMAT (5X*E10.5*F10»2) 

35 FORMAT (5X.4I5.E12.5.E12.5.F10.2.F12.4) 

36 FORMAT !5X.I5*F10.5.I5.F10.5) 
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40 FORMAT ( 1H0 *6X *4HGRAV *3X * 7HSURF PR*3X*8HSTD TEMP *2X *6HRADI US * 2X *9HG 
IAS CONST * IX *4HDELR* 2X *4HIGAS * 1X*4HJ6AS» IX *4HI ATM* 1X*4HJTAN* IX *2HLA 
2) 

41 FORMAT (1H *E10. 4* 2F10 • 3 * F10 .4 * E 1 1 . 4 * F6 • 3 * 2 I 5 * 2 1 4 ) 

50 FORMAT (5F6*2*2I5,3F5«2) 

51 FORMAT (F5*2*F10*1*F10«2) 

52 FORMAT ( 1H *F6 • 2 * F10* 2 ♦ F10. 3 ) 

5601 FORMAT ( 1H1 ) 

5605 FORMAT ( 1H *13) 

C OUTPUT TAPE ON A5 FOR HEIGHT *AN6LE*TR ANSM I TTANCE MATRIX 

NUM1»0 # 

HT 2 = 0 • 0 
READ (2*3) NUM 
20 READ (2*3) JTAN 

READ ( 2 * 5 0 ) RAD *XM INH *DELR*HT*HT4*LAYER1*LAYER2*THLAY1*THLAY2*THLAY3 

READ (2 *10) ANGLE 

READ (2 *3 ) ( WORDS ( I ) * 1 = 1*11 ) 

READ (2*13) G0*P1*P0* XMOL * RD ♦ TO 
READ (2*23) CCAPK * CC 1 *CC2 * JGAS 
READ (2 *23 )CCPK * Cl *C2 * I GAS 
I VERT = — 1 

I F ( IVERT ) 60*60*70 

C HEIGHT VS TEMPERATURE READ IN* IVERT =-l*************************** 

60 READ ( 2 * 3 ) I X A 

READ (2 *4) ( DATA ( I ) *I=1*IXA) 

WRITE (3*28) (DATA(I) *I=1*IXA) 

GO TO 80 

C PRESSURE VS TEMPERATURE READ IN* IVERT =+l ************************ 

70 READ (2*3)1 XA 

READ (2 *31 ) ( PRESS ( I) *1=1*1 XA) 

READ (2 *32) ( TEMP ( I ) *1=1* IXA) 

C CONVERSION FROM PRESSURE TO HE I GHT ******************************** 

HMINd ) = 0.0 

CONST = 8* 317E+07/ ( XMOL*GO > 

DO 100 1=2*1 XA 

TERM= ( CONST* ( TEMP ( I ) +TEMP ( I ~1 ) ) * < PRESS ( I ) -PRESS ( I ~1 ) ) ) /(PRESS { I )+P 
1RESS ( 1-1 ) ) 

TERM=TERM/1 *E+05 
HMIN ( I )=HMIN ( 1-1 )+TERM 

WRITE (3* 33) I *TEMP( I ) *PRESS( I ) *CONST * TERM *HMIN ( I ) 

100 CONTINUE 

DO 115 1 33 1 * I XA 

HH ( I ) =HM I N ( I XA ) -HMIN ( I ) 

DATA ( I ) =TEMP ( I ) 

WRI TE ( 3 ♦ 34 )HH( I ) * DAT A ( I ) 

115 CONTINUE 

DO 120 1=1, IXA 
IXAA=IXA+1-I 
I XA1 =S 2*I —1 
I X A2 = 2*I 

DATA ( I XA1 ) =HH ( IXAA) 

DATA ( IXA2 ) =TEMP ( I XAA ) 

PRESH(I)=PRESS( IXAA) 

PREHH ( I ) =1013.25*PRESH( I ) 

WRITE( 3*35 ) I * IXAA* IXA1 * IXA2 *PRESH( I ) *DATA ( IXA1 ) *DATA ( IXA2 ) * PREH 

1H( I ) 

120 CONTINUE 

I XXX=2*I XA 
IH=HMIN( IXA) 

I H= I H-l 

HT 4 = FL0AT ( I H ) 

HT = HT 4 

WR I TE ( 3 * 36 ) IH * HT4 * I H*HT 

I X A=2* I X A 
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80 READ ( 2 . 3 ) IXB 

READ (2*5) ( DATB (I)»I=1*IXB) 

WRITE(3»28) (DATB(I) »I=1.IXB> 

READ (2*3 ) IXC 

READ (2*5) (DATC(l) .1=1.1 XC) 

WR I TE ( 3 * 28 ) ( DATC (I) .1 = 1. IXC) 

JTA=1 

CAPK=(XMOL>/< 12.45*18. ) 

CCAPK=CCAPK* ( 1 . E+05 J 
30 ANG= ANGLE* .0174 5 3 294 
HT3 = HT 4 
I TA=0 

NUM1 = NUM1 + 1 
CALL SLITECO) 

IDELR1 = ( DELR*1000. /2 . ) +.05 

IDELR=DELR*1000.+.5 

IH=0 

RO= ( RAD+XMINH) *S I N ( ANG ) 

SUM=0. 

SUMBOT=0 .0 

C PRESSURE INTEGRA! ION********************************************** 

WRITE! 3.5601 ) 

N= 1 

xm=o.o 
XI ( 1 )=0.005 
Hl= • 005 
H=0.0 
DEL=0. 

DELL=0.005 
GO TO 608 

600 IF(N-51)601. 602.605 

601 IH= IH+10 
IH1= IH+05 

X(N)=FLOAT( IH)/1000. 

XI ( N ) = FLOAT ( IH1 )/1000. 

DELL=.01 
DEL=DELL 
GO TO 608 

602 IH= IH+10 

I Hl= IH+50 

X ( N ) =FLOAT ( IH ) / 1000 • 

XI ( N ) =FLOAT ( IH1 ) /1000. 

DELL = .05 5 
DEL= .01 
H1=X1 ( N) 

H = X ( N ) 

GO TO 608 
605 I H= I H+I DELR 

I HI = I H+I DELR1 
X ( N ) =FLOAT ( IH) / 1000 • 

X 1 ( N ) = FLOAT (IH1)/1000» 

DELL=FLOAT ( IDELR)/1000. 

DEL=DELL 
HI =X 1 ( N ) 

H=X ( N ) 

608 CALL INTA(Hl.THl) 

CALL INTA(H.THBOT) 

XT 1 ( N ) =TH1 
XT2 ( N ) =THBOT 

GR= G0*{ RAD/ ( H 1+R AD ) )**2 
GRBOT =G0* ( RAD/ ( H+RAD ) ) **2 

SUM = SUM+DELL* < GR*XMOL ) * ( 1 . E+05 ) / ( RD* TH1 ) 

SUMBOT =SUMBOT + DEL *( GRBOT *XMOL )*( l.E+05)/(RD*THBOT) 

XP(N)=P1*(2 .71828182** ( —SUM ) ) 

XPBOT (N)=P1*(2.71828 182** ( -SUMBOT ) ) 



R=RAD+X(N) 

I F ( R— RO )610, 610,612 
610 R=RO 

H= RO— RAD 
I H = H* 1 000 • 

GO TO 600 

6 12 XDEL ( N ) = ( R/SQRT < R**2— R0**2 ) ) *DELL 
IHT= <HT-DELR)*1000. 

N = N+ 1 

I F ( IH-IHT ) 600 ,604 >604 
604 NOL AY=N- 1 
I PR 1 = 0 

WRITE(3,5605) NOL AY 
I PR2= IPR 1+55 
WR I TE ( 3 , 29 ) 

WRITE(3»14)NUM1 

WRITE (3,8) ( WORDS ( I ) , I =1 , 11 ) 

WR I TE ( 3 • 9 ) ANGLE »XMINH»HT2 
WRITE( 3* 6 )CAPK,C1 ,C2 
WR I TE ( 3 * 2 5 ) CCAPK * CC 1 , CC2 
WR I TE ( 3 ,40 ) 

WR I TE ( 3,41 )G0,P1,T0,RAD»RD,DELR, IGAS,JGAS, I ATM, JT AN 
WR I TE ( 3 * 24 ) 

DELU = 0 • 

T SUM 1 = 0 • 

PSUM 1 =0 • 

SUM2 =0 • 

DELUU=0. 

T SUM2 = 0 • 

PSUM2 = 0 • 

SUM2 2 = 0 • 

SUM4=0 • 

SUM44=0 • 

C OPTICAL DEPTH COMPUT AT I ON ************************************* **■»•* 

DO 650 M = 1 , NOL AY 
N=NOLA Y-M+ 1 
H = X ( N ) 

HI »X 1 ( N ) 

TH1=XT1( N) 

PR =XP ( N ) 

DELS=XDEL ( N ) 

CALL I NT B ( H 1 , RHO I ) 

CALL I NT C ( H 1 , RHO J ) 

SUM4=DELS*RH0I* ( (PR/PO )**C1 )* ( < T0/TH1 ) **C2 ) 

DELU1=SUM4*CAPK 

DELU=DELU+DELU1 

T SUM 1 = TH1 *D£LU1+TSUM 1 

CALL PPEF V ( RHOI , PR, I GAS , CAPK , PEFF 1 ) 

PSUM 1 =PEFF 1 *DELU1+PSUM 1 
SUM2 =SUM2+SUM4 
DEPT 1 = CAPK* SUM 2 
AVET1=TSUM1/DELU 
AVEP 1 =PSUM1 /DELU 

SUM44=DELS*RH0J*( (PR/PO ) **CC 1 ) * ( ( T0/TH1 )**CC2 ) 

DELU2 =SUM44*CCAPK 
DELUU=DELUU+DELU2 
TSUM2=TH1*DELU2+TSUM2 

CALL PPE F V ( RHOJ , PR, JGAS, CCAPK ,PEFF2 ) 

5UM22 = SUM22 +SUM44 
DEPT 2=CCAPK*SUM22 
PSUM2=PEFF2*DELU2+PSUM2 
AV ET2 = TSUM2 /DELUU 
AVEP2=PSUM2/DELUU 
PRAT -PR/PO 
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780 I PR I s I PR 1 + 1 

I F ( IPR2-IPR1 ) 782 »782 ,784 
782 I PR2 = I PR 1 + 55 
WR I TE ( 3 ♦ 24 ) 

784 WRITE! 3»26) HI , TH1 ,PR » RHO I .PEFF1 , AVET1 .AVEP1 .DEPT1 .PRAT .PEFF2 ♦ 

1AVET2.AVEP2 .DEPT2 
I F ( H-HT3 >750.750.650 

C OPTICAL DEPTH MATRIX******************************************'**** 

750 I TA= I TA+ 1 

PUNCH 51.H.XT2IN) .XPBOT(N) 

752 XTA( 1 »ITA»JTA)=H 

XTA(2»ITA»JTA)=AVET1 

XTA(3.ITA*JTA)=AVEP1 

XTA<4»ITA»JTA)=DEPT1 

XTA(5.ITA.JTA)=AVET2 

XTA!6.ITA,JTA)=AVEP2 

XTA(7»ITA.JTA)=DEPT2 

IHT3=»HT3 

I F ( IHT 3- 0)706.706. 756 
756 IF( ITA-LAYER1 >762.758.758 
758 I F ( I T A— LAYER2 >764.766.766 
762 HT3=HT3-THLAY1 
GO TO 777 

764 HT3=HT3-THLAY2 
GO TO 777 

766 HT3=HT3-THLAY3 
777 I HT3 =HT 3 

I F ( IHT 3- 0)770, 770. 65 5 
770 HT 3 = 0 . 

655 I F ( HT 2— H >650.706,706 
650 CONTINUE 

706 IF( JTA-JTANJ720, 708.708 
720 JTA= JTA+ 1 

READ (2.10) ANGLE 
GO TO 30 
708 I DONE= 1 
LA= I TA 

WRITE! 5.8 ) (WORDS! I ), 1 = 1 ,11 ) 

WRITE (5, 12) PI. IX A, (DATA! I) ,1 = 1. IXA) 

WRI TE! 5.7 ) HT2 ,NUM , J TAN . L A , ( ( ( XTA ( K T ,MT , NT > ,KT=1 ,7 ) ,MT = 1 .LA) »NT=1 , 

1 JT AN ) 

WRITE! 3.19) 

WR ITEI3.il) HT2.NUM, JTAN.LA 

WRITE<3»18> ( ( (XTA(KT»MT .NT) »KT=1»7) ,MT=1.LA) ,NT=1»JTAN) 

NUM=NUM— 1 

IF(NUM) 710.710,20 

710 WR I TE ( 3 ♦ 52 ) ( X ( N ) ,X T2 ( N ) , XPBOT ( N ) , N = 1 ,NOL AY ) 

STOP 

END 

67 CARDS 

SIBFTC PPEFV M94.XR7.NODECK.. NOLIST, NOREF 

SUBROUTINE PPEFV ( RHO .PR » I GAS . CAPK , PEFF ) 

PBl^l.3 

PB2=6.3 

I F ( I GAS- 2 ) 10,20,30 

10 PEFF=PR+ ( PB1-1 « )*CAPK*10.**(-5)*PR 
GO TO 40 

20 PEFF=PR+ (PB2-1. ) *10. ** ( -3 ) *RHO*12 . 45*CAPK*PR 

GO TO 40 
30 PEFF=PR 
40 RETURN 
END 

12 CARDS 
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Ill 


I ■ III II I I I I I 


$ I BMAP 
INTA 


MOD 

SAVA 

ZA 

I X A 
DATA 

TEM 

INTB 


MODI 

ZB 

I X B 
DATB 

TEM1 

INTC 


MOD 2 

ZC 

IXC 

DATC 

TEM2 


INTA 

SAVE 

100 * NODECK 

SXA 

SAVA * 4 

CLA 

I X A 

ALS 

18 

STD 

MOD+2 

CLA 

4*4 

STA 

TEM 

CLA* 

3*4 

TSX 

TIN1.A 

PZE 

DATA 

PZE 

1 * *** 

AXT 

** * A 

STO* 

TEM 

RETURN 

INTA 

CONTRL 

ZA 

USE 

ZA 

BSS 

1 

BSS 

100 

USE 

PREVIOUS 

BSS 

SAVE 

1 

SXA 

SAVA* A 

CLA 

IXB 

ALS 

17 

STD 

MOD1+2 

CLA 

A* A 

STA 

TEM1 

CLA* 

3*4 

TSX 

T I N 1 * 4 

PZE 

DATB 

PZE 

1* ♦** 

LX A 

SAVA *4 

STO* 

TEM1 

RETURN 

INTB 

CONTRL 

ZB 

USE 

ZB 

BSS 

1 

BSS 

100 

USE 

PREVIOUS 

BSS 

SAVE 

1 

SXA 

SAVA *4 

CLA 

IXC 

ALS 

17 

STD 

MOD2+2 


CLA A* 

4 

STA 

TEM2 

CLA* 

3*4 

TSX 

T I N1 * 4 

PZE 

DATC 

PZE 

1* ♦** 

LX A 

SAVA * A 

STO* 

TEM2 

RETURN 

INTC 

CONTRL 

ZC 

USE 

ZC 

BSS 

1 

BSS 

100 

USE 

PREVIOUS 

BSS 

1 

END 



62 CARDS 
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SIBMAP 

* 

* 

TIN1 


TIN1 

125 *M94 *NODECK 

INTERPOLATION SUBROUTINE 
FAP 

ENTRY 

TIN1 

STO 

TIN1+98 

SXD 

T I N 1 + 87 * 1 

SXD 

T I Nl+88 * 2 

SXD 

T I Nl+89 ,4 

CLA 

2*4 

STO 

TIN1+94 

ADD 

T I N 1+91 

PAX 

0*1 

SXD 

TINl+29,1 

ALS 

1 

STA 

TIN1+95 

CLA 

TIN1+94 

ARS 

17 

SUB 

TIN1+95 

PAX 

0*1 

SXD 

TIN1+31.1 

ADD 

TIN1+95 

PAX 

0*1 

ADD 

1*4 

STA 

T I Nl+23 

STA 

TIN1+41 

STA 

T I N 1 + 43 

CLA 

TIN 1+98 

CAS 

0*1 

T IX 

♦*-1*1*2 

TRA 

TIN1+26 

CLA 

LBT 

T I N 1 + 9 4 

TRA 

T I N 1+4 1 

T I X 

T I N 1 + 3 1*1*0 

LXD 

T I Nl + 90 * 1 

TXL 

T I Nl + 33 * 1 * 0 

LXD 

TIN1+31 *1 

PXD 

0*1 

ARS 

CHS 

18 

ADD 

TIN 1 + 2 3 

STA 

TIN 1 + 54 

ADD 

TIN 1+9 1 

STA 

TIN1+52 

TRA 

TIN1+50 

CLA 

0*1 

T X I 

T I N 1 + 43 * 1 * 2 

FAD 

0*1 

LRS 

35 

FMP 

TIN1+92 

CAS 

T I N 1+98 

TX I 

T I N 1 + 29 * 1 * - 1 

T X I 

TIN1+29*1*-1 

TX I 

TI Nl + 29 * 1 *-~3 

LXD 

T I Nl+90 * 2 

LX A 

T I N 1 + 95 * 1 

CLA 

0*1 

STO 

COM+1 *2 

CLA 

0*1 

FSB 

TIN1+98 

STO 

COM, 2 

TX I 

TINl+58,2 , —2 

T I X 

T I N 1 + 52 * 1 * 2 

CLA 

COM 


TIN1 
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STO 

COM * 2 

LX A 

T I Nl+94 *4 

LXD 

T I Nl + 93 * 2 

TXI 

T I N 1+64 * 2 * — 2 

PXD 

0*2 

PDX 

0*1 

CLA 

COM+2 * 1 

FSB 

COM ♦ 2 

TZE 

TIN1+79 

STO 

T I N 1 + 97 

LDQ 

COM * 2 

FMP 

COM+3 ♦ 1 

STO 

T I N 1+96 

LDQ 

COM+2 *1 

FMP 

COM+ 1*2 

FSB 

T I N 1+96 

FDP 

T I N 1 + 9 7 

ST Q 

COM+3 *1 

TXI 

TIN 1 + 66 * 1 *-2 

CLA 

COM+2 *2 

STO 

COM+2 *1 

T I X 

T I N 1 + 63 * 4 * 1 

CLA 

COM+1 *1 

LXD 

TIN1 + 87 * 1 

LXD 

T I N 1 + 8 8 * 2 

LXD 

T I Nl + 89 *4 

TRA 

3*4 

HTR 


HTR 


HTR 


HTR 


HTR 

f— < 

* 

o 

« 

o 

DEC 

• 5 

HTR 

0*0*2 

HTR 


HTR 


HTR 


HTR 


HTR 


PZE 


PZE 


BSS 

40 

END 



I 

001 

001 

6371. 0000. 00000.10070.00064.. 00003200036002. 0000 1.000. 5 

. 0 

EARTH.ARCT WINT-HIGH LAT WINT COLD. GRD ZA=00.0. QH20=C STRAT 
9. 806650 E+021013.250E+0001013.25E+0028. 9700008. 317E+07273. 18 
.000314 1.0 1.01 

.1292 01.001.02 

076 

00000246.0001.0250.0002.0254.0003.0248.9004.0243.9 
005.0238.8007.0228.6008.0223.6009.0218.5010.0217.2 

014.0211.8015.0210.4015.5210.0024.5210.0025.0209.4 

026.0207.5028.0203.8029.0201.9030.0200.0031.0202.9 

035.0214.5040.0229.0045.0243.5049.0255.1049.5256.6 

050.0258.0050.5258.0055.0258.0055.5257.4056.0256.7 

058.0254.1060.0251.5062.0248.9064.0246.3066.0246.4 
070.0251.8075.0258.6080.0248.6 

032 

000.348 1.0.30 2.0.26 3.0.22 3.7.19 

4.0.17 5.0.1 6.0.059 6.3.05 7.0.026 
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8.0.01 9.0.0051 10. .0024 112.0011 115.0011 

700.0011 

004 

0001.0 80.1.0 

' EOF E SY 

133 CARDS 
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Appendix D 

Atmos 2 - IBM 7094 Fortran IV Program for Computing Quasi-random 
Transmittance for Water Vapor 


SID JOIT V K.UNDE OUTGOING ATMOS RAD SECT 2-TRANS *MI XED LINE.H20 

SPAUSE 

SDATE 010666 

SEXECUTE IB JOB 

SI B JOB GO. MAP 

SIBFTC ATMOS2 M94.XR7 .NODECK » L I ST 
COMMON/AZ/A ( 5 ) .NA(5) 

COMMON/ AA/WORDS (11) .WWAVEI .WAVEF » I NT ER »HT 2 .HT3 .ANGLE .NUM »JTAN . LA 
COMMON /Z A/ I XA.pATAI 200 ) 

COMMON/TT/TRAN < 300 ) 

DIMENSION X I NOT (5).XILINE(5) .XIST (5) . TIS<5). 

1CL A ( 7 ) . CLB ( 7 ) .CPA (4) .CPX ( 4 ) . A I JS ( 5 ) . T I JS ( 5 ) .CLA1 ( 7 ) ,CLB1 (7) 

2>H20( 13000) »TRAND(300) .TRANW ( 300 ) .TRANS (300) .TRANlt 300) 

3 »GEOM (7.56.3) . AA ( 15 ) » AG ( 6 ) 

DIMENSION CLAI 1(7) *CLBI 1(7) 

.999 FORMAT (F6.2.I3. (F10.6) ) 

1000 FORMAT ( 1H0 . 30X > 44HBASI C PARAMETERS FOR OBTAINING TRANSMITTANCE) 

1001 FORMAT (F3. 2 .F5 . 1 » F5 . 2 . F5 . 3 ,F2 . 1 ,F5 . 3 . F5 . 2 . F 12 .6 , F 5 . 1 , F5 . 1 ♦ 1 2 ) 

1002 FORMAT ( 1HQ . 5X » 16HPRESSURE ( ATM ) = 1PE12.6.3X. 17HTEMPERATURE ( K ) = OP 

1F10.2. 3X . 16HDEPTH! PR CM ) = 1PE1 0 . 4 . 3X » 09HHET = 0PF8.2/6X . 15HI 

2NITIAL WAVE = 0PF13 . 1 »3X . 17HF I NAL WAVE = 0PF10.1.3X. 16HWAVE I 

3NCR. = 0PF10.1.3X.09HDWAVE = 0PF8. 1/6X . 16HSTAN. PR ( ATM ) * 0PF1 

42.4.3X.17HSTAN. TEMP. = 0PF10.2.3X. 16HHALF WIDTH = 0PF10.4/6 
5X.16HEMAX = 0PF12.1 .3X.17HCAI J = 1PE10.4) 

1004 F ORMAT (1H0.46X.31HQUASI =RANDOM BAND TRANSMI TTANCE/52X. 15HWATER VAP 
10R GAS/5 IX . 18H LOREN TZ LINE SHAPE) 

1005 FORMAT (1H0.46X.31HQUASI =RANDOM BAND TRANSMI TTANCE /46X . 24HMI XED LOR 
1ENTZ LINE SHAPE/ 52X »15H WATER VAPOR GAS) 

1007 FORMAT ( 1 H .5X.16HCXI = 1 PE 1 2 . 6 . 3X » 17HZETA = IP 

1E10.4.3X.16HZETA2 = 1PE10.4) 

1008 FORMAT (1H1) 

1009 FORMAT! 1H0/ (5X.7HWAVE = F6 . 1 . 5 ( 3X ♦ 1 PE9 . 2 »0PF5 . 1 ) ) / ( 26X . 7HA I J = IP 
15E13.6) ) 

1010 FORMAT ( 1 1 A6 ) 

1017 FORMAT ( 1H0. 10X.30HERROR FROM STATEMENT NUMBER = 15) 

1020 FORMAT ( 1H0.10X. 39HCOMPUTAT I ON OF DIRECT CONTRIBUTION PART) 

1021 FORMAT (20X.6HOM1 = F9.6.3X.5HOM = F9 .6 . 3X . 7HOM11 = F9.6/> 

1028 FORMAT(1HO/(10X»37HCOMPUTATION OF WING CONTRIBUTION PART)/) 

1029 FORMAT ( 1H0 . 10X ♦ 14HCASE NUMBER = 12) 

1030 FORMAT ( 5 ( F10.6 ) ) 

1031 F ORMAT ( 1H .5X.16HBE1 = F12.5. 3X.17HBE2 = FlO 

1 » 5 » 3X > 16HBE3 = 0PF1 0. 5/6X . 16HPB 1 * 0PF12.3.3X. 

217HPB2 = 0PF10.3) 

1032 FORMAT ( 14 ( F4. 3 ) ) 

1033 FORMAT ( 8 ( F8 .8 ) ) 

1040 FORMAT ( 14IF5.3 ) ) 

1041 FORMATUH ./06X.9HDIRECT = 7 ( F5. 3 . IX . F5 . 3 . 3X ) /20X .4 ( F9. 8 . IX . F9 . 8 .3 
IX) ) 

1042 FORMAT ( 1H .03X.9H WING = 7 ( F5. 3 . 1 X . F5 . 3 . 3X ) /20X ♦ 4 ( F9 .8 . IX .F9. 8 . 3X 
1) ) 

1043 FORMAT (1H0.07X.17H WING — MIDPOINT/) 

1048 FORMATUH . 09X . 7H I TEMP= I3.3X.8HI TEMP1= I 3 » 3X * 8HI WAVEU= I10.3X.7HL 
1UPU1= 1 1 0 . 3X . 7HLOWU1 = 1 10 * 3X ♦ 8HIWAVEL= 110) 
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1049 FORMAT ( 1H1 ,40X , 19HDI RECT CONTRIBUTION) 

1050 FORMAT ( 1H , 1X////40X .17HWING CONTRIBUTION) 

1051 FORMAT l 1H *4X , 4HWAVE ,3X * 6HLAMBDA.2 X . 5HTRANS *6X .4HWAVE .3X » 6HLAMBDA . 
12 X » 5HTRANS * 6X , 4HWAVE »3X .6HLAMBDA.2X* 5HTRANS *6X .4HWAVE *3X .6HLAMBDA » 
22X , 5HTRANS.6X .4HWAVE ,3X .6HLAMBDA.2X , 5HTRANS ) 

1052 FORMAT ( 1H , 5 ( 2X >F7. 1 ,F8 . 3 , F9 . 5 ) > 

1053 FORMAT ( 1H . 1X////40X .19HT0TAL TRANSMITTANCE) 

1054 FORMAT ( 1H ,3 X/// 40X.19HFINAL TRANSMITTANCE) 

1058 FORMAT ( 1H0, 10HSTOP AT = 16) 

1061 FORMAT (I3»I3»I3»I3»(E13.7) ) 

1062 FORMAT (£13*7) 

10 63 FORMAT ( 1H1 ♦ 13 * I 3 , 13 , 13 • 7 < E13. 7* 3X ) ) 

1065 FORMAT (5F5.2) 

1066 FORMAT ( 1H1 * 20X * 71HOUTGO I NG ATMOSPHERIC RADIATION SECT2 H20-QUA 

IS I -RANDOM TRANSMITTANCE) 

A5 — GEOMETRY MATRIX TAPE 
A6 — BLANK. TAPE 

A8 — COEFFICIENT TABLE FOR WATER VAPOR-250 WORDS' RECORD-PACKED 
B7 — OUTPUT TAPE FOR THIS PROGRAM-THI S TAPE WILL BE USEDIN ATMOS3 
IDOP=— 1 * LOREN TZ LINE SHAPE**************************************** 
CALL KNUMIX(X*Y *OUT *1) 

READ (2*1030) BE1 . BE2 ,BE3 . PB1 . PB2 

RE AD (2. 1032 ) (CLA(J) * J = 1 » 7 ) » ( CLB ( J ) *J=1.7) 

READ (2. 1033) (CPA(J) * J = 1 *4 ) * ( CPX ( J ) ,J = 1,4) 

READ (2 *1040 ) (CLA1 ( J) *J = 1.7) , ( CLB1U) , J=1 *7) 

RE AD (2. 1001 )HW I DTH* W AVEI * TNOT , PR * DWAVE *PRNOT . TEMPO .PATH tWAVEF * 

1DNU. IWRITE 

C READ TAPE A5 — GEOMETRY MATRIX************************************* 

1080 READ (5*1010) ( WORDS ( I ) . I = 1 , 1 1 ) 

READ (5*0999)P1.IXA. ( DATA ( I ). 1=1*1 XA) 

READ! 5.1061 ) HT 2 *NUM * JTAN .LA * ( < ( GEOM ( KT *MT .NT ) *KT = 1 . 7 ) *MT = 1 *L A ) . 
1NT = 1 . JTAN ) 

WRITE! 3. 1063) I ATM .NUM . JTAN ,LA . ( ( ( GEOM( KT .MT .NT ) .KT = 1.7) »MT = 1 .LA ) » 
1NT=1.JTAN) 

READ (2.1 065 ) ( AG ( JJ ) . JJ=1 . JTAN ) 

LAST = -l 

HT3 = GEOM( 1.1,1 ) 

I DOP=-l 
I BLS=-1 
ICASE=1 
I SE=1 
IER = 1 

C REMEMBER DIMENSION FOR H20(I) IS 13000. 

DO 1096 J = l, 13000 
H20 ( J ) =0 • 

1096 CONTINUE 
1090 MT =LA 

ANGLE=AG( ICASE) 

I XB=MT *2 
I XB1 =MT 
I SEC=-1 

1095 WAVEO=WAVE I 

PATH=GEOM ( 4 ,MT , ICASE ) 

HET=GEOM< l.MT, ICASE) 

TNOT =GEOM ( 2 .MT .ICASE ) 

TNOT=249 . 

PR=GEOM( 3.MT, ICASE) / 101 3. 25 
PATH1 =PATH 

1097 I K =~1 


I WORD=0 
WAVET = WAVEO 
I SHI =WAVEO* 10. 
ISH2=WAVEO 
I SH2= I SH2*10 
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C TEST FOR SHIFTED OR UNSHIFTED WAVE NUMBER ************************ 

I F { ISH1-ISH2) 1870*1098*1099 

1098 I FT 1 =— 1 
I FT2 = 1 

I SB 3 — ! 

SLUMB=0# 

GO TO 1100 

1099 WAVET = WAVEO—2# 5 
SLUMB=1 86 1 • 

I SB= 1 
IFT1-1 
IFT2*— 1 

TEMP =300 - TNOT LARGER THAN 2 50 *********************************** 
TEMP=250* *FOR TNOT WITHIN THE RANGE OF 225#AND 21b. 

TEMP =200 • * TNOT LESS THAN 225. 

1100 I F ( T NOT— 275 • ) 1105*1110*1110 
1105 I F ( T NOT- 225 • ) 1120*1115*1115 
1110 TEMP=300 • 

I T EMP = 1 
GO TO 1300 
1115 TEMP=250 • 

SSW! 300 
1120 TEMP=200. 

I TEMP=-1 
1300 I TT= 1 

C RANGE OF WING CONTR I BUT I ON**************#**********'*********'****** 

MAX=( 4.8000*! 10. )**5*PR*PATH> **< .5)/2. 

EMAX=(MAX+1 )*2 
CX I =PATH /PR 

ZETA= ( 2 • 0/DWAVE ) *HW I DTH* ( PR/PRNOT ) * ( TEMPO/ TEMP ) ** . 5 
ZETA2=ZETA*ZETA 
CAI J=ZETA2*CXI 
WR I TE ( 3 * 1066 ) 

WR I TE ( 3 * 1029 ) ICASE 

WRITE! 3.1010) (WORDS! I ) .1=1*11) 

IF! I DOP) 1220.1220.1225 
1220 WR I T E ( 3 . 1 004 ) 

GO TO 1230 
1225 WR I TE ( 3 . 1 005 ) 

1230 CONTINUE 

WRITE(3.1000> 

WRITE! 3.1002) PR . TEMP .PATH .HET .WAVE I . WAVEF ♦ DNU .DWAVE .PRNOT .TEMPO . 
1HWIDTH.EMAX.CAI J 
WRITE t 3. 1007 ) CXI .ZETA.ZETA2 

WRITE! 3*1041) (CLA! J) »CLB ( J ) . J=1.7) . ( CPA ( J ) .CPX! J) .J=1.4) 

WRITE (3.1042) ! CLA1 ( J ) .CLB1 ( J ) . J=1.7 ) . ( CPA ! J ) .CPX ( J ) . J= 1 . 4 ) 

C CHECK. CORE IF INPUT DATA THERE************************************ 

1294 I WAVEU= ( WAVEF+2.5*EMAX)*10. 

IWAVEL=(WAVEI-2. 5*EMAX ) *10 • 

IF ( IWAVEL-7000 ) 1303.1304.1304 

1303 IWAVEL = 7 000 

1304 IF! IWAVEU-100000) 1306.1306.1305 

1305 I WAVEU= 100000 

1306 I WORDS= ( IWAVEU-IWAVEL)*2/10+10 
INTER= ( WAVEF-WAVE I >/5.0+l. 

WRITE! 3.1048) I TEMP . I TEMPI . IWAVEU.LUPU1 . LOWU1 ♦ IWAVEL 
I F ( I SEC ) 2000.1870.1307 

1307 I F ( I T EMP— I T EMP 1)2000 .1308.2000 

1308 IF! IWAVEU-LUPU1) 1309.1309.2000 

1309 IF ( LOWU1— IWAVEL ) 1310.1310.2000 

C UNSHIFTED OR SHIFTED MESH **************************************** 

1310 IF! ITT-2J1325. 1340.1312 
1312 IER=1310 

GO TO 1870 



n n 


1325 WAVE = WAVEO 
I FT = I FT1 
GO TO 1400 

1340 WAVE=WAVEO+DWAVE/2. 

WAVE I =WAVEI+DWAVE/2» 

WAVEF=WAVEF+DWAVE/2. 

I FT = I FT2 

C PREPARE TO READ INPUT DATA FROM CORE****************************** 

1400 OM 1* 1 • 

I JJ= 1 
I R*1 
1ST 1 =— 1 

IWAVE=WAVE*10.0 
I ST =-l 
GO TO 1900 

C DETERMINATION OF DIRECT CONTRIBUTION TO TRANSMITTANCE************* 

1405 WAVE =WAVE+D WAVE 
I WAVE = WAVE* 10 • 

OM1=1#0 

I F ( WAVE— 1590*0 ) 6301*6301*630 3 
6301 HW IDTH = 0 * 16 
GO TO 6304 

6303 HW IDTH-Q #115 

6304 ZETA= ( 2 • 0/DWAVE ) *HW I DTH* ( PR/PRNOT ) * ( TEMPO/TEMP ) **• 5 
ZETA2=ZETA*ZETA 

CA I J=ZETA2*CX I 
GO TO 1900 
1410 DO 1420 J = 1 *5 

XIST ( J ) =X I NOT < J ) *CA I J 

1420 CONTINUE 

IF( I WRITE) 1415 *1870* 1421 

1421 WRIT£(3*1 009 ) WAVE * < X I NOT ( JN ) *X I L INE ( JN ) * JN” 1 *5)*(XIST( JN ) *JN=1 *5) 
WRITE(3*1020) 

D I RECT CALCULAT I ON************************************************ 
LEG GAUSS QUAD ♦ 7 SUBINT, 8 POINT 

1415 J=1 

1422 CX IST=X I ST { J ) 

IF { IDOP) 1423*1423*1416 

1416 IF ( PR-0.029607 ) 3000 *1423*1423 

C DIRECT CALCULAT I ON * LORENT Z LINE SHAPE**************** 

1423 SUM=0. 

I OM = l 

1424 JOM= 1 

1426 IF(IOM“7> 1430*1430*1428 
1428 GO TO 1480 
1430 CLAI =CLA ( I OM ) 

CL8I=CLBUOM) 

DL.BA=CLBI -CLAI 

PBA=CLBI +CLAI 
SUM1=0# 

1432 IF(J0M“4) 1436*1436*1434 
1434 SUM=SUM+DLBA*SUMl/2 • 

I OM= I OM+ 1 
GO TO 1424 
1436 CPA J= CPA ( JOM) 

CPXJ=CPX< JOM) 

GMIN1 * ( “DLBA*CPX J+PBA ) /2* 

GPUS 1 * ( DLBA*CPX J+PBA ) /2 • 

GPUS=GPUS 1*GPUS1 
GM I NUS=GM I N 1*GMI N1 

SUM1=SUM1+CPAJ* ( EXP ( “CX I ST/ (GPUS+ZETA2 ) )+EXP(— CXI ST/ ( GM I NUS+ZETA2 

1 ) ) ) 

JOM= JOM+1 
GO TO 1432 
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C DIRECT CALC • *M I XED LINE SHAPE. 7 SUBINT. 8 POI NT****************** 

3000 CLBII ( 1) =2.5 
CLAI 1(1) =1.25 
CLBI I (2) =1.25 
CLAI I (2) =0.2500 
CLBI I (31=0.2500 
CLAI 1(3) =0. 1250 
CLBI I (A) =0. 1250 
CLAI I (41=0.0250 
CLBI I (51=0.0250 
CLAII(5)=0.0125 
CLBII (61=0.0125 
CLAI 1(61=0.0025 
CLBII (71=0.0025 
CLAI I (7) =0.0 
SUM=0. 

I0M = 1 
3003 JOM= 1 

BOLTZ=1.38054E-16 
PROTON= 1 .67252E-24 
XM0L=18. 

VEL = 2 .997925E+10 

CONST = SQRT ( 2 .*AL0G ( 2 • 1 *BOL TZ /PROTON 1 /VEL 
DOPHW=CONST*SQRT ( TEMP /XMO 1 *WA VE 
XLORHW=HWIDTH*( PR /PR NOT )*SQRT (TEMPO/TEMP) 

V = ( XLORHW/ DOPHW 1 *SQR T ( ALOG ( 2 . 1 1 
PI=3. 1415927 

QUAN= (PI*Y*CX 1ST 1 /( SORT (P I )*ZETA2 1 
3005 I F ( IOM-7 13015.3015.3010 
3010 GO TO 1480 
3015 CLAI =CLA I I ( IOM) 

CLBI =CLB I I ( IOM 1 
DLBA=CLB I -CLAI 
PBA=CLB I +CLAI 
SUM1 =0 . 

3025 I F ( JOM-4 13035.3035.3030 
3030 SUM=SUM+DLBA*SUMl/5. 

I OM = I OM+ 1 
GO TO 3003 
3035 CPA J=CPA ( JOM 1 
CPXJ=CPX (JOM) 

XP= ( DLBA*CPX J+PBA 1 / 2 . 

XP=XP*SQRT< ALOG (2. 1 1 /DOPHW 
CALL KNUMI X (XP.Y .OUT >2 1 
OUTP=OUT 

XM= (-DLBA*CPX J+PBA 1 /2. 

XM=XM*SQRT ( ALOG (2.1) /DOPHW 
CALL KNUMIXIXM.Y .OUT .2 1 
OUTM=OUT 

XTRANP=EXP( — QUAN*OUTP 1 
XTRANM=EXP ( -QUAN*OUTM) 

SUM1=SUM1+CPAJ* (XTRANP+XTRANM) 

JOM= JOM+ 1 
GO TO 3025 

1480 OM=l .-SUM 

1481 N1=XILINE(J) 

OM 11=1.0 

1485 I F ( N 1-20 1 1486.1486.1490 

1490 OM 1 1 = ( 1 . — OM 1**20 *OM 1 1 

Nl=Nl-20 
GO TO 1485 

1486 OM 1 1 = ( 1 • — OM ) **N 1 *OM 1 1 
0M1=0M1*0M11 

I F ( I WRITE 11488. 1870. 1487 



1487 WR I TE ( 3 * 1 02 1 ) OMl,OM*OMll 

1488 IF(OMl) 1870* 1494 * 1 482 

1482 I F ( J— 5 ) 1483*1494,1870 

1483 J= J+l 

GO TO 1422 

1494 TRAND ( I J J ) = 0M1 

I F ( I JJ-INTER1 1495*1496*1496 

1495 IJJ=IJJ+1 
GO TO 1405 

1496 IF(IWRITE) 1 504 , 1 870 , 1492 
1492 WR I TE ( 3 * 1 049 ) 

WR I TE ( 3 * 105 1 ) 

I RET = 1 
GO TO 1830 

C DETERMINATION OF WING CONTRIBUTION TO TRANSMITTANCE*************** 

1504 IF ( ITT-2 ) 1506*1508*1868 
1506 WAVE I =WA VE I 
I FT = I FT 1 
GO TO 1509 

1508 I FT= I FT2 

1509 IF(IFT) 1501*1870*1502 

1501 TLOW = 700 • 

TUP = 1 0000 • 

GO TO 1510 

1502 TLOW=702.5 

tup = 1 0002 • 5 

C INTERVAL CONTAINING LINES ****************************** 

1510 TWAVE=WAVEI-2.5*EMAX 
DO 1511 JI=1,IJJ 
TRANS ( J I >=1.0 

1511 CONTINUE 
1ST = — 1 

I ST 1 =~1 
GO TO 1518 

1514 TWAVE=TW AVE+DWAVE 

1518 IF(TWAVE-TLOW) 1520 * 1520 , 152 1 

1520 TWAVE=TLOW 

1521 IF(TWAVE-TUP) 1528,1528,1522 

1522 I JJ= I NTER 
GO TO 1759 

1528 LWAVE=TWAVE*10. 

I W AVE=TW A VE* 1 0 . 

IR=-1 

GO TO 1900 

C INCREMENT TO higher wavenumber ************************* 

1530 DNU = 5 0 • 

I DNU=DNU 
ISID€*+1 

IWAVEF=WAVEF*10. 

IWAVEI=WAVEI*10. 

LW AVE= ( L WAVE+ I DNU ) 

IEPSI=IABS< 4*(LWAVE-IWAVE) /100) 

EPS I =FLO AT ( I EPS I ) 

I F ( LWAVE~I W AVEF ) 1535 ,1535,1539 
153 5 I F ( I WAVE I -L WAVE) 1570 ,1570,1530 

1539 LW AVE=TW AVE*10 . 

C INCREMENT TO LOWER WAVENUMBER ************************** 

1540 DNU=~50 • 

I DNU=DNU 
ISIDE=-1 

LW AVE= ( LWAVE+ I DNU ) 

IF (7000-LWAVE) 1542,1542,1514 
1542 I EPS I = I ABS ( 4* ( LW AVE— I WAVE ) /1 00 ) 

EPS I =FLOAT ( I EPS I > 

I F ( I WAVE I — LWAVE ) 1545,1545,1514 
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1545 IF ( LWAVE—IWAVEF ) 1570*1570,1540 

C WING CALCULATION************************************************** 

1570 J=1 

HWIDTH=.10 

ZET A= ( 2 • 0/DWAVE ) *HW I DTH* ( PR/PRNOT ) * ( TEMPO/ TEMP ) **• 5 
ZETA2=ZETA*ZETA 
CAI J=ZETA2*CXI 
T I J 1=1 .0 

C CHECK ON MIDPOINT OR AVERAGE CALC. FOR WING************* 

1580 IF(IBLS) 1584*1870*1582 

C MIDPOINT CALCULATION************************************ 

1582 XIST=XINOT( J) *CA I J 
EPAB=ABS ( EPSI ) 

CORR =EXP ( —BE 2* ( DWAVE/2 • > **BE3*EPAB**BE3+BE2*( 5.0)**BE3) 

TIJ = EXP( — (X 1ST *CORR ) / ( EPAB*EPAB ) ) 

T I JS ( J ) =T I J 
AIJS(J)=X I ST*CORR 
GO TO 1701 

C Leg GAUSS QUAD* 1 INT* 8 POINT************************** 

1584 AIJ=XINOT (J) *CA I J 
AI JS { J ) =AI J 

GO TO 1980 

1585 EPS I = EPS I 
EPSI 1=EPSI 
AI JS( J)=AI J 
CX I ST=A I J 
SUM=0. 

I OM = 1 

1586 JOM = 1 

1588 I F ( IOM-1 >1592*1592*1590 
1590 GO TO 1630 
1592 CL A I -CLA 1 ( I OM ) 

CLB I =CLB 1 ( IOM) 

DLBA=CLBI-CLAI 
PBA=CLBI+CLAI 
SUM 1 = 0. 

1594 I F ( JOM— 4 ) 1598*1598*1596 
1596 SUM=SUM+DLBA*SUMl/4. 

I OM= I OM+ 1 
GO TO 1586 
1598 CPAJ=CPA ( JOM) 

CPXJ=CPX( JOM) 

GMIN1 = EP5I 1- ( -DLBA*CPX J+PBA ) / 2 # 

GPUS1=EPSI 1- ( DLBA*CPX J+PBA ) /2 . 

GPUS=GPUS 1*GPUS 1 
GMINUS=GMIN1*GMIN1 

SUM1=SUM1+CPAJ*( EXP (-CXIST/( GPUS ) ) +EXP { -CX I ST/ ( GM I NUS ) ) ) 

JOM= JOM+1 
GO TO 1594 
1630 TIJ=SUM 

T I JS < J ) =T I J 

1701 N1=XILINE(J) 

TIJ11=1.0 

1702 IF { N 1-2 0 ) 1704 ,1704*1706 
1706 TIJ11=TI J**20*T I J 1 1 

N 1*N1— 20 
GO TO 1702 

1704 TIJ11=TI J**N 1 * T I J 1 1 
T I S ( J ) =T I J 1 1 

1705 TIJ1=TIJ1*TIJ11 
J = J+1 

IF (J-6) 1580 *1720* 1720 
1720 I WAVE I =WAVE 1*10 • 

I JJ* (LWAVE-IWAVEI+50 )/50 
I F ( T I J1-.999) 1722*1754*1754 



1722 TRANW( IJJ)=TIJ1*TRANS( IJJ) 

TRANS ( I J J ) =TRANW ( IJJ ) 

I F ( ISIDE ) 15 AO * 1 870 ♦ 1530 
1754 I F ( ISIDE) 1758*1870*1539 

1758 IF <TWAVE-WAVEF-2.5*EMAX) 1514* 1759* 1759 

1759 I F { I WR ITE )1763 *1870*1760 

1760 WRITE(3*1050) 

WR I TE ( 3 * 1 05 1 ) 

IRET=2 

GO TO 1832 

1763 IF( ITT-2 ) 1766*1788*1870 
1766 DO 1768 IJJ=1*INTER 

TRANK IJJ)=TRAND< IJJ)*TRANS<IJJi 
1768 CONTINUE 

I F ( I WRITE) 1780 *1870 *1770 
1770 WR ITE(3*1053) 

WR I TE ( 3 * 105 1 ) 

I R ET = 3 
GO TO 1834 

1780 IF( IREAD-1 ) 1782,1785 *2005 
1782 IER=1780 

GO TO 1870 
1785 I TT = I TT+ 1 
GO TO 1310 

1788 DO 1790 IJJ=1*INTER 

TRAND ( IJJ ) = TRAND ( IJJ)*TRANS< IJJ) 

1790 CONTINUE 

I F < I WRITE) 1800* 18 70 *1792 
1792 WRITE<3*1053) 

WR I T E ( 3 * 1 05 1 ) 

I RET =4 
GO TO 1830 

C AVERAGE OVER TWO MESHES******************************************* 

1800 DO 1805 I J J = 2 * INTER 

TRAN( IJJ ) = ( TRANK I JJ)+TRAND( I JJ-1 )+TRAND( IJJ) ) /3 • 

1805 CONTINUE 

WR I TE ( 3 * 1054 ) 

WR I TE ( 3 * 1 05 1 ) 

I R ET = 5 
GO TO 1845 

C WRITE INSTR FOR TRANS MATRIX**************************** 

1830 DO 1831 JI=1,IJJ 
TRAN ( JI ) =TR AND ( J I ) 

1831 CONTINUE 
GO TO 1845 

1832 DO 1833 JI=1,IJJ 
TRAN(JI)=TRANS(JI) 

1833 CONTINUE 
GO TO 1845 

1834 DO 1829 JI=1*IJJ 
TRAN (JI ) =TRAN1 (JI ) 

1829 CONTINUE 

GO TO 1845 
1845 I U= I JJ/5 
IU=IU+1 

DO 1847 I A = 1 * I U 
J J = 1 

DO 1846 11=1*5 
I AA= I A+I U*( II-l ) 

IFURET-5) 1852*1850 *1870 
1850 WAVE I = WA VEO 

WAVEF=WAVEO+< FLOAT ( I NTER ) -1 • ) *5 • 

1852 AA ( JJ ) =WAVE I+FLOAT ( IAA*5)-5* 

JJ=JJ+1 
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AA( JJ) *10000. /AA(JJ-l) 

AA< JJ)=TRAN( IAA) 

JJ*JJ+ 1 

1846 CONTINUE 

WR I TE ( 3 * 1 052 ) ( AA ( JJ> *JJ=1*15 ) 

1847 CONTINUE 

I F ( I RET-2 ) 1504*1763*1848 

1848 I F { I RET— 4 ) 1780 ♦ 1800 * 186 8 
1868 I SE® 1 

C WRITE ON TAPE6 OF TRANSMISSION VALUES FOR A GIVEN LAYER ********** 

TRAN ( 1 )=FLOAT ( I ATM) 

WRITE! 6* 1062) TRAN C 1 ) * ( TRAN ( I J J ) * I JJ = 2 * I NTER ) 

MT*MT~1 

C TEST FOR TOP LAYER IN GEOMETRY MATRIX ****************** 

IF(O-MT) 1095 * 1861 *1861 
1861 WRITE(3tl029) ICASE 
MT«MT+1 

C TEST FOR ANOTHER ANGLE ********************************* 

IF< ICASE-JTAN) 2130 * 2128 * 2 128 
1864 I ST OP= 1864 
REWIND 6 
CALL READ2 
WRITEOt 1058) ISTOP 
STOP 

1870 WR I TE < 3 * 1017 ) IER 
STOP 

C READ INPUT DATA FROM CORE ****************************** 

1900 IF! 1ST >1905*1870*1952 
1905 I D=DWAVE*10 • 

IW«1 

I F ( IFT) 1910*1870*1920 
1910 I WAVES = I WAVE 1 
GO TO 1930 
1920 IWAVES=I WAVE 2 

19 3 0 IF! I WAVE- 1 WAVES) 1932 *1950*1940 
1932 I ER® 1930 

GO TO 1870 

1940 I W AVES = I WAVES+ 1 D 
IW=IW+1 
GO TO 1930 
1950 IWN-! IW— 1 )*5 + l 

IF! I BOS I ) 1952 *1870*1956 
1952 IF! I FT ) 1956*1870*1954 

1954 IFUST) 1955*1870*1956 

1955 I WN® I WN+6000 

1956 DO 1960 JST=1*5 
A ( JST ) =H2 0 ( IWN) 

I WN = I WN+ 1 
1960 CONTINUE 

CALL UNPAK 
DO 1962 JST = 1*5 
XINOT( JST ) -A ( JST ) 

XILINE(JST) =NA ( JST ) 

1962 CONTINUE 

IF! I ST1) 1976*1870*1974 
1974 IWN® I WN— 10 
1976 1ST * l 

I F ( IR) 1530*1870*1410 

1980 IF ( A I J ) 1978*1983*1981 

1978 I ER- 1980 

GO TO 1870 

1981 AIJA®-AI J/(EPSI*EPSI) 

SUM® EXP ( A I J A ) 



T I J = SUM 

XIL1-XILI NE ( J ) 

SUMN1=EXP ( A I JA*XIL1 ) 

IF (SUMN1-1.00) 1585 * 1982*1982 

1982 T I J 1 1 =SUMN 1 

T I JS ( J ) “A I JA 
GO TO 1705 

1983 TIJ11-1. 

T I JS ( J ) *T I J 11 
GO TO 1705 

C READ INPUT DATA FROM T APE***************************************** 

2000 I F ( I WORDS- 13000) 200 2 *2004 *200 A 

2002 I BOS I =— 1 

EPDN= ( WAVE T— 700 • )*( .4) 

I R EAD = 1 
GO TO 2060 

C SINGLE MESH READ IN ************************************ 

2004 I BOS 1 = 1 

I WAVE1 = 7 000 
I W AVE2 = 7 025 
I READ=2 
GO TO 2006 

2005 I READ = I R EAD-1 
I TT = I TT+ 1 

I F ( X SB ) 2001,1870*2003 

2001 SLUMB = 1861 • 

GO TO 2006 

2003 SLUMB = 0 • 

2006 LUMB = 0 
NUMB = 0 
CALL REWI 

I F ( I TEMP ) 2007 *2 008 *2010 

2007 XLUMB=LUMB 

LUMB =X LUMB +7444 • +SLUMB 
GO TO 2028 

2008 XLUMB=LUMB 

LUMB = XLUMB+3722 * + SLUMB 
GO TO 2028 
2010 XLUMB=LUMB 

LUMB=XLUMB+SLUMB 
2028 IF(LUMB) 1870*2032*2030 
2030 CALL RTAW 

NUMB=NUMB+ 1 

IF ( NUMB-LUMB) 2030*2032*1870 
2032 1=1 

EP =0 • 

2034 CALL RTAW 

DO 2036 J = 1 * 5 
H2 0 ( I)=A(J) 

1 = 1 + 1 

2036 CONTINUE 
EP=EP+2 • 

I F ( EP-372 0 • ) 2034*20 38*2038 
2038 ISEC=-1 

GO TO 1310 

C PARTIAL READ IN OF TWO MESHES ************************** 

2 060 IF ( EPDN-EMAX) 2 062 *20 62 *2090 
C START READ IN FROM BEG OF TAPE************************** 

20 62 I F C I TEMP) 2064 *2 066 *2 068 
2064 LUMB=7444 
GO TO 2070 
2066 LUMB=3722 
GO TO 2070 
2068 LUMB = 0 
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2070 IWAVE1=7000 
IWAVE2=7025 
NUMB=0 
CALL REWI 

IF(LUMB) 1870,2074,2072 
2072 CALL RTAW 

NUMB=NUMB+1 

IF (NUMB-LUMB) 2072.2074,1870 
2074 1=1 

EP = 0. 

2076 CALL RTAW 

DO 2078 J = 1.5 
H20( I )=A( J) 

1 = 1 + 1 

2078 CONTINUE 
EP=EP+2 • 

IF ( 2400.-EP ) 2080.2080,2076 
2080 LOWU 1= IWA VE1 

LUPU1=60000+LOWU1 
I TEMP 1 * I T EMP 
CALL REW I 
LUMB=LUMB+1861 
NUMB=0 

2082 CALL RTAW 

NUMB=NUMB+1 

I F { NUMB-LUMB ) 2082 , 2084 , 1 870 

2084 EP=0 

2085 CALL RTAW 

DO 2086 J=1 .5 
H20 1 I )=A( J) 

1 = 1+1 

2086 CONTINUE 
EP=EP+2. 

IF(24 00. -EP) 2088,2088,2085 
2088 L0WS1= IWAVE2 

LUPS1 =60000+LOWS 1 

ISEC=+1 

GO TO 1310 

C START READ IN NOT FROM BEG OF TAPE********************** 

2090 I F ( I TEMP ) 2092,2094,2096 
2092 LUMB= (WAVET-700,— 2,5*EMAX ) /5.+7444. 

GO TO 2097 

2094 LUMB= ( WAVET-700 ,-2 • 5*EMAX ) /5.+3722. 

GO TO 2097 

2096 LUMB= ( WAVET-700 • -2. 5*EMAX ) /5. 

2097 I WAVE1 = ( WAVET-2 • 5*EMAX ) *10 • 

I W AVE2 = I WAVE 1 + 2 5 

NUMB=0 
CALL REW I 

2098 CALL RTAW 
NUMB=NUMB+1 

I F ( NUMB-LUMB ) 2098,2100,1870 
2100 EPl=EPUP+600. 

1 = 1 
EP-0 

2102 CALL RTAW 

DO 2104 J=1 ,5 
H20 ( I J =A ( J ) 

1 = 1 + 1 

2104 CONTINUE 
EP=EP+2. 

2106 IF(2400.-EP) 1870,2108,2102 
2108 L0WU1= IWAVE1 

LUPU1=L0WU 1+6 0000 
ITEMP1-ITEMP 
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LUMB=LUMB+1861 
NUMB=0 
CALL REW I 
2110 CALL RTAW 

NUMB=NUMB+1 

I F ( NUMB-LUMB ) 2110,2112*1870 
2112 EP = 0 
2114 CALL RTAW 

DO 2116 J = l,5 
H20< I )=A ( J) 

1 = 1 + 1 

2116 CONTINUE 
EP=EP+2. 

2118 I F ( 2400 • -EP ) 1870 . 2 120 * 2 1 14 
2120 L0WS1= I WA VE2 

LUPS1= I WAVE 2 +6 0000 

ISEC=+1 

GO TO 1310 

FOR EACH NEW ANGLE OF VIEW REPROCESS TAPE6 TO MAKE NEW TABLE OF 
HEIGHT VS TRANSMISSIONS FOR EACH WAVE NUMBER » AND WRITE ON TAPE6 
2128 LAST = 1 
2130 I GAS = 1 

I CASE= I CASE + 1 
FREQ=WAVEO+DWAVE 
MT= I XB 1 
K 1 =0 
K= 1 

K1=K1+INTER 
REWIND 6 

HEIGHT VS TRANSMISSIONS FOR EACH WAVE NUMBER 

READ ENTIRE TRANSMISSION VALUES FOR A GIVEN ANGLE INTO H20(I) 

DO 2140 J=1 *MT 

READ (6*1 062 ) < H20 ( I ) » I =K »K1 ) 


K=K1+1 


K1=K1+INTER 
2140 CONTINUE 

C PRINT OUT OF HEIGHT VS TRANSMITTANCE 

WR I TE ( 3 * 1066 ) 

WRITE(3»1029)ICASE 

WR I TE ( 3 » 1 01 0 ) (WORDSU) *1 = 1,11) 

I F ( I DOP) 4220,4220,4225 
4220 WR I TE ( 3 , 1 004 ) 

GO TO 4230 
4225 WRITE! 3* 1005) 

4230 CONTINUE 

WR I TE < 3 , 1 000 ) 

WRITE (3, 1002) PR, TEMP, PATH, HET, WAVE I 
1HWIDTH,EMAX,CAI J 
WRITE (3, 1007) CXI ,ZETA,ZETA2 
WRITE (3, 1041) (CLA t J ) ,CLB ( J ) » J=1 ,7 ) , 
WRITE <3* 1042) (CLAlt J) »CLB1 ( J ) , J=1 , 7 
NN 1 1 =0 
NNN 1 1 = 0 
1 1 NTER = I NTER-1 
DO 2330 J = 1 » 1 1 NTER 
I ADD= J*5 
ADD = FLOAT ( I ADD) 

TRAND! J) =WAVEI+ADD 
2330 CONTINUE 
2335 XW1=TRAND ( NN 11+1 ) 

XW 2= TRAND ( NN1 1+2 ) 

XW3=TRAND( NN11+3 ) 

XW 4= TRAND ( NN 1 1+4 ) 


MATRIX ********************** 


,WAVEF ,DNU,DWAVE,PRNOT .TEMPO, 


(CPA! J) ,CPX ( J ) , J=1 ,4 ) 

) , ( CPA ( J ) »CPX ( J ) » J = 1 ,4 ) 
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non 


XW5*TRAND(NNll+5) 

XW6=TRAND ( NN11+6 ) 

XW7=TRAND (NN11+7 ) 

XW8=TRAND (NN11+8 ) 

XW9=TRAND (NN11+9 ) 

XW10=TRAND( NN11+10 ) 

XW11 = TRAND< NN 11+11 ) 

XW12=TRAND( NN11+12 ) 

XW13=TRAND( NN11+13 ) 

XW14=TRAND( NN1 1 +14 ) 

XW15=TRAND( NN11+15 ) 

WRITE! 3*1008) 

WR I TE ( 3 * 2 340 ) XW1 . XW2 »XW3 . XW4 .XW5.XW6.XW7.XW8.XW9.XW10.XW11.XW12.X 
1W13.XW14.XW15 

2340 FORMAT ( 1H . 10X . 15 < 2X .F6.1 ) ) 

2310 DO 2320 J=1.MT 

TRAN ( J ) =GEOM( 1 * J ♦ 1 ) 

NN1 = (MT- J ) * I NTER+NN1 1 
H201=H20(NNl+2 ) 


H202=H20 ( NN1+3 ) 
H203=H20 ( NN1+4 ) 
H204=H20(NNl+5 ) 
H205=H20(NNl+6) 
H206=H20(NNl+7) 
H207=H20 ( NN 1+8 ) 


H208=H20( NN1+9 ) 

H209=H20(NN1+10) 

H2010=H20< NN1+1 1 ) 

H2011=H20(NN1+12 ) 

H2012=H20( NN 1+1 3 ) 

H2013-H2 0 ( NN1+14 ) 

H2014=H20( NN1+1 5 ) 

H2015=H20( NN1+16 ) 

WRITE (3 *2312) TRAN( J ) *H201 .H202 *H203 *H20 A .H2Q5 .H206 .H207 ,H208 , H209 
1.H2010.H2011.H2012.H2013.H2014.H2015 
2312 FORMAT ( 1H . 3X .F7 . 2 . 1 5 ( 2X . F6 .4 ) ) 

2320 CONTINUE 

NNN 1 1 =NNN1 1+1 
NN 1 1 =NNN 11*15 

IF (NN11- INTER 12335. 2142 *2142 
2142 REWIND 6 
KK = 2 

WWAVE I=WAVEI+5* 

2144 CALL RITE1! IXA) 

2145 K=KK 

REARRANGE DATA IN H20(I).AND CONSTRUCT TABLE OF HEIGHT VS 
TRANSM I SS I ON***RECORD SEQUENCE FOR EACH SPECTRAL INTERVAL IS WAVE 
NUMBER .HEIGHT * TRAN .HE I GHT * TRAN*** ******* 

MT= I XB1 
TRAN ( 1 ) =FREQ 
J=2 

DO 2150 1=1.1X81 
TRAN ( J ) =GEOM( l.MT.l) 

J = J + 1 

TRAN ( J ) =H20 ( K ) 


MT =MT-1 
K=K+ I NTER 
J = J+1 

2150 CONTINUE 

C READ TABLE. THUS CONSTRUCTED. INTO TAPE B7 FOR EACH WAVE NUMBER. 

1 1 = I XB+1 
CALL RITEB7I II ) 



i 

i. 


KK«KK+1 

FREQ=FREQ+OWAVE 
IF(KK-INTER) 2145 *2145*2155 
2155 IF ( LAST ) 1090,1090,2160 

2160 I F ( NUM-1 ) 1864 *1864*1080 
END 

830 CARDS 


SIBMAP RITEB7 100 »L I ST * XR 7 *M94 * NODECK 

* RITEB7 WRITES BUFFER TRAN ( 1 1 ) ON TAPE B7. THIS INCLUDES TABLES OF 

* TRANSMISSIONS AND OTHER . PARAMETERS FOR INPUT TO ATMOS3. 


OUT 

FILE 

*B < 3 ) »B( 3 ) *BLK = 3 00*B IN*HOLD * MOUNT ♦ OUT PUT 

RITEB7 

SAVE 

( A*2»l *3*5*6*7 ) I 


CLA* 

3>4 


ALS 

18 


STD 

I I 


TSX 

• OPEN 9 4 


PZE 

OUT 


TSX 

• WR I TE * 4 


PZE 

OUT 

II 

I OCD 

TRANf *300 


RETURN 

RITEB7 

TT 

CONTRL 

TT 


USE 

TT 

TRAN 

BSS 

300 ALL PARAMETERS FOR OUTPUT 


USE 

PREVIOUS 

RITE1 

SAVE 

1*2*3*4*5*6*7*I 


CLA* 

3*4 


ADD 

=1 ADD 1 FOR IXA ITSELF 


ALS 

18 


STD 

JJ 


TSX 

• OPEN * 4 


PZE 

OUT 


TSX 

•WRITE*4 


PZE 

OUT 


IOCP 

BUF 1 * * 20 

JJ 

I OCD 

BUF2 * * 200 


RETURN 

RITE1 

AA 

CONTRL 

AA 


USE 

AA 

BUF 1 

BSS 

2 0 WORDS (11) * WW AVE I * WAVEF * I N TER * I GAS * LA * 


USE 

PREVIOUS HT3 * ANGLE *NUM, JTAN 

ZA 

CONTRL 

ZA 


USE 

ZA 

BUF2 

BSS 

200 I XA*DATA( 200 ) 


USE 

PREVIOUS 

READ2 

SAVE 

1*2*3*4*5*6*7*I 


TSX 

•CLOSE. 4 


PZE 

OUT 


RETURN 

END 

READ2 


43 CARDS 

$ I BFTC KNUMIX M94 ,XR7 , L I ST . REF *NODECX 
SUBROUTINE K.NUM I X ( X I N , Y I N »OUT , I 1 ) 

DIMENSION A(42) ,HH( 10) »XX( 10) 

DIMENSION RA( 32 ) »CA( 32 ) »RB<32 ) » CB < 32 ) » B ( 44 ) » AK ( 5 ) » AM ( 5 ) »DY ( 4 ) 

GO TO (400,401 ) , 1 1 

400 READ (2,710) ( HH ( I ) , I = 1 , 1 0 ) , ( XX ( I ) , I = 1 , 1 0 ) , ( A ( I ) , I = 1 , 42 ) 

710 FORMAT ( 5E14.8/5E14. 8/5F14.8/5F14. 8/ ( 5E14. 8 ) ) 

RETURN 

401 X=X I N 

Y = YIN 


94 



X2 * 

X*X 




Y 2 = 

Y*Y 




IF (X 

-10. ) 

200*201 

.201 

200 

IF (Y 

-1.) 202*202* 

203 

203 

RA ( 1 ) 

= 0. 




CA ( 1 ) 

* 0. 




RB ( 1 ) 

= 1. 




CB ( 1 ) 

* 0. 




R A ( 2 ) 

= X 




CA ( 2 ) 

= Y 




RB ( 2 ) 

= .5- 

X2+Y2 



CB ( 2 ) 

* -2. 

*X*Y 



CB1 = 

CB ( 2 ) 




UV1 = 0 

. 




DO 250 J=2 ♦ 31 
JMINUS = J-l 
JPLUS = J+l 
FLOAT J = JMINUS 
RBI = 2.*FLOATJ+RB(2) 

RA1 = -FLOATJ*(2.*FLOATJ~l.)/2. 

R A ( JPLUS ) =RB1*RA ( J ) -CB 1 *CA < J ) +RA 1*R A ( JM I NUS ) 

CA( JPLUS >=RB1*CA( J ) +CB 1*RA < J > +RA 1*C A (JMINUS) 

RB ( JPLUS ) =RB1*RB ( J ) -CB1*CB ( J ) +RA1*RB ( JMINUS ) 

CB (JPLUS ) =RB1*CB ( J ) +CB1+RB ( J ) +RA1+CB ( JMINUS ) 

UV= ( CA ( JPLUS ) *RB ( JPLUS ) -RA ( JPLUS ) *CB ( JPLUS ) ) / ( RB C JPLUS ) *RB < JPLUS ) + 
1CB ( JPLUS ) *CB ( JPLUS ) ) 

IF ( ABS ( UV-UV1 > -l.E-6 ) 251*250*250 

250 UV1=UV 

251 OUT = UV/ 1 * 772454 
RETURN 

202 IF ( X-2# ) 301*301*302 
301 AINT * 1. 

MAX = 12.+5.*X2 
DO 303 L=1 • MAX 


AJ * MAX-L+1 

303 AINT = AI NT* ( — 2.*X2 ) / ( 2 .*AJ+1 • ) + l • 

U = -2 .*X*AI NT 

GO TO 30 A 

302 IF (X-4.5) 305*306*306 

305 B<43)=0. 

B(44) = 0 . 

J * 42 

DO 307 K = 1*42 
B C J ) = * 4*X*B ( J + l ) -B ( J+2 ) +A ( J ) 

307 J = J-l 

U = B(3)-B( 1) 

GO TO 304 

306 AINT = 1*0 
MAX » 2 .+40 « /X 
AMAX = MAX 

DO 308 K=1*MAX 

AINT = AINT*(2.*AMAX-1#)/<2.*X2)+1« 

308 AMAX * AMAX -1* 

U « -AINT/X 

304 V- 1 .772454*EXP ( — X2 ) 

H = .02 

JM = Y/H 

IF (JM) 310*311*310 


311 H= Y 
310 Z = 0. 

L = 0 

DY ( 1 ) * 0. 

312 DY < 2 ) * H/2 « 
DY ( 3 ) = DY ( 2 ) 
DY ( 4 ) = H 
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318 AK ( 1 ) = 0. 

AM ( 1 ) = 0. 

DO 313 J*l*4 
YY = Z+DYU) 

UU * U+.5*AK(J) 

VV = V+.5*AM(J) 

AKCJ+1) * 2 •* ( YY*UU+X*VV ) *H 
AM ( J + l ) = ~2 •* ( 1 *+X*UU— YY *VV ) *H 
IF (J-3) 313 » 314 * 313 
314 AK ( 4 ) =2 • *AK ( 4 ) 

AM ( 4 ) = AM ( 4 ) +AM ( 4 ) 

313 CONTINUE 
Z=Z+H 
L = L+l 

U * U+* 1666667* ( AK( 2 )+2 **AK { 3 )+AK ( 4 )+AK ( 5 ) ) 

V = V+# 1666667* t AM ( 2 ) +AM ( 3 ) +AM ( 3 ) +AM ( 4 ) +AM (5)1 

IF(JM) 315* 320 * 3 15 

315 IF (L-JM) 318*317*320 

317 AJM = JM 

H * Y-A JM*H 
GO TO 312 

320 OUT = V/ 1 *772454 
RETURN 
201 FI = 0. 

DO 330 J«lflO 

330 F1 = F1+HH< J) / ( Y2+(X“XX( J) )*(X-XX( J) ) )+HH( J) / ( Y2-MX+XXC J ) )*(X+XX ( J) ) 

1) 

OUT ** Y*F 1/3*1415927 

RETURN 

END 

105 CARDS 


SIBMAP 

RTAW 

lOOfNODECK 

RTAW 

SAVE 

1*2*3*4*5*6*7*I 

IN 

FILE 

* A < 4 ) * A ( 4 ) * MOUN 


TSX 

• OPEN * 4 


PZE 

IN 


TSX 

•READ *4 


PZE 

IN 


PZE 

EOFA * ♦ *-2 


I OCT 

A- 5 * ♦ 5 

EOFA 

RETURN 

RTAW 

REWI 

SAVE 

1*2*3*4*5*6*7*I 


TSX 

• CLOSE » 4 


PTW 

IN 


RETURN 

REWI 

UNPAK 

SAVE 

1 


AXT 

5*1 


CAL 

A*1 


ANA 

*0377 


SLW 

NA * 1 


CAL 

A * 1 


ANA 

*0777777777400 


SLW 

A ♦ 1 


T IX 

*-6*1*1 


RETURN 

UNPAK 

AZ 

CONTRL 

AZ 


USE 

AZ 

A 

BES 

5 

NA 

BES 

5 


USE 

PREVIOUS 


END 
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•46224367E 0 .28667551E 0 .10901721E 0 .24810521E-1 .32437733E-2 

» 22833864E-3 .78025565E-5 .10860694E-6 .43993410E-9 .22939360E-12 

.24534071 .73747373 1.2340762 1.7385377 2.254974 

2.7888061 3.3478546 3.944764 4.6036824 5.3874809 

. OOOOOOOOE 0 . 19999999E 0 .OOOOOOOOE 0-0 . 18400000E 0 .OOOOOOOOE 0 

• 15583999E 0 .OOOOOOOOE 0-0. 12166400E 0 .OOOOOOOOE 0 .87708159E-1 

.OOOOOOOOE 0-0.58514124E-1 .OOOOOOOOE 0 .36215730E-1 .OOOOOOOOE 0 

-0.20849765 E-l .OOOOOOOOE-O .11196011E-1 .OOOOOOOOE 0-0 . 5623 1896E-2 
.OOOOOOOOE 0 • 26487634E-2 .OOOOOOOOE 0-0 . 1 1732670E-2 .OOOOOOOOE 0 

. 48995 199E— 3 .OOOOOOOOE 0-0 . 19336308E-3 .OOOOOOOOE 0 .72287745E-4 

. 00000000 EOO-O. 2565 55 12E-4 .OOOOOOOOE 0 .86620736E-5 .OOOOOOOOE 0 

-0.27876379E-5 .OOOOOOOOE 0 .85668736E-6 .OOOOOOOOEOO-0.25184337E-6 

.OOOOOOOOE 0 • 70936022E-7 

00011367600000067500000070000000013000000006300000 
05000100005000100005000100001000050001000050001000050001 
3626837831370665222381031012285418343464525532417966664896028986 
-100000000-0300-0600-0900-0950-1000010000050000000-0300-0600-0900-0950 
0100995030000010005001000300000000010000002000000050-1 
00000 
70. 

EOF E SY 
52 CARDS 
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Appendix E 

Atmos 3 - IBM 7094 Fortran IV Program for Computing the Solution 
to the Radiative Transfer Equation 


SID JOIT V KUNDE OUTGOING ATMOS RADIATION SECT 3-RAD TRANS EON 

SPAUSE 

SDATE 020866 

SEXECUTE IB JOB 

SIB JOB GO .NOLOGIC, NOMAP 

SIBFTC ATMOS3 M94 »XR7 »NOL I ST ♦ NOREF .NODECK 

COMMON/AA/GO»XMINH.PO,RAD»CAPK.*C1.C2.DELR»DELHO»HT,XMAXH.HC,TO 
COMMON/BB /WORDS! 11) . WAVEI .WAVEF. INTER. HT2 ,HT3 , ANGLE ,NUM . IG. I XB 
1,1 WAVE.WAVEO. I8C 

COMMON / CC /HE IGHT » DT AUDH .OLIVER. TR ANU » XX LOG 
COMMON / PP / NOP E » NOPP » SF 
COMMON/RR/PLQT ( 300 ) , IZERO( 125) 

COMMON/SS/RAD1 ( 1200 ) »RAD2 ( 1200 ) ,RAD3 ( 1200 ) »RAD4( 1200 ) ,RAD5(1200 ) » 
1TWAVE ( 15 ) 

2. STO1I350) ,ST02< 350 ) ,ST03 ( 350) .STO4(350) ,ST05(350) ♦ 

3ST 06 ( 350 ) ,ST07( 350) ,ST08( 350) ,ST09 ( 350 ) .ST010 ( 350 ) ,ST01 1(350) . 
5TBW( 350) ,TBWO( 350) ,TRANG( 350) .RATIO! 350) » IST0K350) .IST02C350) .1ST 
603(350) 

COMMON/TT /SUMRAT < 100 ) .ALT I (100) » UMAX IM , SPWORD (11) ,HH< 100) .TEMPI 100 
1),IIXA 

COMMON/ZA/IXA, DATA (200) 

COMMON/ ZB1 /WAVO 1 ,DATB( 200) , I XE , A7M, FT SW , I ATM1 
COMMON/ZB2/WAV02 .DATBB (200 ) ♦ I XF » A8M ♦ GSW, I ATM2 

COMMON/ZCl/WK 11 ) , WAV1 .WAVFl . INTI .HT21 »HT31 .ANGL1 , NUM1 » JTAN1 »LA1 
COMMON/ ZC2/W2 (11) » WAV2 » WAVF2 , I NT2 , HT22 ,HT 32 » ANGL2 ,NUM2 , JTAN2 ,LA2 
COMMON/ ZC 3 / IVGK.IATM 
INTEGER A7M.A8M.FTSW.GSW 

COMMON /DAT AC /TAB 15 ( 1 00 ) /DATAD/TABLE ( 100) 

DATA CA.CB/1.1909000E-12, 1.43879/ 

DIMENSION ALT I T ( 100 ) .SUMLA (100) » SUML(IOO) 

DIMENSION SUMLAY ( 1200) .XH1I1200) 

4 FORMAT (1H , 04X » 35HATMOS CONTRIB, FRACTION OF TOTAL = FS.4) 

5 FORMAT (1H0.05X.14HWAVE NUMBER = F8.3) 

6 FORMAT (10X.11A6) 

7 FORMAT ( 1H0.4X.24HGROUND RADIATION WITH = El 5. 6 » 13X , 29HATM0SPHER I C 
1RADIATI0N WITH = E15.6) 

8 FORMAT ( 1H1 » 7X »4HW AVE .5X .6HLAMBDA » 2X » 8HGRD TRAN) 

9 FORMAT (5X.27HTOTAL RADIATION WITH-OUT = E15 .6 , 10X , 23HTOTAL RAD I AT I 
ION WITH * E15.6) 

10 FORMAT ( 2F6 « 0 ) 

11 FORMAT ( 1H0 » 5I26H WAVE NO. PHI )) 

12 FORMAT ( 1H .10F13.4) 

13 FORMAT (313. 4F 6.2) 

14 FORMAT (1H0///5I26H WAVE NO. EPSILON)) 

15 FORMAT ( 5F6. 1 ) 

16 FORMAT (12) 

17 FORMAT (414) 

26 F0RMATI1H »0PF5.2.F8.5.1PE10.2»0PF8.5,1PE10.3,2X, 

1 0PF5.2.F8.5, 1PE10 . 2.0PF8.5.1PE10.3.2X, 

2 OPF5.2»F8.5»1PE10.2»OPF8.5»1PE10.3,2X ) 

27 FORMAT ( 1H . 3 ( 43HHE I GHT DTAU-DH RADW TRANSU XXLOG )) 

28 FORMAT (1H1.F7.2) 
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29 


FORMAT ( 1H » 5X47H ( W I THOUT WEIGHTING FUNCTION) TOTAL RADIATION = E12 
1 • 5 * 5X28HRAD I AT I ON FROM ATMOSPHERE = E12 • 5 /6X47H (WITH WEIGHTING FUN 
2CTION) TOTAL RADIATION = E 12 . 5 . 5X 28HRAD I AT I ON FROM ATMOSPHERE = 
3 E12.5) 

32 FORMAT (1H1) 

34 FORMAT (1H .12011) 

46 FORMAT ( 1H1.7H . 12X .42H0UTG0ING ATMOSPHERIC RADIATION SECTI 

ION 3/20X .47H0UTG0ING INFRARED RADIANCE. SPHERICAL ATMOSPHERE//) 

48 FORMAT ( 1HO ♦ 19X . 8HANGLE = F6.2.8H DEGREES . 3X . 8HXMI NH = F6.2.12HKM 

1 /20X . 19HUPPER INT. LIMIT = F6.2.12HKM .3X.19HLOWE 

2R INT. LIMIT = F6.2.12HKM ) 

52 FORMAT { 1H . 19X . 46HRADI ANCE = WATT/ (( CM**2 ) *STER*5 (WAVE NUMBERS))// 
1) 

54 FORMAT (1H .34X97HW / 0 WEIGHTING FUNCTION*** 

1WITH WEIGHTING FUNCTION* *//2X4HWAVE ♦ 5X5HLAM 
2DA.4X3HPHI .5X4HEPSI » 7X3HGRD . 1 OX5HATMOS .9X5HTOT AL . 7X2HTB .9X3HGRD. 10 
3X5HATMOS .9X5HTOTAL.7X2HTB) 

56 FORMAT ( 1H .0PF7 . 1 . F8 .3 . F8 . 3 . F8. 3 . 1 PE 14.4 . E14. 4 . E 14.4 ,0PF8 . 2 . 

1 1PE14.4.E14.4.E14.4.0PF8.2 ) 

57 FORMAT ( 1H . F10 . 1 . F 10 .3 . F10 .6 ) 

60 FORMAT (6E12. 5 ) 

61 FORMAT ( 1H . 10X . F 10 . 3 »3X . E 12 . 6 . 3X . F 10 .5 ) 

62 FORMAT ( 1H . 10X . F 10. 3 .3X » E12 . 6 ) 

NOPP NUMBER OF PAIRS OF POINTS IN PHI VS. WAVE NUMBER TABLE (13) 
NOPE=NUMBER OF PAIRS OF POINTS IN EPSILON VS. WAVE NUMBER TABLE (13) 

SF = PHI. AND IS FOUND IN SUBROUTINE PLANK. 

A7M IS TAPE A 7 MULTIPLICITY SWITCH 
A8M IS TAPE A8 MULTIPLICITY SWITCH 
FTSW IS FIRST TIME SWITCH FOR SUBROUTINE READT 
GSW IS GAS MULTIPLICITY SWITCH 
INTER — NUMBER OF INTERVAL 

IXB — NUMBER OF PAIRS OF DATA MULTIPLIED BY TWO FOR HEIGHT VS TRANS 
MISSION TABLE 

IG — NUMBER OF ANGLES OF VIEW 
NUM — NUMBER OF ATMOSPHERIC MODELS 
IBC = 1 
FTSW = 0 
XMINH=0. 

I VGK=— 1 

READ (2.17) NOPP. NOPE .A7M.A8M 
C READ IN OF SPECTRAL RESPONSE***************************** 

read (2.6) Tspword< r ) . 1 = 1 . 11 ) 

I A a 2*N0PP 

READ (2.10) ( T AB 15 ( I) . 1 = 1. 1 A) 

WR I TE ( 3 * 32 ) 

WR I TE ( 3 . 1 1 ) 

WRITE! 3.12) ( TAB 15 ( I ) .1=1* I A) 

C READ IN OF SPECTRAL EM ISS I V I TY*************************** 

I A =2*NOPE 

READ (2*10) ( TABLE (I).I=1»IA) 

WRITE! 3. 14) 

WR I TE ( 3 * 12 ) (TABLE ( I ) . 1 = 1, IA) 

C GSW=GAS SWITCH (0.1 =A7 ONLY , 2.+ = A7 AND A8 ) ********************* 
READ ( 2 ♦ 13 ) I ATM1 . I ATM2 .GSW.TM INI .TMAX 1 .TMIN2 .TMAX2 
C TWAVE(I) — WAVE NUMBER FOR PRINT OUT OF VERTICAL PARAMETERS*** 

READ (2.15) (TWAVE(I). 1=1.15) 

NNUM=0 
DELR=. 1 
19 I ANG=0 
21 JTA=0 

DO 50 1=1.125 
I ZERO ( I ) =0 
50 CONTINUE 

DO 1003 1=1.1200 


100 


1003 SUMLAYd )=0.0 
I A* 1 

I WAVET = TWAVE ( IA)*10. 

WR I TE ( 3 * 32 ) 

C MODEL ATMOSPHERE TABLE FROM TAPE ( S ) ********************************* 
CALL READF 
GO TO 39 
22 CALL READT 
39 I F ( I WAVE— I WAVET ) 42>40»42 
AO IWP=1 

I A= I A+l 

I WAVET = TWAVE ( IA)*10, 

GO TO 44 
42 IWP=-1 
44 J=1 

C TOP LAYER************************************************************ 
RUPW=0. 

RUPWO=0* 

DTAUDH = 0*0 
DW AVE=5 • 

H=HT3 

DELRl=DELR/2.0 

C ATMOSPHER I C LAYER S**** «■*******##********•*■**********■**-»****#******#** 

626 H1=H-DELR 1 
XH1 t IAA) “HI 
H=H-DELR 

CALL INTB(H.TRAN) 

I F ( TRAN ) 627 .628 » 628 

627 TRAN=0.0 

628 TRANL'TRAN 

CALL INTA ( HI *TH1 ) 

CALL PLANK t WAVEO .DWAVE *TH1 .RADW.RADWO ) 

HEIGHT =H1 

C TRANSFER EQUATION**************************************** 

OL IVER=( T RANU-TRANL ) *RADW 
RUPW=RUPW + OLIVER 
RUPWO = RUPWO+( TRANU-T RANL ) *RADWO 
DTAUDH = ( TRANU-T RANL > /DELR 

XXL0G=7. 7 

IF(IWP) 665 » 722 * 6A9 
6A9 CALL STORAD ( IAA ) 

C SUMMATION OF RADIANCE FOR EACH LAYER FOR ENTIRE SPECTRUM************ 
665 OLIVER' ( TRANU-TRANL ) *RADW 

SUMLAY ( IAA) “SUMLAYt IAA) + OLIVER 
TRANU'TRANL 
LAYTOT “IAA 
I AA= I AA+ 1 

I F t H-HT2 J708.708.626 
708 JTA= JTA+ 1 

C GROUND RADIATION***************************************************** 
CALL INTA ( H »TH1 ) 

CALL PLANK < WAVEO .DWAVE .TH 1 .RADW »R ADWO ) 

RW=RUPW 

RWO'RUPWO 

CALL INTRPDIWAVEO.EPSI .NOPE) 

RADGW=TRANU*EPSI*RADW 
RADGWO“TRANU*EPSI *R ADWO 
TRANG ( JTA ) “TRANU 

C TOTAL UPWARD RADIATION********************************************** 
RUPW=RUPW+RADGW 
RUPWO=RUPWO+RADGWO 
IF(IWP) 709.722.800 


709 PLOT ( JTA ) =RUPWO 
STOl ( JTA ) =WAVE0 

ST02( JTA 1=10000. /STOl(JTA) 

ST03 ( JTA ) = SF 
ST 04 (JTA) « EPSI 
ST06 ( JTA ) =RADGWO 
ST07(JTA)=RW0 
ST08 ( JTA ) =RUPWO 
ST09 (JTA 1 =RADGW 
ST010 ( JTA ) =RW 
ST011 ( JTA)=RUPW 

TBW( JTA) =CB*WAVEO/ALOG< 5.0*CA*WAVE0**3/RUPW +1.0) 

TBWO( JTA )=CB*WAVE0/AL0G( 5. 0*CA*WAVE0**3 /RUPWO +1.0) 

C TEST FOR ANOTHER WAVE NUMBER**************************************** 
WRITE(3*63) WAVEO 
63 FORMAT ( 1H .F15.2) 

IF ( JTA-INTER) 22.710*710 

C PLOT OF UPWARD RADIATION******************************************** 

710 PLOT 1 =PLOT ( 1 ) 

DO 715 1=1. JTA 

PLOT ( I ) =PLOT ( I ) /PLOT 1 
PLOT ( I ) =PLOT ( I >*100. 

715 CONTINUE 
WRITE<3.32) 

CALL PWD ( 4 ) 

DO 720 1=1. JTA 
I PLOT =PLOT ( I ) 

I PLOT = I PLOT + 1 
I F ( I PLOT-120 ) 718.718.716 

716 I PLOT = 120 

718 WR I TE ( 3 » 34 ) ( I ZERO ( J ) . J = 1 . IPLOT ) 

720 CONTINUE 
GO TO 850 

C TEST FOR ANOTHER ANGLE********************************************** 

721 IANG=IANG+1 
IF(IANG-IG) 21.724.724 

C TEST FOR ANOTHER ATMOSPHERIC MODEL********************************** 
724 NNUM=NNUM+1 
IVGK=-1 

I F ( NNUM-NUM ) 19.722.722 
722 CALL READ2 

CALL MESAGE ( 9 ) 

C PRINT OUT OF VERTICAL PARAMETERS FOR GIVEN WAVE NUMBER*************** 
800 I AB 1 =0 
I AB2 =0 
I AB= I AA 

810 WRITE13. 28) WAVEO 
WRITEI3.27) 

IF(IAB-180) 830.820.820 
820 I AB 1 = I AB2 + 1 
I AB2=lAB2+60 

WRITE (3. 26) (RADI (IB) »RAD2 ( IB) .RAD3( IB) .RAD4( IB) .RAD5( IB) .RADI ( I B+ 
160 ) .RAD2 ( IB+60) .RAD3 ( IB+60) »RAD4( I B+60 ) »RAD5( IB+60) .RADI ( IB+120) . 
2RAD2 (IB+120 ) .RAD3 ( IB+120) .RAD4 ( IB+120) .RAD5 ( IB+120) , I B= I AB1 . I AB2 ) 

I AB2= IAB2+120 
I AB= I AB- 180 
GO TO 810 

830 II=FLOAT ( IABJ/3.+1. 

IAA1 = I I 
I AA2 = I I + I I 
I AB1 = I AB2 + 1 
I B= I AB1 
I Bl= IB+I AA1 
I B2= I B+ I AA2 
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DO 840 1*1 t II 

WR I TE ( 3 .26 ) RADI ( IB) *RAD2 < IB) .RAD3f IB) .RAD4! IB) .RAD5I IB) * 

1RAD1! 181 ) »RAD2 ( IB1 ) ♦ RAD3 ( I B1 ) ,RAD4( IB1) .RAD5! IB1 ) ♦ 

2RAD1 t IB2 ) *RAD2 ( 1B2 ) » RAD3 ( IB2 ) »RAD4( IB2) »RAD5( IB2 ) 

IB* IB+1 
IB1=IB1+1 
IB2= IB2+ 1 
840 CONTINUE 
K=1 

1 1 =LAYTOT— 9 
1 1 I=LAYTOT 

841 SUMLA1=0 • 0 

DO 842 I*II.III 

842 SUMLA1=SUMLA1 + RAD3 (I) 

ALTITU=XH1< II+5) + DELR1 
ALT IT ( K) *ALTI TU 
SUMLA(K) *SUMLA1 

K*K+1 
11*11-10 
1 1 1 = 111-10 
IF ( 1 1 ) 8,43.843.841 

843 WRITEI3.32) 

KMAX IM=K— 1 

DO 845 1*1. KMAX IM 
JJJJJJ=KMAX IM— 1+1 
ALT I ( I ) =ALT I T ( JJJJJ J ) 

SUML ( I ) =SUMLA ( JJJJJJ ) 

SUMRAT ( I ) =SUML ( I ) /RUPW 
845 WR I TE ( 3 . 6 1 ) ALT I (I) . SUML ( I ) , SUMRAT ( I ) 

WTMAX=0.05 

WTMIN=0.0 

CALL OCPLOT IWTMAX .WTMI N .SUMRAT . 3 ) 

WRITE (3. 5) WAVEO 
GO TO 709 

C PRINT OUT OF ALL SPECTRAL INTERVALS********************************** 
850 WR I TE ( 3 » 46 ) 

CALL PWD ( 4 ) 

WRITE! 3. 6) (SPWORD! I ) . 1 = 1 .11 ) 

WR I TE ( 3 » 48 ) ANGLE »XM I NH ,HT 3 *HT2 
WR I TE ( 3 . 17 ) INTER 
WR I TE ( 3 ♦ 52 ) 

IAB1=0 
I AB2=0 
I AB= I NTER 
860 WR I TE ( 3 » 54 ) 

IF ( IAB-47 ) 880.870.870 
870 I A81= IAB2+1 

I AB2 = I AB2 + 47 

WRITE (3. 56) ( STOl (IB) .ST02 ( IB) .ST03J IB) .ST04 ( I B ) . ST06 ( IB) .ST07( IB) 
1 »ST08 ( IB) » TBWO ( IB) .ST091 IB) .ST010( IB) .STOll ( IB) »TBW( IB) .IB = IAB1. 
2IAB2 ) 

IAB = IAB - 47 
872 WR I TE ( 3 . 32 ) 

WR I TE ( 3 ♦ 54) 

IF ( I AB-5 5 ) 880.874.874 
874 IAB1=IAB2+1 
I AB2= I A82+55 

WRITE! 3.56) (STOl! IB) .ST02 ( IB) .ST03! IB) .ST04! IB) »ST06( IB) .ST07 (IB) 
1.ST08! IB) .TBWO! IB ) . ST09 l IB ) .ST010 ( IB) .STOll ( IB) *TBW( IB) .IB=IAB1. 
2IAB2) 

I AB= IAB— 55 
GO TO 872 
880 IAB1=IAB2+1 
IB=I AB1 
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DO 890 1=1*1 AB 

WRITE <3* 56) (ST01( IBJ .ST02 ( IB) »ST03( IB) *ST04< IB) *ST06( IB ) * ST 07 < IB) 
1.ST08I IB) »TBWO( IB) *ST09( IB) *ST010 ( IB) *ST0X1 ( IB) »TBW( IB) ) 

IB= IB+1 
890 CONTINUE 


SSUMW=0.0 

SSUMWO=0.0 


SUMATM=0.0 
SUMGRD=0.0 
DO 925 1=1* INTER 
SUMATM=SUMATM+ST010< I ) 

SUMGRD=SUMGRD+ST09( I ) 

SSUMWO=SSUMWO + ST08 ( I ) 

925 SSUMW=SSUMW + STOll(I) 

A TMCON =SUMA TM/SSUMW 
WR I TE ( 3 * 7 ) SUMGRD*SUMATM 
WR I TE ( 3 ♦ 9 ) SSUMWO . SSUMW 
WRITE (3*4) ATMCON 
WRITEI3.8) 

WRITE (3. 57) ( ST01 ( I ) *ST02< I ) . TRANG ( I ) *1=1. INTER) 


K= 1 

1 1 =LAYT0T-9 
1 1 I=LAYTOT 

1009 SUMLA1=0 « 0 

DO 1010 I=II*III 

1010 SUML A 1 =SUMLA 1 + SUMLAY(I) 
ALT I TU=XH1 ( 1 1 + 5 ) + DELR1 
ALT I T ( K) =ALTITU 

SUMLA ( K ) =SUMLA 1 
K=K + 1 


11=11-10 

111=111-10 

IF ( 1 1 ) 1021*1021.1009 
1021 WR I TE ( 3 * 32 ) 

KMAXIM=K-1 
DO 1033 1=1* UMAX I M 
JJJJJJ=KMAXIM-I+1 
ALTI ( I ) =ALT I T ( JJJJJJ ) 

SUML ( I ) =SUMLA ( JJJJJJ ) 

SUMRAT ( I ) =SUML ( I ) /SSUMW 
1033 WR ITE(3*61) ALT 1(1) *SUML t I ) *SUMRAT t I ) 

CALL EOF A3 

C PLOT WAVE NUMBER VS TBW********************************************* 
CALL OCPLOT(TMAXl.TMINl*TBWO*l) 

TRMAX =1.0 
TRMIN = 0.0 

C PLOT WAVE NUMBER VS TRANG******************************************* 
CALL OCPLOT ( TRMAX ♦ TRMI N . TRANG *2 ) 

C PLOT WEIGHTING FUNCTIONS******************************************** 


C 


WTMAX=0.05 

WTMIN=0.0 

CALL OCPLOT (WTMAX .WTMI N .SUMRAT *3) 

WRI TE ( 3 ♦ 7 ) SUMGRD.SUMATM 
WR I TE ( 3 * 9 ) SSUMWO* SSUMW 
WRITE (3. A) ATMCON 

PLOT TEMPERATURE VS HEIGHT****************************************** 


TMAX3=330 • 
TMIN3=130. 

I IXA=IXA/2 
DO 1041 N=1 .II XA 


JI J=IXA-N*2+2 
JI JI=IXA-N*2+1 
TEMP(N) = DATA( JI J) 
HH(N)=DATA( JI JI ) 
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1041 CONTINUE 

CALL OCPLOTITMAX3.TMIN3.TEMP.4) 

C PROGRAM IS FINISHED AT 722****************************************** 

GO TO 721 
END 

377 CARDS 


SIBFTC OCPLOT M94.XR7 »NOL 1ST . NODECK »NOREF 
SUBROUT I NE OCPLOT ( XMAX » XM I N * VS » NAME ) 

CONSTANT NAME = 1(IS WAVE NO. VS TBB).=2(IS WAVE NO. VS TRANG ) *********** 
DIMENSION NSCALE t 5 ) 

COMMON/BB/WORDS 1 11 ) . WAVE I .WAVEF * I NT ER .HT2 »HT3 * ANGLE *NUM » I G. I XB 
1*1 WAVE tWAVEO * I BC 
COMMON /OP /TMAX.TM IN 
COMMON/RR/PLOT ( 300 ) *IZERO(125) 

COMMON/SS/A t 6000 ) * 
lTWAVEt 15 ) 

2* STOH350) *ST02( 350 ) *ST03( 350) .ST04 < 350 ) *ST05 ( 350 ) * 

3ST06 t 350 ) *ST07 ( 350 ) .ST08I350) »ST09(350> .ST010 < 350 > .ST01 1 ( 350 > . 

5TBW ( 350 ) *TBW0( 350) . TRANG( 350) *RAT I Ot 350) . ISTOl t 350 ) ♦ I ST02 ( 350) ♦ 1ST 
603(350) 

COMMON/TT /SUMRAT ( 100 ) ♦ ALT I ( 100 > . KMAX I M .SPWORD < 11 ) *HH( 100) *TEMP( 100 
1 ) * I IXA 

COMMON /ZB2 /WAV02 *DATBB(200)*IXF*A8M*GSW*IATM2 

COMMON/ZC1/W1 ( 11 ) * WAV 1 *WAVF1 »INT1 .HT21 *HT31 *ANGL1 *NUM1 ♦ JTAN1 .LAI 
COMMON /ZC2/W2 ( 11 ) . WAV2 . WAVF2 . I NT2 * HT22 *HT32 * ANGL2 *NUM2 * JTAN2 *LA2 
INTEGER A7M.A8M.FTSW.GSW 

DATA BCD/0546060606060/. ( NSCALE ( I ) . I = 1 ♦ 5 ) /32767 .0 . 2 .0 . 3/ 

DATA HMIN.HMAX/O. 0.70.0/ 

DATA W0MIN.W0MAX/500. 0.2000.0/ 

WRITE (3.30) 

30 FORMAT (1H1) 

GO TO (50.55). GSW 
50 WR I TE ( 3 ♦ 1 0 ) <W1 ( I ) .1 =1.11 ) 

WRITEI3.10) (SPWORD( I). 1=1. 11) 

10 FORMAT (30X.11A6) 

GO TO 60 

55 WRITE13.10) (Will). 1=1. 11) 

WRITEI3.10) (W2( I ) .1=1.11) 

WRITEI3.10) (SPWORDI I ) ♦ 1=1 .11 ) 

60 CONTINUE 

GO TO ( 200.200.210.220.205 ) .NAME 
200 CALL PLOT 1 < NSCALE .12.25.10.10) 

C THE 2ND ARG IS NO. OF BLOCKS DESIRED ALONG Y 

C THE 3RD ARG IS NO. OF POINTS DESIRED IN EACH BLOCK ALONG Y 

C THE 4TH ARG IS NO. OF BLOCKS DESIRED ALONG X 

C THE 5TH ARG IS NO. OF SUB-INTERVALS DESIRED IN EACH BLOCK ALONG X 

CALL PLOT2I A .XMAX .XM I N . WOMAX . WOMI N ) 

CALL PLOT 3 ( BCD .VS.ST01.INTER) 

C THE 1ST ARG IS THE BCD CHARACTER DESIRED TO REPRESENT THE POINTS 

CALL FPLOT4I52.52HW AVE NUMBE 

1 R ) 

GO TO (201.202.203.204.205 ) .NAME 

201 WRITEI3.21) 

21 FORMAT!// 42X.52HT EMPERATUR 
1 E ) 

RETURN 

202 WR I TE ( 3 . 22 ) 

22 FORMAT!// 42X.52HG ROUND TRANSMISSION VALU 
2E S ) 

203 RETURN 

210 CALL PLOT 1 ( NSCALE .14.10*10.10) 

CALL PLOT 2IA.XMAX.XMIN. HMAX *HMI N ) 
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CALL PL0T3 ( BCD » VS. ALTI .KMAXIM) 

CALL FPL0T4 ( 32 ♦ 52H HEIGHT*KM 

1 ) 

WRITE (3. 26) 

26 FORMAT!// 42X . 52HFRACT IONAL OUTGOING EFFECTIVE RADIANCE / KM 
1 ) 

RETURN 

220 CALL PLOTKNSCALE.14.10.10.10 ) 

CALL PLOT 2 ( AtXMAX *XM IN ♦ HMAX *HMIN ) 

CALL PLOT 3 ( BCD*VS*HH .IIXA) 

CALL FPLOT 4(52* 52H HEIGHT*KM 

1 ) 

WR I TE ( 3 *21 ) 

204 RETURN 

205 RETURN 
END 

72 CARDS 


SIBFTC PLANK M94 *XR7 *NOL 1ST .NODECK .NOREF 

SUBROUT I NE PLANK ( WAVEO * DWAVE ♦ TH1 . R ADW . RADWO ) 

DIMENSION PLK ( 5 ) 

COMMON /PP /NOPE .NOPP *SF 

C SF = PHI. AND IS FOUND IN SUBROUTINE PLANK 
DATA CA .CB/1. 1909000E-12 . 1.43879/ 

RADW=0. 

RADWO=0. 

DW=WAVE0-2.0*DWAVE/5.0 
DO 21 1=1.5 

PLK ( I ) = CA*DW**3/(EXP(CB*DW/TH1)-1.0) 

17 CALL INTRPC(DW.SF.NOPP) 

RADW = RADW+PLK < I ) *SF 

18 RADWO=RADWO+PLK ( I ) 

DW=DW+DWAVE/5. 

21 CONTINUE 
RETURN 
END 

SIBFTC STORAD M94.XR7 .NOL 1ST .NODECK. NOREF 
SUBROUTINE STORAD ( IAA) 

C THIS SUBROUTINE STORES HEIGHT .DTAUDH. R 1PW . TRANU .RUPWO IN RADX(IAA). 

C WHERE X=1 5 . 

COMMON/CC /HE IGHT. DTAUDH. OLIVER. TRANU .XXLOG 

COMMON/ SS/ RADI ( 1200 ) *RAD2( 1200) .RAD3 ( 1200) *RAD4( 1200) .RAD5(1200) . 

1 TWAVE ( 15 ) 

2. STOl ( 350 ) .ST02I 350 ) .ST03 (350) .ST04( 350 ) *ST05( 350) . 

3ST06I 350 ) ,ST07< 350) .ST08( 350) .ST09 ( 3 50 ) .ST010 ( 350 ) .STOll ( 350) . 

5TBW( 350 ) . TBWO ( 350) .TRANG< 350) .RATIOt 350 ) . ISTOl ( 350 ) . I ST02 ( 350 ) . I ST 
603 ( 350 ) 

RADI ( IAA ) =HEIGHT 
RAD2 ( I AA ) = DTAUDH 
RAD3 ( IAA ) =OLI VER 
RAD4 ( IAA ) =TRANU 
RAD5 ( IAA ) =XXLOG 
RETURN 
END 

36 CARDS 


SIBFTC MATMUL M94 .XR7 .NOL I ST . NODECK . NOREF 
SUBROUTINE MATMUL (NOOM) 

C SUBROUTINE MATMUL MULTIPLIES MATRICES DATB AND DATBB . OR EXCHANGES 
C DATBB FOR DATB ACCORDING TO WHETHER NOOM IS .EQ. 1 OR 2 • 

C ALL ARGUMENTS( EXCEPT NOOM) ARE TRANSMITTED THROUGH COMMON STO-AGE 
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C AREAS. SUBROUTINE REAOT CALLS MATMUL AND DIRECTS ITS USAGE. 

COMMON/BB/WORDS (11). WAVE I .WAVEF. I NTER .HT2 .HT3 .ANGLE .NUM * IG. I XB 
1*1 WAVE .WAVEO * I BC 
COMMON /Z A/ 1 XA * DATA ( 200 ) 

COMMON/ZB 1/WAV01 *DATB(200) » I XE. A7M. FTSW . I ATM1 
COMMON/ZB2/WAV02 .DATBB ( 200 ) , I XF . A8M *GSW . I ATM2 

C0MM0N/ZC1/W1 ( 11 ) , WAV1.WAVF1 .INT1.HT21 *HT 31 * ANGL 1 *NUM1 » JTAN1 * LAI 
COMMON/ZC2/W2 ( 11 ) * WAV2 .WAVF2 * INT2 *HT22 *HT 32 * ANGL 2 * NUM 2 * JTAN2 *LA2 
COMMON/ZC3/ I VGK * I ATM 
COMMON/ TS/T RAM ( 300 ) 

INTEGER A7M.A8M.FTSW.GSW 
DATA C1/6HMUTAPE/ 

C A7M IS TAPE A7 MULTIPLICITY SWITCH 
C A8M IS TAPE A8 MULTIPLICITY SWITCH 
C FTSW IS FIRST TIME SWITCH FOR SUBROUTINE READT 
C GSW IS GAS MULTIPLICITY SWITCH 

C NOOM IS THE NO. OF OPERATING MODE ( ^MULTIPLICATION OF DATB*DATBB 
C AND 2=EXCHANGE DATBB FOR DATB ) 

GO TO ( 90 * 250 ) *IBC 
90 GO TO (100*200) *GSW 
100 ANGLE=ANGL1 
WAVE I =WAV1 
WAVEF=WAVF 1 
HT2=HT2 1 
HT 3 = HT 3 1 
INTER=INT1-1 
NUM=NUM1 
IG“JTAN1 
IBC»2 

WRITE (3* 500) WAV1»WAVF1.INT1.HT21*HT31*ANGL1»NUM1*JTAN1*LA1 
WRITE! 3.500) WAVE I * WAVEF . I NTER . HT2 . HT3 . ANGLE »NUM , I G » I VGK. 

5 00 FORMAT ( 1H0 . 10X *2F 10 . 1 . I 5 . 3F10 .2 . 3 1 5 ) 

GO TO 250 

200 INT78 = (WAVF2-WAV1 ) /5. + 1. 

I F ( INTI— 1-INT78) 11*12*12 

11 INTER=INT78 
WAVE I =WAV1 
WAVEF=WAVF2 
GO TO 13 

12 INTER® INT 1-1 
WAVE I = WAV 1 
WAVEF = WAVF1 

13 NUM=NUM1 
IG = JTAN1 
ANGLE = ANGL1 
HT 3 = HT 3 1 

HT 2 = HT 2 1 
I BC®2 

WRITE (3* 500) WAV1*WAVF1»INT1.HT21.HT31*ANGL1 *NUM1 *JTANl»LAl 
WRITE( 3.500 ) WAV2.WAVF2*INT2*HT22*HT32.ANGL2.NUM2. JTAN2.LA2 
WRITE (3. 500) WAVEI.WAVEF. INTER. HT2.HT3. ANGLE. NUM * I G* I VGK 
250 CONTINUE 

IXB = LA1*2 
IX = IXB/2 
J = 2 

GO TO (1.2.3) .NOOM 

C TWO GASES*********************************************************** 

1 DO 10 1 = 1. IX 

DATB ( J ) =DATB ( J ) *DATBB ( J ) 

10 J= J+2 

C ONE GAS .A7********************************************************** 

3 I WAVE=WAV01*10 . 0 
I ATM=ATM 1 
WAVEO = WAV01 
GO TO A 
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C ONE GAS»A8*«******»******** a*************************************** 

2 DO 20 I=1*IX 

DATBt J)=DATBB< J) 

20 J = J+2 

IWAVE=WAV02*10.0 
I ATM=ATM2 
WAVEO = WAV02 
4 GO TO (60,30) *GSW 

30 IF(IVGK) 31*31,35 

31 IVGK=+1 
W1(11)=C1 
INTER=INTER+1 
CALL RITE1( IXA) 

INTER=INTER-1 

WR I TE ( 3 , 6280 ) WAVEO 

6280 FORMAT ( F 1 0,3 ) 

WR I TE ( 3 * 628 1 ) (DATB(J) *J=1*IXB) 

6281 FORMAT ( 1 0F10, 5 ) 

35 TRAM ( 1 ) =WAVE0 

I FB= I XB+ 1 
DO 37 1=2, IFB 
37 TRAM ( I ) =DAT8 ( I -1 ) 

CALL RITEB7 ( IFB) 

60 RETURN 
END 


95 CARDS 


5IBFTC MESAGE M94 *XR7 *NOL 1ST * NODECK! * NOR EF 
SUBROUTINE MESAGE (NOM) 

COMMON/ZBl/WAVOl, DATS (200) , I XE , A7M , FTSW , IATM1 
COMMON /ZB2 /WAV02 ,DATBB (200 ) , I XF , A8M ,GSW , I AT M2 
INTEGER A7M , A8M , FTSW *GSW 
C A7M IS TAPE A7 MULTIPLICITY SWITCH 
C A8M IS TAPE A8 MULTIPLICITY SWITCH 
C FTSW IS FIRST TIME SWITCH FOR SUBROUTINE READT 
C GSW IS GAS MULTIPLICITY SWITCH 

C SUBROUTINE MESAGE PRINTS ON LINE ERRORS, TAPE HANDLING AND SNAFUES 
C THAT OCCURS DURING THE READING OF TAPES A7 AND A8. MOST OF THE 
C MESSAGES ORIGINATE IN SUBROUTINE READT 

C NOM IS THE NUMBER OF THE MESSAGE THAT IS TO BE PRINTED 

C 

GO TO (1, 2, 3, 4, 5, 6, 7, 8, 9) , NOM 

1 PRINT 10 

10 FORMAT (1H0.30X50HSWI TCH TAPES A7 AND A8 PRESS ST ART ( AC HAS HPR 1 

1)//////////) 

PAUSE 1 
RETURN 

2 REWIND 8 

21 REWIND 7 

22 PRINT 20 

20 FORMAT ( 1H0 , 30X89HTHI S MESS IS IMPOSSIBLE, DATA IS OUT OF FOCUS 

1STOP — (AC HAS HPR 2), PRESS START TO EXIT//////////) 

PAUSE 2 
CALL EXIT 

3 PRINT 30 

30 FORMAT ( 1HO, 30X53HWAVEQ1 IS .GT. WAVE02*NO MORE CORRELATION IS POS 
1SIBLE) 

GO TO 2 

4 REWIND 7 
PRINT 40 

40 FORMAT ( 1H0 , 30X48HM0UNT NEXT A7 TAPE ON TAPE UNIT A7(AC HAS HPR 4)/ 
1 //////////) 

PAUSE 4 
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return 

5 REWIND 8 
PRINT 50 

50 FORMAT < 1H0 ♦ 30X48HMOUNT NEXT A8 TAPE ON TAPE UNIT A8IAC HAS HPR 5)/ 

1 / //////// ) 

PAUSE 5 
RETURN 

6 PRINT 60 

60 FORMAT ( 1H0 • 30X38HTAPE A7 HAS EOF INI T I ALLY *STOP PROGRAM) 

GO TO 21 

7 PR I NT 70 

70 FORMAT ( 1H0.30X54HTAPES (A7.A8) ARE OUT OF DATA TABLES AT THE SAME 
1 TIME) 

GO TO 2 

8 PRINT 80 

80 FORMAT ( 1H0 » 30X38HTAPE B7 HAS EOF I N I T I ALLY * STOP PROGRAM) 

GO TO 2 
9 PRINT 90 

90 FORMAT (1H0.30X48HTHIS RUN IS FINISHED* NOTE NEXT INSTRUCTION ) 

GO TO (21.2) .GSW 
END 

57 CARDS 


SIBFTC PWD M94.XR7. NOLIST. NODECK. NOREF 

SUBROUTINE PWD(LL) 

COMMON/ZA/ 1 XA * DATA ( 200 ) 

COMMON/ZB1 /WAVO 1.DATBI200 ) * I X E . A7M . FTSW ♦ I ATM1 
COMMON/ZB2 / WAV02 .DATBB (200 ) * I XF.A8M.GSW. I AT M2 

COMMON/ZC1/W1 ( 11 ) * WAV1 * WAVF 1 *INT1*HT21 .HT31 * ANGL1 .NUM1 . JTAN1 .LAI 
COMMON / Z C2 / W2 ( 1 1 ) * WAV2 .WAVF2 . INT2.HT22 .HT32 * ANGL2 .NUM2 . JTAN2 . LA2 
COMMON/ZC3/ IVGK . I ATM 

6 FORMAT ( 1H0.3X.11HHT BOTLAY = F9.5.3X. 20HATMOSPHERE NUMBER = I2.3X.2 
11HHEIGHT OF TOP LAYER = F6.2.12HICM /4X.15HANGLE OF VIEW = F 

25.2.8H DEGREES. 3X.20HWAVE NUMBER RANGE =<F6.1.3H TOF6.1.6H) CM-1 ) 

10 FORMAT l 9X * 1 1 A6 ) 

11 FORMAT (5 ( 3X.F6.2.2X.F8.2 ) ) 

GO TO (100.200.100.400) *LL 

100 WR I TE ( 3 * 1 0 ) (W1(I ) *1=1*11) 

WRITEI3.il) ( DATA ( I ) * I = 1 * I XA ) 

WR I TE ( 3 * 6 ) HT2 1 . IATM1.HT31 . ANGL 1.WAV1.WAVF1 
GO TO (300.300.200) ♦ LL 
200 WRITE13.10) (W2(I).I =1.11) 

WR I TE ( 3 . 1 1 ) ( DATBB ( I ) . I =1.IXA> 

WRITEI3.6) HT22 * I ATM2 . HT32 . ANGL2 . WAV2 . WAVF2 
300 RETURN 

400 WRI TE ( 3 . 10 ) (W1(I). 1 = 1.11) 

RETURN 

END 

25 CARDS 


$ I BMAP EOFA3 10 .NODECIC .M94 »NOL I ST .NOREF 

OUTPT FILE .OU1 .OU2 .READY .OUTPUT .BLX*25 
EOFA3 SAVE ( 4 .2 . 1 .3 .5 . 6 .7 ) I 

TSX .OPEN. 4 

MZE OUTPT 

TSX .CLOSE. 4 

MZE OUTPT 

RET RETURN EOFA3 
END 

9 CARDS 
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SIBMAP REACT 200 .NOL I ST * XR7 »M94 * NODECK . NOREF 

* READT SUBROUTINE READS TAPE A7 * AND TAPE A8 IF REQUIRED 

* CHECKS FOR INITIAL WAVE NUMBER CORRELATION BETWEEN A7 AND 

* A8.APPR0PIATELY CALLS FOR MATMUL (MATRIX MULTIPLICATION) 

* OF DATB X DATBB AND DIRECTS THE PROPER MAINTENANCE OF 

* DATB AT ALL TIMES. ALL PARAMETERS ARE TRANSMITTED 

* THROUGH COMMON //. DATB IS UPDATED ONCE DURING EACH 

* CALL READT 


INP1 

FILE 

* A ( 3 ) * A ( 3 ) 

*8LK»300*BIN*HOLD *MOUNT 

I NP2 

FILE 

* A (4 ) * A ( 4 ) 

.BLK=300. BIN. HOLD. DEFER 

READT 

SAVE 

14*2*1*3*5 

*6*7)1 


CLA 

GSW 

ARE THERE TWO GASES 


CAS 

= 2 

OR MORE 


TRA 

READ3 

GSW GREATER THAN 2 


TRA 

READ3 

GSW EQUAL TO 2 


TRA 

READ1 

GSW LESS THAN 2 

* THIS 

PART IS 

FOR SINGLE 

TAPE (A7) ONLY 

READ1 

TSX 

READA *4 

READ TAPE A7 


CALL 

RETURN 

MATMUL (=3 ) 
READT 

SET FOR DATB ONLY 

* THIS 

PART IS 

FOR NTH TIME OF DOUBLE TAPE (A7.A8) READ 

READ3 

CLA 

BUF1 

WAVE01 


FSB 

BUF4+12 

WAV1 


TNZ 

CLA 

*+3 

*1 

TEST A8 STATUS 


STO 

FSW1 

TURN A 7 FILL SWITCH TO ON 


CLA 

FSB 

BUF2 

BUF5+12 

WAVE02 


TNZ 

CLA 

*+3 

CONTINUE 


STO 

FSW2 

TURN A8 FILL SWITCH TO ON 


NZT 

FSW1 

IS A 7 BEING FILLED 


TRA 

READ31 

NO 


NZT 

FSW2 

IS A8 BEING FILLED 


TRA 

CALL 

READ41 
MESAGE ( -3 ) 

NO 

READ31 

NZT 

FSW2 

IS A8 BEING FILLED 


TRA 

READ42 

NO 

READ4 

TSX 

READA t 4 

READ A 7 ONLY 


CLA 

BUF 1 

WAVE01 


CAS 

BUF2 

COMPARE WAVE02 


TRA 

TRUBL 

WAVEOl.GT .WAVE02 


TRA 

ALINE 

WA VEO 1 • EQ • WAVE02 


TRA 

BNETH 

WAVE01.LT. WAVE02 

READ41 

TSX 

READB ♦ 4 

READ A8 ONLY 


CALL 

RETURN 

MATMUL ( =2 ) 
READT 

SET FOR DATBB ONLY 

READ42 

TSX 

TSX 

READA *4 
READB *4 

READ A 7 AND A8 


CLA 

BUF 1 

WAVE01 


CAS 

BUF2 

COMPARE WAVE02 


TRA 

EOFA 

WAVEOl.GT. WAVE02 


TRA 

READ43 

WAVE01.EQ.WAVE02 


TRA 

EOFA+2 

WAVE01.LT. WAVE02 

READ43 

CALL 

RETURN 

MATMUL ( =1 ) 
READT 

MULTIPLY DATB AND DATBB MATRICES 

ALINE 

STZ 

TRA 

FSW2 

READ43 

SET A8 FILL SWITCH TO 0FF=1 

BNETH 

CALL 

RETURN 

MATMUL ( *3) 
READT 

SET FOR DATB ONLY 


• THIS SUBROUTINE READS NTH RECORD FROM TAPE A7 


READA SX A KEEPA * A 

TSX .READ. 4 READ TABLES FROM 

PZE INP1 TAPE A 7 
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KEEPA 

* THIS 
READB 

KEEPB 

EOFA 

EOFB 

EOFC 

* THIS 

* THIS 
REAOF 

* THIS 


READR 


PZE 

EOFA 9 9 *— 2 


I ORT 

BUFl * .** 


AXT 

**.4 


TRA 

1*4 


SUBROUT I 

NE READS NTH RECORD FROM TAPE 

SX A 

KEEPB * A 


TSX 

• READ * A 

READ TABLES FROM 

PZE 

INP2 

TAPE B7 

PZE 

EOFB t * 2 


I ORT 

BUF2 • .** 


AXT 



TRA 

1 1 A 


ZET 

FSW2 

IS A 8 BEING FILLED 

TRA 

EOFC 

YES 

CLA 

= 1 

NO 

STO 

FSW1 

TURN A 7 FILL SWITCH TO ON =1 

TRA 

READ41+1 


ZET 

FSW1 

IS A 7 BEING FILLED 

TRA 

EOFC 

YES 

CLA 

= 1 

NO 

STO 

FSW2 

TURN A8 FILL SWITCH TO ON *1 

TRA 

BNETH 


CALL 

MESAGE ( =2 ) 

OUT OF FOCUS-STOP 

SUBROUT I 

NE READS. DIRECTS TAPE HANDLING OF A7 

SUBROUTINE READS *DI RECTS TAPE HANDLING OF A8 

SAVE 

1.2*3. 4.5.6 * 7 . 

I 

CLA 

GSW 


CAS 

*=2 

IS THERE MORE THAN 1 GAS 

TRA 

*+3 

YES 

TRA 

*+2 

YES 

TRA 

RE ADM 

NO 

PART IS 

FOR FIRST TIME 

OF DOUBLE TAPE (A7.A8) READ 

ZET 

FT SW 

IS THIS FIRST TIME 

TRA 

READG 


TSX 

• OPEN *4 


PZE 

INP1 


TSX 

• OPEN * 4 


PZE 

I NP2 


ST Z 

FSW1 

SET FILL SWITCHES 

STZ 

FSW2 

TO =* 0 

CLA 

= 1 


STO 

FT SW 

SET FIRST TIME SWITCH TO OFF=l 

TSX 

•READ. 4 


PZE 

INP1 


PZE 

EOF 1 . . *— 2 


IOCP 

BUF4 * * 20 

READtWl ( 11 ) .ETC... ) FOR GAS1 

I OR T 

BUF3* .** 

READ( I XA .DAT A ( 200 ) ) FOR GAS 1 

TSX 

• RE AD * 4 


PZE 

I NP2 


PZE 

EOF2. .*-2 


IOCP 

BUF5* .20 

READ ( W2 (11). ETC..#) FOR GAS2 

I ORT 

BUF2..** 

READt I XA ♦ DAT A ( 200 ) ) FOR GAS 2 

CALL 

PWD ( =3 ) 


TSX 

RE ADA . 4 


TSX 

READB *4 


CLA 

BUFl 

WAVE 01 

CAS 

BUF2 

COMPARE WAVE 02 

TRA 

MES1 

WAVE01.GT .WAVE02 

TRA 

MULTI 

WAVE01.EQ.WAVE02 

NOP 


WAVE01.LT. WAVE02 

CLA 

= 1 

LEAVE DATS AS IS 

STO 

FSW2 

SET FILL SWITCH FOR A8 TO ON = 1 ......... 

CALL 

MATMUL ( *3 ) 

SET FOR DATB ONLY 

RETURN 

READF 

/ 
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MULTI CALL 

MATMUL ( 3 1 ) 

MULTIPLY MATRICES DATB AND DATBB 

RETURN 

READF 


* THIS PART HANDLES THE SWITCHING OF TAPES WHEN A7 STARTS WITH WAVEOl.., 

* ADPATrD THAN 

LI AV/FTk? FRA, 


m vjK 1 t K i nr%n 

MES 1 TSX 

• REW * 4 

REWIND TAPE A7 

PZE 

INP1 


TSX 

•REWt 4 

REWIND TAPE A8 

PZE 

I NP2 

, 

CALL 

MESAGE ( =1 ) 

WRONG TAPE ON A7 

TRA 

READR 


READG NZT 

FSWl 

IS A 7 BEING FILLED 

TRA 

READH 

NO 

NZT 

FSW2 

IS A8 BEING FILLED 

TRA 

READJ 

NO 

CALL 

MESAGE (=3) 


* THIS PART IS 

FOR A 7 READ 

WHILE A8 IS BEING FILLED 

READH NZT 

FSW2 

IS A 8 BEING FILLED 

TRA 

READK 

NO 

READI TSX 

• READ * 4 

READ TAPE A7 ONLY 

PZE 

INP1 


PZE 

EOF3 * **— 2 


IOCP 

BUF4 * * 20 

W1 ETC#®# 

I ORT 

BUF3* t** 

READt I XA ♦ DATA < 200 ) ) FOR GAS 1 

CALL 

PWD ( = 3 ) 


TSX 

RE ADA 9 4 


CLA 

BUF 1 

WAVEOl 

CAS 

BUF2 

COMPARE WAVE02 

TRA 

TRUBL 

WAVEOl. GT. WAVE02 

TRA 

align 

WAVEOl. EQ.WAVE02 

TRA 

BELOW 

WAVE01.LT. WAVE02 

* THIS PART IS 

FOR TAPE A8 

READ ONLY IF A7 IS BEING FILLED.. 

READJ TSX 

• READ *4 

READ TAPE A8 ONLY 

PZE 

I NP2 


PZE 

E0F4* .*-2 


IOCP 

BUF5 * * 20 

W2 ETC... 

I ORT 

BUF2 * 

READ! IXA«DATA( 200) ) FOR GAS 2 

CALL 

PWD ( =3 ) 


TSX 

READS * 4 


CALL 

MATMUL (=2) 

SET FOR DAT88 ONLY 

RETURN 

READF 


• THIS PART READS A8 AFTER 

A 7 READS IN STEP ( CORRELATED 

READK TSX 

• READ 1 4 

TAPE A 7 

PZE 

INP1 


PZE 

E0F5 * * *— 2 


IOCP 

BUF4 ♦ * 20 

W1 *ETC# • • 

I ORT 

BUF3 * >** 

READ( I XA a DAT A ( 200 ) ) FOR GAS 1 

READL TSX 

• READ « 4 

TAPE A8 

PZE 

INP2 


PZE 

E0F6 ♦ **— 2 


IOCP 

BUF5 * * 20 

W2t ETC... 

I ORT 

BUF2 * 

READI I XA » DATA < 200 7 ) FOR GAS 2 

CALL 

PWD ( =3 ) 


TSX 

READA *4 


TSX 

RE ADB *4 


CLA 

BUF 1 

WAVEOl 

CAS 

BUF2 

COMPARE WAVE02 

TRA 

TRUBL 

WAVE01.GT •WAVE02 

TRA 

ALIGN+ 1 

• EQ. 

TRA 

TRUBL 

• LT# 

* THIS PART IS 

USED WHEN A 7 

READ FIRST MATCHES ••••••••.«• 

ALIGN STZ 

FSW2 

SET A8 FILL SWITCH TO OFF = 0 

CALL 

MATMUL ( =1 ) 

MULTIPLY MATRICES DATB AND DATBB 

RETURN 

READF 


BELOW CALL 

MATMUL ( “3 ) 

SET FOR DATB ONLY 

RETURN 

READF 



112 


laaamiui i 





TRUBL 

CALL 

MESAGE ( =2 ) 

OUT OF FOCUS - - STOP 

EOF 1 

CALL 

MESAGE ( =6 ) 

EXIT (FIRST TIME) 

EOF 2 

CALL 

MESAGE (=8) 

EXIT (FIRST TIME) 

EOF3 

TRA 

TRUBL 


EOF* 

TRA 

TRUBL 


EOF 5 

CLA 

= 1 



STO 

FSW1 

SET A 7 FILL SWITCH TO ON =1 


TRA 

READJ 


EOF 6 

CLA 

= 1 

NO 


STO 

FSW2 

SET A 8 FILL SWITCH TO ON =1 


CALL 

PWD < = 3 ) 



TSX 

READA** 



TRA 

BELOW 


* this 

PART IS 

FOR THE NTH 

TIME OF SINGLE TAPE READ 

READM 

ZET 

FTSW 



TRA 

READN 



CLA 

= 1 



STO 

FTSW 

SET FIRST TIME SWITCH TO OFF * 1 


TSX 

•OPENt* 



PZE 

INP1 


READN 

TSX 

.READ. A 



PZE 

INP1 



PZE 

EOFl * .*-2 



IOCP 

BUF* * ♦ 20 

W1 9 ETC* • • 


I OR T 

BUF3 * *** 

READ( IXA . DATA ( 200 ) ) FOR GAS 1 


CALL 

PWD ( = 1 ) 

PRINT WORDS (11) AND DATA (200) 


TSX 

READA** 



TRA 

BELOW 


FSW1 

BSS 

1 

FILL SWITCH FOR GAS 1 

FSW2 

BSS 

1 

FILL SWITCH FOR GAS 2 

ZA 

CONTRL 

ZA 



USE 

ZA 


BUF3 

BSS 

202 

WAVEO .DATA (200) * I ATM 


USE 

PREVIOUS 


ZB 1 

CONTRL 

ZB 1 



USE 

ZBI 

TABLE FOR WAV01 

BUF 1 

BSS 

BSS 

BSS 

201 

1 

1 

AND DATB ( 200 ) 

I X E 
A7M 

TAPE A 7 MULTIPLICITY SWITCH 

FTSW 

BSS 

1 

FIRST TIME SWITCH ( I N I T I ALLY=0 ) 


USE 

PREVIOUS 


ZB 2 

CONTRL 

ZB2 



USE 

ZB2 

TABLE FOR WAV02 

BUF2 

BSS 

201 

AND DATBB (200) 

I X F 

BSS 

1 


A8M 

BSS 

1 

TAPE A8 MULTIPLICITY SWITCH 

GSW 

BSS 

1 

GAS MULTIPLICITY SWITCH 


USE 

PREVIOUS 


ZC1 

CONTRL 

ZC1 



USE 

ZC1 

(W1 ( 11 ) . WAV1 .WAVF1 ♦ INTI » IGAS1 . 

BUF* 

BSS 

20 

(HT31 .ANGL 1 . NUM1 .JTAN1 .LAI 


USE 

PREVIOUS 


ZC2 

CONTRL 

ZC2 



USE 

ZC2 

W2 l 1 1 ) .WAV2.WAVF2.INT2.1GAS2. 

BUF5 

BSS 

20 

HT32 » ANGL2 » NUM2 .JTAN2.LA2 


USE 

PREVIOUS 



END 




SIBMAP 

INTA 

1 00 * NODECK * NOL I ST * NOREF 

INTA 

SAVE 



SX A 

SAV* ♦ * 


CLA 

IXA 


ALS 

17 


250 CARDS 
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STD 

MOD+2 



CLA 

4 * A 



STA 

TEM 



CLA* 

3*4 


MOD 

TSX 

T I N1 * 4 



PZE 

DATA 



PZE 

3 * *** 


SAVA 

AXT 

** * A 



STO* 

TEM 



RETURN 

INTA 


TEM 

BSS 

1 


INTB 

SAVE 




SX A 

SAV4*4 



CLA 

IXB 



ALS 

17 



STD 

MOD1+2 



CLA 

A t A 



STA 

TEM1 



CLA* 

3*4 


MODI 

TSX 

T I N1 * 4 



PZE 

• DATS 



PZE 

1* *** 



LX A 

SAVA*A 



STO* 

TEM1 



RETURN 

INTB 


TEM1 

BSS 

1 


BB 

CONTRL 

BB 



USE 

BB 


BUF6 

BSS 

19 


I X B 

BSS 

1 


I WAVE 

BSS 

1 


WAVEO 

BSS 

1 



USE 

PREVIOUS 


ZA 

CONTRL 

ZA 



USE 

ZA 


I X A 

BSS 

1 


DATA 

BSS 

200 


I ATM 

BSS 

1 



USE 

PREVIOUS 


ZB 1 

CONTRL 

ZB 1 



USE 

ZB 1 

TABLE FOR WAV01 

WA VO 1 

BSS 

1 


DATS 

BSS 

200 

AND DATS! 200 ) 

IXE 

BSS 

1 


A7M 

8SS 

1 

TAPE A 7 MULTIPLIC 

FTSW 

BSS 

1 

FIRST TIME SWITCH 


USE 

PREVIOUS 



END 




SIBMAP 

INTRPC 

50 * NODECK * NOL I ST ♦ NOR EF 

I NT RPC 

SAVE 

(1*2*4) I 


CLA 

A * A 


STA 

CSF 


CLA* 

5. A 


ALS 

18 


STD 

PAIRS 


CLA* 

3*4 


STO 

CWAVE 


TSX 

TIN1* A 


PZE 

TAB15 

PAIRS 

PZE 

1*0*** 


STO* 

CSF 


RETURN 

INTRPC 


I 


TY SWITCH 
( INI TI ALLY=0 ) 


53 CARDS 
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CWAVE 

PZE 

0 

CSF 

PZE 

0 

* 

TABLE OF 

PHI VS. WAVE NUMBER 

DATAC 

CONTRL 

DATAC 


USE 

DATAC 

TAB15 

BSS 

100 


USE 

PREVIOUS 

I NTRPD 

SAVE 

(1*2*4*3*5*7>I 


CLA 

4*4 


STA 

EPSI 


CLA* 

5*4 


ALS 

18 


STD 

PARES 


CLA* 

3*4 


TSX 

T I N1 *4 


PZE 

TABLE 

PARES 

PZE 

1*0*** 


STO* 

EPSI 


RETURN 

I NTRPD 

EPSI 

PZE 

0 

* TABLE OF EPSILON VS. WAVE NUMBER 

DATA D 

CONTRL 

DA TAD 


USE 

DA TAD 

TABLE 

BSS 

100 


USE 

PREVIOUS 


END 



SIBMAP 

TIN1 

125 *M9A * NODECK *NOL I ST *NOREF 

* 

INTERPOLATION SUBROUTINE TIN1 

* 

CARDS 

COLUMN 

* 

FAP 



ENTRY 

TIN1 

TIN1 

STO 

TIN1+98 


SXD 

T I N 1 + 87 * 1 


SXD 

T I Nl + 88 * 2 


SXD 

T I N 1 + 89 *4 


CLA 

2*4 


STO 

TIN1+94 


ADD 

T I N 1+9 1 


PAX 

0*1 


SXD 

T I N 1+29 * 1 


ALS 

1 


STA 

T I N 1+95 


CLA 

T I N 1+94 


ARS 

17 


SUB 

T I N 1+95 


PAX 

0*1 


SXD 

TIN1+31.1 


ADD 

T I N 1+95 


PAX 

0*1 


ADD 

1*4 


STA 

T I N 1 + 2 3 


STA 

T I N 1+4 1 


STA 

T I Nl+43 


CLA 

TIN1+98 


CAS 

0*1 


T IX 

*-1*1 *2 


TRA 

TIN1+26 


CLA 

TIN1+94 


LBT 



TRA 

T I Nl+41 


T I X 

T I Nl+31 *1*0 


40 CARDS 
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LXD 

T I Nl+90 * 1 

TXL 

T I Nl+33 *1*0 

LXD 

T I Nl+3 1*1 

PXD 

0*1 

ARS 

18 

CHS 


ADD 

TIN1+23 

STA 

T I N 1+54 

ADD 

TIN1+91 

STA 

TIN1+52 

TRA 

TIN1+50 

CLA 

0*1 

TXI 

T I N 1+43 * 1 * 2 

FAD 

0*1 

LRS 

35 

FMP 

T I N 1+92 

CAS 

TIN1+98 

TXI 

T I Nl + 29 * 1 *- 1 

TXI 

T I Nl + 29 * 1 *- 1 

TXI 

T I Nl+29 * 1 * - 3 

LXD 

T I Nl + 90 * 2 

LX A 

T I N 1 + 9 5 * 1 

CLA 

0*1 

STO 

COM+ 1*2 

CLA 

0*1 

FSB 

T I N 1+98 

STO 

COM ♦ 2 

TXI 

TINl+58*2*-2 

T IX 

TIN1+52*1*2 

CLA 

COM 

STO 

COM *2 

LXA 

T I N 1+9 A * 4 

LXD 

T I N 1 + 9 3 * 2 

TXI 

T I N 1+64 * 2 * -2 

PXD 

0*2 

PDX 

0*1 

CLA 

COM+2 * 1 

FSB 

COM *2 

TZE 

T I N 1+79 

STO 

T I N 1+97 

LDQ 

COM *2 

FMP 

COM+3 * 1 

STO 

TIN1+96 

LDQ 

COM+2 *1 

FMP 

COM+ 1*2 

FSB 

T I N 1+96 

FDP 

T I N 1+97 

STQ 

COM+3 *1 

TXI 

T I N 1 + 66 ♦ 1 • ~2 

CLA 

COM+2 *2 

STO 

COM+2 *1 

T IX 

T I N 1+63 *4*1 

CLA 

COM+1 *1 

LXD 

TIN1+87 • 1 

LXD 

T I N 1 + 88 * 2 

LXD 

T I N 1 + 89 * 4 

TRA 

3*4 

HTR 


HTR 


HTR 


HTR 


HTR 

1*0*0 

DEC 

.5 

HTR 

0*0*2 
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U1 u>o 


COM 


AO 


HTR 
HTR 
HTR 
HTR 
HTR 
PZE 
PZE 
BSS 
END 

108 CARDS 


SIBMAP RITEB7 1 00 . L I ST *XR7 ♦ M9A . NODECK 

* RITEB7 WRITES BUFFER TRAN < 1 1 ) ON TAPE B7 . THIS INCLUDES TABLES OF 

* TRANSMISSIONS AND OTHER PARAMETERS FOR INPUT TO ATMOS3 * 

OUT FILE * B ( 3 ) » B ( 3 ) *BLK=3Q0. BIN. HOLD * MOUNT * OUTPUT 

RITEB7 SAVE < A *2 ♦ 1 *3 *5 . 6 * 7 ) I 

CLA* 3. A 

ALS 18 

STD II 

TSX .OPEN. A 

PZE OUT 

TSX .WRITE . A 

PZE OUT 

II I OCD TRAM . * 300 

RETURN RITEB7 

TS CONTRL TS 

USE TS 

TRAM BSS 300 ALL PARAMETERS FOR OUTPUT 

USE PREVIOUS 

RITE1 SAVE 1 .2.3.A.5.6 *7. I 

CLA* 3 * A 

ADD *1 ADD 1 FOR IXA ITSELF 

ALS 18 

STD JJ 

TSX .OPEN . A 

PZE OUT 

TSX .WRITE.A 

PZE OUT 

IOCP BUF 1 ♦ t 20 

JJ I OCD BUF2 • • 2 00 

RETURN RITE1 

ZC1 CONTRL ZC1 

USE ZC1 

BUF 1 BSS 20 WORDS ( 11 > .WWAVEI .WAVEF. INTER. IGAS. LA. 

USE PREVIOUS HT3. ANGLE. NUM.JTAN 

ZA CONTRL ZA 

USE ZA 

BUF 2 BSS 200 I X A . DAT A < 2 00 ) 

USE PREVIOUS 

READ2 SAVE 1 . 2 .3 . A . 5 . 6 .7 . I 

TSX .CLOSE. A 

PZE OUT 

RETURN READ2 

END 

i 

2 2 11 

SPECTRAL RESPONSE =UN I T Y 


§ 88 : 81:8 

00 . 001.0 

3000.01.0 

2205.00305.00105.00305.00 


EOF E SY 


55 CARDS 
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"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof” 

—National Aeronautics and Space Act op 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 
important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless of 
importance as a contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: Information receiving limited distribu- 
tion because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Scientific and technical information generated 

under a NASA contract or grant and considered an important contribution to 
existing knowledge. 

TECHNICAL TRANSLATIONS: Information published in a foreign 
language considered to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information derived from or of value to NASA 
activities. Publications include conference proceedings, monographs, data 
compilations, handbooks, sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION PUBLICATIONS: Information on tech- 
nology used by NASA that may be of particular interest in commercial and other 
non-aerospace applications. Publications include Tech Briefs, Technology 
Utilization Reports and Notes, and Technology Surveys. 


Details on the availability of these publications may be obtained from: 


SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



